Ultrafast optically induced tunneling in narrow metallic gaps from the time dependent density functional perspective

1. Problema e Contexto

O tunelamento de elétrons através de barreiras de potencial é um efeito quântico fundamental com aplicações em eletrônica e microscopia de tunelamento (STM). Recentemente, a combinação da resolução espacial atômica do STM com a resolução temporal de pulsos ópticos curtos deu origem à Microscopia de Tunelamento Ultrafast (USTM).

O desafio central abordado neste trabalho é compreender e modelar quantitativamente o transporte de elétrons induzido opticamente em gaps metálicos estreitos (na escala sub-nanométrica a nanométrica). Enquanto teorias perturbativas (como a teoria de Tien-Gordon) funcionam bem para campos ópticos fracos e gaps largos (onde o transporte é dominado por emissão sobre a barreira), elas falham em descrever regimes de campos fortes e pulsos ultracurtos (de um único ciclo ou poucos ciclos). Além disso, muitas abordagens existentes dependem de parâmetros ajustáveis ou simplificam excessivamente a interação elétron-elétron, dificultando uma análise quantitativa robusta sem ajustes experimentais.

O objetivo é identificar os canais de tunelamento assistido por fótons (PAT), entender a transição entre regimes de tunelamento e emissão de campo óptico, e explicar dados experimentais recentes sem utilizar parâmetros de ajuste.

2. Metodologia

Os autores empregam uma abordagem multiescala, combinando métodos de muitos corpos com modelos analíticos e numéricos simplificados:

• Teoria do Funcional da Densidade Dependente do Tempo (TDDFT):

  • É o método principal, utilizado para cálculos ab initio (sem parâmetros de ajuste).
  • O sistema modelado consiste em um dímero de nanofios cilíndricos metálicos paralelos (raio $R=5$ nm) separados por um gap estreito ($d_{gap}$ variando de 0,8 nm a 2 nm).
  • Utiliza-se o modelo de "jellium" (elétrons livres) para descrever a densidade eletrônica, com aproximação de densidade local adiabática (ALDA).
  • O sistema é excitado por pulsos ópticos polarizados (ciclo único e vários ciclos) e, em alguns casos, submetido a um viés de corrente contínua (DC) para quebrar a simetria.
  • A dinâmica é resolvida via equações de Kohn-Sham dependentes do tempo, calculando a densidade de corrente eletrônica e a transferência líquida de elétrons.

• Abordagem de Um Elétron Ativo (WPP - Wave-Packet Propagation):

  • Uma solução numérica da equação de Schrödinger dependente do tempo (TDSE) em 1D para um único elétron ativo.
  • Utiliza potenciais modelo que incluem interações de imagem clássicas (que a TDDFT padrão com ALDA não captura corretamente em longas distâncias).
  • Serve para fornecer insights intuitivos sobre a dinâmica de partículas individuais e espectroscopia de estados finais.

• Teoria de Campo Forte Semiclássica (SFT):

  • Um modelo analítico baseado na aproximação de campo forte (SFA) e no método do ponto de sela.
  • Fornece uma descrição semi-clássica dos processos de emissão e tunelamento, permitindo a definição de parâmetros como o parâmetro de Keldysh ($\gamma$) e um novo parâmetro $\zeta$ para distinguir regimes de transporte.

• Consistência entre Métodos:

  • Os autores ajustam o trabalho de saída ($\Phi$) nos cálculos TDDFT para compensar a falta do potencial de imagem de longo alcance na ALDA, garantindo que a altura da barreira de tunelamento corresponda à realidade física e seja consistente com os modelos WPP e SFT.

3. Contribuições Principais

• Identificação de Canais de Tunelamento: Mapeamento claro dos canais de tunelamento assistido por fótons (PAT) de 1, 2 e múltiplos fótons em função da intensidade do campo óptico e da largura do gap.
• Transição de Regimes: Descrição detalhada da transição do regime perturbativo (tunelamento assistido por fótons) para o regime de emissão de campo óptico (tunelamento dominado pelo campo) e a transição para processos "over-the-barrier" (acima da barreira).
• Modelo Livre de Parâmetros: Demonstração de que a TDDFT, quando configurada corretamente, pode reproduzir dados experimentais quantitativos sem necessidade de parâmetros ajustáveis, validando o mecanismo físico subjacente.
• Papel do Viés DC: Análise de como um viés DC quebra a simetria do sistema, permitindo a observação do canal de tunelamento de um único fóton ($\ell=1$) em gaps simétricos, o que é crucial para a comparação com experimentos de STM.

4. Resultados Chave

Dinâmica do Transporte e Regimes

• Regime Perturbativo (Campos Fracos): Para gaps simétricos sem viés DC, a transferência líquida de elétrons é zero para absorção de um único fóton devido à simetria. O canal dominante é a absorção de dois fótons ($\ell=2$). A dependência da corrente com o campo segue $N \propto E_g^{2\ell}$.
• Efeito do Viés DC: A aplicação de um viés DC quebra a simetria, permitindo que o canal de um fóton ($\ell=1$) domine em campos baixos. Isso explica resultados experimentais onde a não-linearidade da corrente muda de $n \approx 1$ para $n \approx 2$ ou superior à medida que o campo aumenta.
• Parâmetro $\zeta$: Os autores utilizam o parâmetro $\zeta$ (relacionado ao tempo imaginário de emissão) para definir a fronteira entre o regime de tunelamento perturbativo ($\zeta > 1$) e o regime onde excitações sobre a barreira se tornam possíveis ($\zeta < 1$).
• Ordem Eficiente de Fótons ($n$): Em regimes de múltiplos fótons, a ordem efetiva de fótons $n$ extraída da dependência de potência da corrente converge rapidamente para o limiar de fotoemissão sobre a barreira ($n \approx W_{tb}/\hbar\omega$), onde $W_{tb}$ é a altura da barreira de tunelamento. Isso significa que, em campos fortes, a corrente não reflete mais um único canal de absorção, mas uma média sobre vários canais.

Dependência da Geometria e Campo

• Largura do Gap: Aumentar a largura do gap suprime o tunelamento de baixa energia (canais de 1 e 2 fótons), tornando o transporte dominado por processos de múltiplos fótons com energias próximas ao topo da barreira.
• Controle de CEP (Carrier Envelope Phase): Em pulsos de ciclo único, a transferência líquida de elétrons depende do CEP, mas essa dependência desaparece no regime de um único fóton em sistemas simétricos e só é significativa no regime de múltiplos fótons ($\ell \ge 2$).

Comparação com Experimentos

• Antenas Bowtie de Ouro: Os resultados TDDFT reproduzem qualitativamente e semi-quantitativamente os dados experimentais de Luo et al. (Nature Communications, 2024) sobre antenas bowtie. O modelo prevê corretamente a transição de $n=1$ para $n=2$ e a dependência com o viés DC.
• Junções STM de Prata: Os resultados também concordam com experimentos de Lin et al. (ACS Photonics, 2023) em junções STM de prata, reproduzindo a variação da corrente com o viés e o tamanho do gap.
• Discrepâncias: Pequenas discrepâncias quantitativas (ordem de magnitude) são atribuídas a efeitos não capturados pelo modelo, como o "ringing" de plásmons (que aumenta a duração efetiva do pulso no gap) e relaxação de elétrons quentes (termo "lento" de emissão), que estão além do escopo da ALDA-TDDFT.

5. Significado e Conclusão

Este trabalho estabelece uma base teórica robusta e livre de parâmetros para entender a eletrônica de petahertz e a microscopia de tunelamento ultrafast.

• Validação de Mecanismos: Confirma que o transporte de elétrons em gaps metálicos estreitos é um processo coerente e rápido, onde a excitação e o tunelamento não podem ser separados temporalmente.
• Guia para Experimentos: Fornece diretrizes claras sobre como controlar o transporte de elétrons variando o tamanho do gap, o viés DC e a forma de onda do pulso óptico.
• Limitações e Futuro: O estudo destaca que, embora a TDDFT capture a componente "rápida" e coerente do transporte, a descrição completa de processos experimentais exige futuramente a incorporação de efeitos de relaxação (espalhamento elétron-elétron e elétron-fônon) e estrutura de bandas, que foram negligenciados no modelo de elétrons livres.

Em suma, o artigo demonstra que a TDDFT é uma ferramenta poderosa para prever e explicar fenômenos complexos de tunelamento óptico em nanoestruturas, fechando a lacuna entre teorias analíticas simplificadas e dados experimentais complexos.