Resumos
A plasmônica é uma disciplina emergente que investiga a interação entre a luz e os elétrons quase livres em sistemas metálicos em uma interface metal/dielétrico. Esta interação dá origem aos denominados Plásmon-Poláritons de Superfície SPPs (da sigla em inglês Surface Plasmon-Polaritons), que são um tipo de quasipartícula resultante do acoplamento entre o fóton e o plásmon de superfície. A relação de dispersão dos SPPs é um conceito chave no entendimento da interação da luz com a matéria em estruturas nanométricas, com aplicações diversas nos campos da fotônica e da bioengenharia. Neste trabalho apresentamos em detalhes a dedução da relação de dispersão para o plásmon-poláriton de superfície para uma interface plana metal/dielétrico. Para as partes real e imaginária do vetor de onda ao longo da interface, deduzimos expressões gerais cujos casos particulares reproduzem resultados conhecidos na literatura. Em seguida é realizada uma discussão sobre o significado físico dos resultados obtidos. Por fim, utilizando o modelo de Drude, apresentamos uma expressão explícita para a relação de dispersão na forma frequência versus vetor de onda.
Palavras-chave:
plasmônica; plásmon de superfície; plásmon-poláriton de superfície
Plasmonics is an emergent area of the nanophotonics that studies interaction between electromagnetic fields and free electrons in a metal on a metal/dielectric interface. Due to this interaction, the called surface plasmon polaritons (SPP) are generated. SPPs are electromagnetic excitations coupled to the collective oscillation of conduction electrons, which propagate along the interface between a conductor and a dielectric medium, with a broad range of applications, from photonics to bioengineering. One of the basic concepts related to SPPs is the dispersion relation, basic for understanding of the interaction between light and metallic nanostructures. In this paper we derive the dispersion relation of surface plasmon-polariton at planar interface and its fundamental properties are outlined.
Keywords:
plasmonic; surface plasmon; surface plasmon polariton
1. Introdução
A Física da Matéria Condensada, por sua grande amplitude de temas investigados, é considerada a mais ampla área da Física [1[1] Marvin L. Cohen, Phys. Rev. Lett. 101, 250001 (2008).]. De um ponto de vista fundamental, os sistemas pesquisados em matéria condensada constituem problemas de muitos corpos de difícil solução, o que requer a utilização conjunta dos conceitos básicos de mecânica quântica, física estatística e eletromagnetismo. De um ponto de vista prático, muitos dos avanços tecnológicos que fazem parte do nosso cotidiano tiveram origem nos estudos em pesquisa básica nessa área. Podemos citar como exemplos a invenção do transistor, do laser de estado sólido, aparelhos de DVD, telas de cristal líquido, dentre muitos outros. Seus impactos não se limitam à engenharia de materiais, mas são responsáveis também pela inserção de novas tecnologias nas ciências da vida, como na regeneração de tecidos, terapias gênicas e antitumorais, e diversos outros avanços em bioengenharia [2[2] Ado Jorio, Bioengineering Applications of Carbon Nanostructures (Springer Verlag, Berlin, 2016).].
A intensa investigação, tanto teórica quanto experimental, das propriedades eletrônicas, óticas, magnéticas, etc, que ocorrem na Física da Matéria Condensada possibilitou o desenvolvimento de novas disciplinas como a fotônica e a spintrônica. A plasmônica é uma subárea da nanofotônica que estuda a interação entre a radiação eletromagnética e os elétrons quase livres (elétrons de condução) em sistemas metálicos em uma interface metal/dielétrico. De forma mais pormenorizada, podemos dizer que a plasmônica lida com a geração, propagação e detecção de ondas de plasma ou ondas plasmônica em sistemas metálicos [3[3] M. Dragoman and D. Dragoman, Progress in Quantum Electronic 32, 1 (2008).]. Vejamos o que exatamente significa essa afirmação.
Muitas propriedades físicas dos metais e semicondutores dopados podem ser entendidas em um contexto teórico no qual o comportamento dos elétrons de condução desses materiais se assemelha à de um plasma. Um plasma é um conjunto eletricamente neutro de partículas carregadas que interagem via forças coulombianas, no qual pelo menos um tipo é móvel e que exibe um comportamento coletivo [4[4] Ricardo S. Oyarzabal, José D. Szezech Junior, Antonio M. Batista, Iberê L. Caldas, Ricardo L. Viana e Kelly C. Iarosz, Revista Brasileira de Ensino de Física 37, 2308 (2015).]. Por analogia foi cunhada a expressão Plasma em Estado Sólido para descrever o comportamento dos elétrons de condução nesses materiais [5[5] P.M. Platzman and P.A. Wolff, in: Solid State Physics, Suppl. 13, edited by F. Seitz, D. Turnbull and H. Ehrenrich (Academic Press, New York, 1973).]. A neutralidade elétrica do sistema é garantida pela carga positiva da rede, cujo arranjo periódico discreto pode ser ignorado em uma primeira aproximação e tomada como contínua.
Esses elétrons apresentam uma oscilação coletiva cuja frequência, denominada frequência de plasma, é dada por (no sistema cgs) [6[6] Lukas Novotny and Bert Hecht, Principles of Nano-Optics (Cambridge University Press, Cambridge, 1999), 2ª ed, p. 5, p. 371.]
onde é a carga eletrônica, é a massa efetiva do elétron a qual incorpora a interação elétron-rede e é a concentração dos elétrons de condução.
Segundo a teoria quântica de campos é possível quantizar um campo oscilatório e associar ao mesmo um quantum elementar ao qual se pode atribuir propriedades de partícula [7[7] Rodrigo Gonçalves Pereira e Eduardo Miranda, Revista Brasileira de Ensino de Física 24, 237 (2002).]. Por exemplo, os fótons são as partículas (quantum elementar) associadas ao campo eletromagnético (uma interessante discussão sobre o conceito de fóton é encontrada na Ref. [8[8] Indianara Silva, Revista Brasileira de Ensino de Física 37, 4204 (2015).]). O quantum associado às oscilações coletivas de um gás de elétrons é denominado de plásmon. O plásmon é um exemplo das denominadas quasipartículas (ou excitações elementares) as quais, de forma simplificada, podem ser definidas como entidades fictícias que descrevem de modo relativamente simples o comportamento de sistemas de muitos corpos em matéria condensada [9[9] Richard D. Mattuck, A Guide to Feynman Diagrams in the Many-body Problem (McGraw-Hill, New York, 1967).,10[10] S. Nakajima, Y. Toyozawa and R. Abe, The Physics of Elementary Excitations (Springer-Verlag, Berlin, 1980).]. A energia de um plásmon é calculada pela relação
onde é a constante de Planck normalizada. Para metais, a concentração dos elétrons de condução é , o que implica em uma energia do plásmon da ordem de . Para um semicondutor, de grande volume (bulk), considerando , teremos da ordem de (este valor pode mudar significativamente em materiais de baixa dimensionalidade). Classicamente, o plásmon corresponde a uma oscilação na densidade de carga eletrônica e esta oscilação pode ocorrer no interior ou na superfície do material. No primeiro caso falamos em plásmons de volume (bulk plasmons) ou simplesmente plásmons. No segundo caso falamos em plásmons de superfície (surface plasmons).
Outro exemplo de quasipartícula que nos interessa neste trabalho é o poláriton, o qual resulta do acoplamento entre o fóton e uma excitação elementar do sólido. Há vários tipos de poláritons, a depender da quasiparticula que se acopla com o fóton. Fônon-poláriton que resulta do acoplamento entre o fônon transversal ótico e o fóton; o éxciton-poláriton que resulta do acoplamento entre o éxciton (par elétron-buraco ligado) e o fóton e, finalmente, o plásmon-poláriton de volume - BPP (da acrossemia Bulk Plasmon Polariton) que resulta da interação entre os fótons e os elétrons livres no material. No presente trabalho estamos interessados neste último. Vale a pena observar que a denominação plásmon-poláriton não é muito adequada pois pode induzir o leitor a imaginar que exista um acoplamento entre o plásmon e o fóton. Conforme é discutido no Apêndice A Apêndice A Como vimos na Introdução, as variadas propriedades físicas de um sólido podem ser descritas em termos do conceito fundamental de quasipartículas ou excitações elementares. A interação entre diferentes tipos de excitações elementares dá origem a novas quasipartículas que resultam do acoplamento entre estas (os modos híbridos ou modos acoplados). Como exemplo típico podemos citar os diferentes tipos de poláritons. Contudo, é importante enfatizar que só pode haver acoplamento entre modos de vibração paralelos, i.e., para haver acoplamento ambos os modos devem ser longitudinais ou ambos devem ser transversais. Assim, o fônon-poláriton resulta do acoplamento entre o fóton e o fônon transversal ótico. Não é possível, por exemplo, haver acoplamento entre o fóton e o fônon longitudinal ótico. As excitações elementares podem, em grande parte, ser entendidas a partir da análise das equações de Maxwell em um meio material, escrita em termos de componentes de Fourier (ver Apêndice B). A lei de Gauss, equação (B.16) pode ser escrita como (A.1) ϵ ( k → , ω ) k → ⋅ E → ( k → , ω ) = 0 . Da equação (A.1) temos duas soluções possíveis: (A.2) ϵ ( k → , ω ) = 0 e (A.3) k → ⋅ E → ( k → , ω ) = 0 . A primeira solução, equação (A.2), está associada as oscilações longitudinais coletivas do plasma que, na linguagem de segunda quantização, constituem os plásmons que são determinados pelas raízes de ϵ(k→,ω). Usando o modelo de Drude e desprezando amortecimento (Γ=0), vemos da equação (51) que a condição ϵ(k→,ω)=0 ocorre para ω=ωpl/ϵ. Como os plásmons de volume correspondem a modos longitudinais, não pode haver acoplamento entre estes e os fótons. A segunda solução, equação (A.3), define a condição de transversalidade. Das equações de Maxwell, na ausência de fontes, obtemos a seguinte equação (A.4) ∇ × ( ∇ × E → ) = - 1 c 2 ∂ 2 ∂ t 2 ( ϵ E → ) , a qual, em termos de componentes de Fourier, pode ser escrita como (A.5) ( k → ⋅ E → ( k → , ω ) ) k → - k 2 E → ( k → , ω ) = - ω 2 c 2 ϵ ( k → , ω ) E → ( k → , ω ) . Utilizando (A.3) em (A.5) obtemos (A.6) k 2 = ω 2 c 2 ϵ ( k → , ω ) , que corresponde à relação de dispersão para o fóton se propagando em um meio material. Vemos que essa propagação é mediada pela interação entre o fóton e os elétrons de condução do material. O fóton revestido por essa interação dá origem a uma nova quasipartícula denominada plásmon-poláriton de volume. , não pode haver acoplamento entre o fóton e o plásmon de volume já que este representa uma oscilação longitudinal, enquanto o fóton repesenta uma oscilação transversal. Contudo, em uma interface em que incida uma onda eletromagnética p-polarizada, a polarização pode ser decomposta em uma componente perpendicular à interface e outra componente paralela à interface. A componente paralela torna-se uma onda longitudinal e, desse modo, pode acoplar-se com o plásmon de superfície. Esse acoplamento dá origem ao Plásmon-Poláriton de Superfície (SPP da acrossemia Surface Plasmon Polaritons).
Do ponto de vista da física clássica, os Plásmon-Poláritons de Superfície podem ser visualizados como uma onda eletromagnética superficial que se propaga na interface entre o metal e o dielétrico e que decai exponencialmente com a distância na direção perpendicular à interface, provocando o confinamento da luz na interface metal/dielétrico. Esse efeito de confinamento pode ser utilizado em uma grande variedades de aplicações tais como na fabricação de LEDs mais eficientes, no aumento da resolução de microscópios e na melhora da sensibilidade de detectores químicos e biológicos [11[11] Harry A. Atwater, Scientific American 296, 56 (2007).]. Uma das mais interessantes aplicações da plasmônica ocorre na denominada microscopia de campo próximo. Com o advento do Microscópio Ótico de Campo Próximo SNOM (da sigla em inglês Scanning Near-Field Optical Microscope) é possivel obter imagens com uma resolução muito maior do que a estabelecida pelo critério de Rayleigh [12[12] Francisco E.G. Guimarães, Revista Brasileira de Ensino de Física 18, 56 (1996).], o qual, devido ao fenômeno de difração da luz [13[13] Daniel M. Reis, Edson M. Santos e A.V. Andrade-Neto, Revista Brasileira de Ensino de Física 37, 2312 (2015).], determina o limite da resolução espacial de um sistema ótico convencional a aproximadamente metade do comprimento de onda da luz utilizada como prova. A grande vantagem do SNOM é que ele combina a grande resolução dos Microscópios de Varredura por Sonda SPM (da sigla em inglês Scanning Probe Microscope) com a versatilidade das técnicas óticas [6[6] Lukas Novotny and Bert Hecht, Principles of Nano-Optics (Cambridge University Press, Cambridge, 1999), 2ª ed, p. 5, p. 371.,14[14] Thiago L. Vasconcelos, Bráulio S. Archanjo, Benjamin Fragneaud, Bruno S. Oliveira, Juha Riikonen, Changfeng Li, Douglas S. Ribeiro, Cassiano Rabelo, Wagner N. Rodrigues, Ado Jorio, et. al., ACSNano 9, 6297 (2015).].
O estudo dos Plásmon-Poláritons de Superfície fornece importante informações sobre o comportamento da luz em superfícies metálicas de tamanhos nanométricos. Esse comportamento pode ser investigado a partir da análise das equações de Maxwell e suas apropriadas condições de contorno. Este é o caminho que seguiremos neste trabalho. Em particular, deduziremos a relação de dispersão para os Plásmon-Poláritons de Superfície em uma interface plana metal/dielétrico e exploraremos as informações físicas que podem ser obtidas a partir dessa relação.
2. Dedução da relação de dispersão dos plásmon-poláritons de superfície em uma interface metal-dielétrico plana
As propriedades dos plásmon-poláritons de superfície na interface entre um metal e um dielétrico podem ser obtidas a partir das equações de Maxwell em cada meio e das condições de contorno associadas ao campo eletromagnético na superfície de descontinuidade entre os dois meios.
As condições de contorno sobre os vetores do campo eletromagnético nos informam como esses se comportam na interface entre dois meios com propriedades materiais diferentes. Este é um problema abordado nos livros textos de eletromagnetismo [ver, por exemplo, Ref. [15[15] John R. Reitz, Frederick J. Milford e Robert W. Christy, Fundamentos da Teoria Eletromagnética. (Editora Campus, São Paulo, 1982), p. 95, p. 199.]]. Mostraremos a seguir que, partindo das soluções das equações de Maxwell e das condições de contorno, podemos obter a relação de dispersão para o SPP em uma interface plana metal dielétrico e dessa relação inferir as propriedades físicas do sistema.
Vamos considerar dois meios semi-infinitos separados por uma interface plana. O meio é um metal caracterizado por uma função dielétrica complexa dependente da frequência , onde e representam as partes real e imaginária, respectivamente. Por sua vez, para o meio , meio dielétrico, assumiremos que a função dielétrica é uma grandeza real e constante. Vamos escolher um sistema de coordenadas cartesianas onde é o plano de incidência e a interface coincide com o plano , conforme esquematizado na Figura 1.
Esquema geométrico para uma interface metal/dielétrico no plano . O meio dielétrico ocupa o semi-espaço enquanto o metal o semi-espaço .
Os campos de uma onda plana monocromática podem ser representados como superposições de duas polarizações lineares perpendiculares que, por conveniência, são escolhidas como perpendicular ao plano de incidência (denominado s-polarizado) e outra paralela ao plano de incidência (denominado p-polarizado). É possível mostrar [16[16] Stefan Alexander Maier, Plasmonics: Fundamentals and Applications (Springer, New York, 2007), p. 10-27.] que a componente perpendicular (onda s-polarizada) não produz carga superficial. Assim, podemos considerar apenas ondas p-polarizada, cujas expressões para os campos elétrico e magnético podem ser escritas como [17[17] H. Raether, Surface Plasmons on Smooth and Rough Surfaces and on Gratings (Springer-Verlag, Berlin, 1998), Appendix I.]
para
para , onde é a componente da amplitude do campo elétrico, é a componente do vetor de onda e representa o vetor unitário na direção do eixo . Idem para as outras componentes. As equações (3-6) podem ser escritas de forma mais compacta e mais adequada para manipulação como
onde o índice identifica o meio: indica o meio metálico e o meio dielétrico com
e
com o sinal para e o sinal para .
As equações de Maxwell em um meio material não-magnético, na ausência de fontes, em termos de componentes de Fourier são escritas como (ver detalhes no Apêndice B Apêndice B As equações de Maxwell em um meio material, utilizando o sistema cgs, são escritas como (B.1) ∇ × H → ( r → , t ) = 1 c ∂ D → ( r → , t ) ∂ t + 4 π c J → c o n ( r → , t ) (B.2) ∇ × E → ( r → , t ) = - 1 c ∂ B → ( r → , t ) ∂ t (B.3) ∇ ⋅ D → ( r → , t ) = 4 π ρ e x t ( r → , t ) (B.4) ∇ ⋅ B → ( r → , t ) = 0 , onde D→(r→,t) é o vetor deslocamento elétrico, J→con(r→,t) é a densidade de corrente de condução e ρext(r→,t) é a densidade volumétrica de carga. Muitas vezes é conveniente escrever as equações de Maxwell em termos de componentes de Fourier. Assumindo um meio sem dispersão espacial podemos escrever o campo E→(r→,t) na seguinte forma (B.5) E → ( r → , t ) = 1 2 π ∫ - ∞ ∞ E → ( r → , ω ) e - ı ω t d ω , onde E→(r→,ω) é a transformada de Fourier de E→(r→,t) dada por (B.6) E → ( r → , ω ) = ∫ - ∞ ∞ E → ( r → , ω ) e ı ω t d t , com expresões semelhantes para os outros campos. Diferenciando o campo elétrico temos que (B.7) ∇ × E → ( r → , t ) = 1 2 π ∫ - ∞ ∞ ∇ × E → ( r → , ω ) e - ı ω t d ω (B.8) ∇ ⋅ E → ( r → , t ) = 1 2 π ∫ - ∞ ∞ ∇ ⋅ E → ( r → , ω ) e - ı ω t d ω (B.9) ∂ E → ( r → , t ) ∂ t = 1 2 π ∫ - ∞ ∞ ( - ı ω ) E → ( r → , ω ) e - ı ω t d ω Usando as equações acima nas equações de Maxwell, estas ficam (B.10) ∇ × H → ( r → , ω ) = - ı ω c D → ( r → , ω ) + 4 π c J → c o n ( r → , ω ) (B.11) ∇ × E → ( r → , ω ) = ı ω c B → ( r → , ω ) (B.12) ∇ ⋅ D → ( r → , ω ) = 4 π ρ e x t ( r → , ω ) (B.13) ∇ ⋅ B → ( r → , ω ) = 0 Para um meio não-magnético (B→=H→) e na ausência de fontes temos que (B.14) ∇ × H → ( r → , ω ) = - ı ω c ϵ ( ω ) E → ( r → , ω ) (B.15) ∇ × E → ( r → , ω ) = ı ω c H → ( r → , ω ) (B.16) ∇ ⋅ ϵ ( ω ) E → ( r → , ω ) = 0 (B.17) ∇ ⋅ H → ( r → , ω ) = 0 onde usamos que D→(r→,ω)=ϵ(ω)E→(r→,ω). )
onde, mais uma vez, o índice ou identifica o meio.
Devemos agora indicar como os vetores de campo variam na interface entre os dois meios. A componente tangencial do campo elétrico é contínua na interface entre dois meios, bem como a componente normal do vetor deslocamento se não há cargas superficiais livres, como acontece no presente caso. Do mesmo modo, a componente normal do campo magnético é contínua através da interface [15[15] John R. Reitz, Frederick J. Milford e Robert W. Christy, Fundamentos da Teoria Eletromagnética. (Editora Campus, São Paulo, 1982), p. 95, p. 199.]. Utilizando essas condições de contorno obtemos
Utilizando a equação (17) vemos que a componente do vetor de onda é contínua, i.e.,
Vamos utilizar as equações de Maxwell e as condições de contorno para deduzir a relação de dispersão para o plásmon-poláriton de superfície. Aplicando o rotacional na equação (8) temos que
Utilizando as equações (7) e (21) em (13) temos que
Utilizando as equações (9-11) em (22) obtemos explicitamente para os meios 1 e 2:
Para o meio
Para o meio
Das equações (17), (24) e (26) obtemos
a qual, junto com a Eq. (18), forma o sistema
Para que o sistema acima tenha solução não trivial o determinante tem que ser igual a zero. Assim, obtemos
De forma semelhante, partindo da equação (14) obtemos após alguns cálculos
Das equações (29) e (30) obtemos a relação de dispersão
A equação (31) expressa a relação entre o vetor de onda ao longo da interface e a frequência angular e é denominada relação de dispersão para os plásmon-poláritons de superfície em uma interface entre um metal e um dielétrico. Essa equação é fundamental para o entendimento da interação entre a luz e a matéria em uma interface metal-dielétrico, que dá origem ao plásmon-poláriton de superfície. Suas implicações físicas serão discutidas nas próximas Seções.
Deve ser observado que a expressão (31) foi deduzida utilizando apenas as equações de Maxwell e, portanto, ela é geral. Como é uma função complexa, é também uma função complexa, i.e., , cujas expressões gerais serão deduzidas a seguir.
3. Vetor de onda dos plásmon-poláritons de superfície: caso geral
Considerando que a equação (31) pode ser reescrita como
onde
Da equação (32) obtemos que
Para os metais nobres em geral, temos que a parte real da função dielétrica é negativa e a parte imaginária é positiva, i.e., e . Para os dielétricos . Nessas condições vemos da Eq. (34) que é sempre uma grandeza real positiva. Por sua vez, pode ser positiva, nula ou negativa. Se , . Se , isso implica em . Finalmente, se , vemos que . Contudo, qualquer que seja o caso, e são números reais positivos.
As expressões gerais (35) e (36) são suficientementes complicadas para que valha a pena analisar alguns casos limites em que elas se reduzem a formas mais simples conhecidas da literatura. Em geral é considerada a situação em que a parte imaginária da função dielétrica do metal é muito menor que o módulo da parte real, i.e., . Neste limite vemos que
Então obtemos
Do mesmo modo obtemos da equação (36)
As equações (38) e (39) são as expressões normalmente encontradas na literatura [6[6] Lukas Novotny and Bert Hecht, Principles of Nano-Optics (Cambridge University Press, Cambridge, 1999), 2ª ed, p. 5, p. 371.].
4. Significado físico dos resultados obtidos
Vamos agora iniciar uma discussão sobre o significado físico dos resultados obtidos nas seções anteriores.
4.1. Comprimento de onda e comprimento de propagação dos plásmon-poláritons de superfície
Utilizando (20), vemos das expressões (3) e (5) que o campo elétrico na interface é da forma
onde e são dados, respectivamente, pelas equações (35) e (36) no caso geral, ou (38) e (39) no caso particular . A equação (40) representa um campo elétrico progressivo amortecido.
A parte real do vetor de onda, , determina o comprimento de onda do SPP, , o qual no caso limite é dado por
onde é o comprimento de onda da luz no vácuo. Nas regiões do visível e infravermelho, para metais nobres e um dielétrico típico, temos que . Como pode ser visto da equação (41), isso implica que
Vemos então que é sempre menor que o comprimento de onda da radiação que excita o sistema. Esse fato é a base que possibilita a construção de dispositivos óticos (circuitos plasmônicos) em tamanhos menores que o comprimento de onda da luz visível.
Já a parte imaginária, , determina o comprimento de propagação do SPP quando este se propaga ao longo da interface. Vemos da equação (40) que produz um decaimento exponencial da amplitude do campo elétrico na interface. Isso nos permite definir um comprimento de propagação para o SPP, que é a distância para a qual o campo elétrico decai do seu valor inicial, ou seja, (ou para a intensidade). Usando a equação (39) obtemos
O comprimento de propagação define um comprimento máximo para um circuito plasmônico. A Figura 2 mostra o comprimento de propagação do SPP em função do comprimento de onda, calculado pela equação (43), para uma interface ar/prata. Para calcular a função dielétrica utilizamos o modelo de Drude, ver equações (52) e (53) adiante, com os parâmetros e . Vemos que o SPP pode percorrer uma grande distância (considerando dispositivos na escala nanométrica) antes de se dispersar.
Comprimento de propagação para uma interface ar/prata, usando o modelo de Drude para o cálculo da função dielétrica com os parâmetros e .
4.2. Profundidade de penetração dos plásmon-poláritons de superfície
Na direção perpendicular à interface, o campo elétrico decai exponencialmente com a distância a interface nos dois meios. De fato, se substituimos a equação (31) na equação (30) obtemos
Como para metais típicos , vemos que é uma grandeza imaginária. Para o caso limite em que obtemos para o meio
onde
Para o meio temos que
onde
Utilizando as equações (45) e (47), as equações (5) e (3) podem ser rescritas como
que representa uma onda que se propaga na direção e decai exponencialmente na direção . Esse tipo de onda é denominada onda evanescente.
A profundidade de penetração é definida como a distância para a qual o campo eletromagnético se reduz a do seu valor na interface. Assim
Para uma interface ar/prata, para comprimentos de onda no intervalo entre a , a profundidade de penetração no ar, , está situada na faixa entre a . Por sua vez, para o mesmo intervalo de comprimentos de onda, a profundidade de penetração na prata será , sendo praticamente independente do comprimento de onda [3[3] M. Dragoman and D. Dragoman, Progress in Quantum Electronic 32, 1 (2008).]. Isso mostra que a profundidade de penetração no meio metálico é muito menor que no meio dielétrico.
A Figura 3 mostra de forma esquemática a oscilação de carga elétrica na interface metal/dielétrico. Deve ser observado que as cargas positivas indicam região na qual há uma diminuição de carga negativa e não um excesso de carga positiva. Também são representadas as linhas de força do campo elétrico nas proximidades da interface (linhas curvas com setas) bem como a profundidade de penetração nos dois meios (em vermelho).
Esquema da distribuição de carga na interface metal/dielétrico. Também é esquematizado (em vermelho) o decaimento exponencial dos campos elétricos nos meios metálico e dielétrico, perpendicularmente à superfície da interface. As linhas curvas com setas indicam as linhas de força do campo elétrico.
4.3. Relação de dispersão na forma
Uma relação de dispersão é representada, em geral, como uma relação entre a energia, expressa em termos da frequência angular , e o momento linear, expresso em termos do vetor de onda . Queremos obter, para o SPP, uma equação mais explícita do tipo . Para isso, vamos utilizar o modelo de Drude, o qual, apesar de sua simplicidade, fornece uma boa descrição das propriedades óticas e de transporte de metais e semicondutores dopados. Neste modelo a interação entre os elétrons de condução do material é negligenciada e os mesmos são espalhados de maneira aleatória pelos íons da rede com uma frequência característica , onde é um tempo de relaxação médio [19[19] A.V. Andrade-Neto, Revista Brasileira de Ensino de Física 39, e2304 (2017).]. O modelo de Drude fornece a seguinte expressão para a função dielétrica de um gás de elétrons [19[19] A.V. Andrade-Neto, Revista Brasileira de Ensino de Física 39, e2304 (2017).]
onde é a constante dielétrica a altas frequência do material, é a frequência de plasma do gás de elétrons dada pela Eq.(1) e , como definido acima, representa uma constante de amortecimento. Separando a equação acima nas suas partes real e imaginária obtemos
Para metais típicos temos que . Para a prata, por exemplo, temos que e [18[18] William L. Barnes, J. Opt. A: Pure Appl. Opt. 8, S87 (2006).]. Isso significa que podemos considerar apenas a parte real da função dielétrica para a obtenção da relação de dispersão. Usando a equação (52) na equação (38) temos que
Podemos reescrever a equação (54) na seguinte forma
Se o meio dielétrico for o vácuo () e fazendo e a equação (55) torna-se
que representa a relação de dispersão para o SPP. De forma semelhante, podemos obter a relação de dispersão para o plásmon-polaríton de volume. Da condição de transversalidade da onda eletromagnética, obtemos das equações de Maxwell a seguinte expressão, conforme equação (A.6) do Apêndice A Apêndice A Como vimos na Introdução, as variadas propriedades físicas de um sólido podem ser descritas em termos do conceito fundamental de quasipartículas ou excitações elementares. A interação entre diferentes tipos de excitações elementares dá origem a novas quasipartículas que resultam do acoplamento entre estas (os modos híbridos ou modos acoplados). Como exemplo típico podemos citar os diferentes tipos de poláritons. Contudo, é importante enfatizar que só pode haver acoplamento entre modos de vibração paralelos, i.e., para haver acoplamento ambos os modos devem ser longitudinais ou ambos devem ser transversais. Assim, o fônon-poláriton resulta do acoplamento entre o fóton e o fônon transversal ótico. Não é possível, por exemplo, haver acoplamento entre o fóton e o fônon longitudinal ótico. As excitações elementares podem, em grande parte, ser entendidas a partir da análise das equações de Maxwell em um meio material, escrita em termos de componentes de Fourier (ver Apêndice B). A lei de Gauss, equação (B.16) pode ser escrita como (A.1) ϵ ( k → , ω ) k → ⋅ E → ( k → , ω ) = 0 . Da equação (A.1) temos duas soluções possíveis: (A.2) ϵ ( k → , ω ) = 0 e (A.3) k → ⋅ E → ( k → , ω ) = 0 . A primeira solução, equação (A.2), está associada as oscilações longitudinais coletivas do plasma que, na linguagem de segunda quantização, constituem os plásmons que são determinados pelas raízes de ϵ(k→,ω). Usando o modelo de Drude e desprezando amortecimento (Γ=0), vemos da equação (51) que a condição ϵ(k→,ω)=0 ocorre para ω=ωpl/ϵ. Como os plásmons de volume correspondem a modos longitudinais, não pode haver acoplamento entre estes e os fótons. A segunda solução, equação (A.3), define a condição de transversalidade. Das equações de Maxwell, na ausência de fontes, obtemos a seguinte equação (A.4) ∇ × ( ∇ × E → ) = - 1 c 2 ∂ 2 ∂ t 2 ( ϵ E → ) , a qual, em termos de componentes de Fourier, pode ser escrita como (A.5) ( k → ⋅ E → ( k → , ω ) ) k → - k 2 E → ( k → , ω ) = - ω 2 c 2 ϵ ( k → , ω ) E → ( k → , ω ) . Utilizando (A.3) em (A.5) obtemos (A.6) k 2 = ω 2 c 2 ϵ ( k → , ω ) , que corresponde à relação de dispersão para o fóton se propagando em um meio material. Vemos que essa propagação é mediada pela interação entre o fóton e os elétrons de condução do material. O fóton revestido por essa interação dá origem a uma nova quasipartícula denominada plásmon-poláriton de volume. ,
Utilizando (52) em (57), e considerando , obtemos a relação de dispersão para o plásmon-poláriton de volume.
Na Figura 4 são mostradas as relações de dispersão para o plásmon-poláriton de volume (BPP) e para o plásmon-poláriton de superfície (SPP) em uma interface metal dielétrico. O eixo da frequência está normalizado pela frequência de plasma, enquanto o eixo horizontal está em unidades de . A relação de dispersão para o BPP está representada pela curva em vermelho (obtida da equação (58)). Por sua vez, a relação de dispersão do SPP, obtida da equação (56), é a curva em preto. A relação de dispersão para a luz no vácuo (linha verde) é dada pela relação linear, e é denominada linha da luz. Finalmente, a linha azul representa a relação de dispersão para o plasmon de superfície () e é denominada linha do plasmon de superfície.
Relações de dispersão do plásmon-poláriton de volume (curva vermelha), para a luz no vácuo (linha verde), o plásmon de superfície (linha azul) e para o plásmon polaríton de superfície (curva em preto).
Muitas das propriedades físicas dos SPPs podem ser entendidas a partir da análise da Figura 4, conforme faremos a seguir.
Fisicamente, o SPP pode ser intepretado como sendo o resultado do acoplamento entre as oscilações coletivas dos elétrons na superfície (cujos quanta são os plásmons de superfície) e uma onda eletromagnética superficial (cujos quanta são os fótons). Na interface metal/dielétrico ocorre uma forte interação entre os modos fotônico e plasmônico representados, respectivamente, pela linha de luz e linha de plasmon de superfície na Figura 4. Essa interação elimina o cruzamento entre esses modos dando origem a um novo modo que possui uma natureza híbrida fóton-plásmon o qual, por essa razão, é denominado Plásmon-Polarítons de Superfície.
Podemos destacar duas regiões na Figura 4. Para valores pequenos de vetor de onda, , a curva de dispersão do SPP está muito proxima da linha da luz. Isto significa que nesta região os SPPs possuem um caráter fotônico. A medida que aumenta, a curva de dispersão se afasta cada vez mais da linha de luz e se aproxima assintoticamente da frequência quando .
Na região oposta, , como já dito, a curva de dispersão tende assintoticamente para o valor que define a linha do plásmon de superfície. Assim, na região o SPP possui um caráter plasmônico. É interessante notar que podemos obter o mesmo resultado a partir da equação (54). Dessa equação vemos que, para , quando a frequência assume o valor
a qual, para , assume o valor .
5. Conclusões
Neste trabalho realizamos uma introdutória apresentação pedagógica da emergente disciplina da plasmônica, a qual estuda a interação entre uma onda eletromagnética e os elétrons de condução de uma superfície metálica na inteface entre um metal e um dielétrico, incluindo nanoestruturas metálicas. Sob certas condições, essa interação produz o confinamento da luz em dimensões da ordem ou menores que seu comprimento de onda. Aqui, partindo das equações de Maxwell, deduzimos em detalhes a relação de dispersão para o plásmon-poláriton de superfície na interface plana entre um metal e um dielétrico. Conseguimos obter relações gerais, as quais recuperam expressões conhecidas da literatura. Também obtivemos uma expressão explícita para a relação de dispersão, versus , utilizando a expressão de Drude para a função dielétrica de um gás clássico de elétrons. Adicionalmente, procuramos fazer uma discussão detalhada sobre o significado físico dos resultados obtidos.
Apêndice A
Como vimos na Introdução, as variadas propriedades físicas de um sólido podem ser descritas em termos do conceito fundamental de quasipartículas ou excitações elementares. A interação entre diferentes tipos de excitações elementares dá origem a novas quasipartículas que resultam do acoplamento entre estas (os modos híbridos ou modos acoplados). Como exemplo típico podemos citar os diferentes tipos de poláritons. Contudo, é importante enfatizar que só pode haver acoplamento entre modos de vibração paralelos, i.e., para haver acoplamento ambos os modos devem ser longitudinais ou ambos devem ser transversais. Assim, o fônon-poláriton resulta do acoplamento entre o fóton e o fônon transversal ótico. Não é possível, por exemplo, haver acoplamento entre o fóton e o fônon longitudinal ótico.
As excitações elementares podem, em grande parte, ser entendidas a partir da análise das equações de Maxwell em um meio material, escrita em termos de componentes de Fourier (ver Apêndice B Apêndice B As equações de Maxwell em um meio material, utilizando o sistema cgs, são escritas como (B.1) ∇ × H → ( r → , t ) = 1 c ∂ D → ( r → , t ) ∂ t + 4 π c J → c o n ( r → , t ) (B.2) ∇ × E → ( r → , t ) = - 1 c ∂ B → ( r → , t ) ∂ t (B.3) ∇ ⋅ D → ( r → , t ) = 4 π ρ e x t ( r → , t ) (B.4) ∇ ⋅ B → ( r → , t ) = 0 , onde D→(r→,t) é o vetor deslocamento elétrico, J→con(r→,t) é a densidade de corrente de condução e ρext(r→,t) é a densidade volumétrica de carga. Muitas vezes é conveniente escrever as equações de Maxwell em termos de componentes de Fourier. Assumindo um meio sem dispersão espacial podemos escrever o campo E→(r→,t) na seguinte forma (B.5) E → ( r → , t ) = 1 2 π ∫ - ∞ ∞ E → ( r → , ω ) e - ı ω t d ω , onde E→(r→,ω) é a transformada de Fourier de E→(r→,t) dada por (B.6) E → ( r → , ω ) = ∫ - ∞ ∞ E → ( r → , ω ) e ı ω t d t , com expresões semelhantes para os outros campos. Diferenciando o campo elétrico temos que (B.7) ∇ × E → ( r → , t ) = 1 2 π ∫ - ∞ ∞ ∇ × E → ( r → , ω ) e - ı ω t d ω (B.8) ∇ ⋅ E → ( r → , t ) = 1 2 π ∫ - ∞ ∞ ∇ ⋅ E → ( r → , ω ) e - ı ω t d ω (B.9) ∂ E → ( r → , t ) ∂ t = 1 2 π ∫ - ∞ ∞ ( - ı ω ) E → ( r → , ω ) e - ı ω t d ω Usando as equações acima nas equações de Maxwell, estas ficam (B.10) ∇ × H → ( r → , ω ) = - ı ω c D → ( r → , ω ) + 4 π c J → c o n ( r → , ω ) (B.11) ∇ × E → ( r → , ω ) = ı ω c B → ( r → , ω ) (B.12) ∇ ⋅ D → ( r → , ω ) = 4 π ρ e x t ( r → , ω ) (B.13) ∇ ⋅ B → ( r → , ω ) = 0 Para um meio não-magnético (B→=H→) e na ausência de fontes temos que (B.14) ∇ × H → ( r → , ω ) = - ı ω c ϵ ( ω ) E → ( r → , ω ) (B.15) ∇ × E → ( r → , ω ) = ı ω c H → ( r → , ω ) (B.16) ∇ ⋅ ϵ ( ω ) E → ( r → , ω ) = 0 (B.17) ∇ ⋅ H → ( r → , ω ) = 0 onde usamos que D→(r→,ω)=ϵ(ω)E→(r→,ω). ).
A lei de Gauss, equação (B.16) pode ser escrita como
Da equação (A.1) temos duas soluções possíveis:
e
A primeira solução, equação (A.2), está associada as oscilações longitudinais coletivas do plasma que, na linguagem de segunda quantização, constituem os plásmons que são determinados pelas raízes de . Usando o modelo de Drude e desprezando amortecimento (), vemos da equação (51) que a condição ocorre para . Como os plásmons de volume correspondem a modos longitudinais, não pode haver acoplamento entre estes e os fótons. A segunda solução, equação (A.3), define a condição de transversalidade.
Das equações de Maxwell, na ausência de fontes, obtemos a seguinte equação
a qual, em termos de componentes de Fourier, pode ser escrita como
Utilizando (A.3) em (A.5) obtemos
que corresponde à relação de dispersão para o fóton se propagando em um meio material. Vemos que essa propagação é mediada pela interação entre o fóton e os elétrons de condução do material. O fóton revestido por essa interação dá origem a uma nova quasipartícula denominada plásmon-poláriton de volume.
Apêndice B
As equações de Maxwell em um meio material, utilizando o sistema cgs, são escritas como
onde é o vetor deslocamento elétrico, é a densidade de corrente de condução e é a densidade volumétrica de carga.
Muitas vezes é conveniente escrever as equações de Maxwell em termos de componentes de Fourier.
Assumindo um meio sem dispersão espacial podemos escrever o campo na seguinte forma
onde é a transformada de Fourier de dada por
com expresões semelhantes para os outros campos.
Diferenciando o campo elétrico temos que
Usando as equações acima nas equações de Maxwell, estas ficam
Para um meio não-magnético () e na ausência de fontes temos que
onde usamos que .
Agradecimentos
A.V. A-N agradece à Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) pelo suporte financeiro para o desenvolvimento deste trabalho (proc. 23038.003304/2016-72). A. R-N. agradece ao CNPq pela bolsa de iniciação científica. A. J. agradece financiamento do CNPq (552124/2011-7) e FINEP (01.13.0330.00).
Referências
-
[1]Marvin L. Cohen, Phys. Rev. Lett. 101, 250001 (2008).
-
[2]Ado Jorio, Bioengineering Applications of Carbon Nanostructures (Springer Verlag, Berlin, 2016).
-
[3]M. Dragoman and D. Dragoman, Progress in Quantum Electronic 32, 1 (2008).
-
[4]Ricardo S. Oyarzabal, José D. Szezech Junior, Antonio M. Batista, Iberê L. Caldas, Ricardo L. Viana e Kelly C. Iarosz, Revista Brasileira de Ensino de Física 37, 2308 (2015).
-
[5]P.M. Platzman and P.A. Wolff, in: Solid State Physics, Suppl. 13, edited by F. Seitz, D. Turnbull and H. Ehrenrich (Academic Press, New York, 1973).
-
[6]Lukas Novotny and Bert Hecht, Principles of Nano-Optics (Cambridge University Press, Cambridge, 1999), 2ª ed, p. 5, p. 371.
-
[7]Rodrigo Gonçalves Pereira e Eduardo Miranda, Revista Brasileira de Ensino de Física 24, 237 (2002).
-
[8]Indianara Silva, Revista Brasileira de Ensino de Física 37, 4204 (2015).
-
[9]Richard D. Mattuck, A Guide to Feynman Diagrams in the Many-body Problem (McGraw-Hill, New York, 1967).
-
[10]S. Nakajima, Y. Toyozawa and R. Abe, The Physics of Elementary Excitations (Springer-Verlag, Berlin, 1980).
-
[11]Harry A. Atwater, Scientific American 296, 56 (2007).
-
[12]Francisco E.G. Guimarães, Revista Brasileira de Ensino de Física 18, 56 (1996).
-
[13]Daniel M. Reis, Edson M. Santos e A.V. Andrade-Neto, Revista Brasileira de Ensino de Física 37, 2312 (2015).
-
[14]Thiago L. Vasconcelos, Bráulio S. Archanjo, Benjamin Fragneaud, Bruno S. Oliveira, Juha Riikonen, Changfeng Li, Douglas S. Ribeiro, Cassiano Rabelo, Wagner N. Rodrigues, Ado Jorio, et. al., ACSNano 9, 6297 (2015).
-
[15]John R. Reitz, Frederick J. Milford e Robert W. Christy, Fundamentos da Teoria Eletromagnética (Editora Campus, São Paulo, 1982), p. 95, p. 199.
-
[16]Stefan Alexander Maier, Plasmonics: Fundamentals and Applications (Springer, New York, 2007), p. 10-27.
-
[17]H. Raether, Surface Plasmons on Smooth and Rough Surfaces and on Gratings (Springer-Verlag, Berlin, 1998), Appendix I.
-
[18]William L. Barnes, J. Opt. A: Pure Appl. Opt. 8, S87 (2006).
-
[19]A.V. Andrade-Neto, Revista Brasileira de Ensino de Física 39, e2304 (2017).
Datas de Publicação
-
Publicação nesta coleção
2017
Histórico
-
Recebido
14 Nov 2016 -
Revisado
13 Fev 2017 -
Aceito
14 Fev 2017