First-principles electron-phonon scattering in real-time TDDFT

Resumo Técnico: Espalhamento Elétron-Fônon de Primeiros Princípios em TDDFT em Tempo Real

Enunciado do Problema
A teoria do funcional da densidade dependente do tempo em tempo real (rt-TDDFT) é uma ferramenta estabelecida de primeiros princípios para simular a dinâmica coerente de elétrons em sólidos excitados por laser. No entanto, sua formulação padrão baseia-se em uma evolução unitária dos orbitais de Kohn-Sham, que descreve um sistema quântico fechado. Consequentemente, a rt-TDDFT convencional não consegue capturar processos irreversíveis, como o relaxamento de portadores, a termalização e a decoerência decorrentes do acoplamento de elétrons excitados ao seu ambiente. Em materiais reais, as interações elétron-fônon (e-ph) são o mecanismo dominante para esses processos dissipativos, contudo, integrar essas interações em um arcabouço de primeiros princípios em tempo real, mantendo o acesso a observáveis espectroscópicos, tem sido um desafio central. Abordagens existentes frequentemente dependem de modelos fenomenológicos ou equações de Boltzmann semiclássicas que rastreiam apenas populações de portadores, falhando em descrever o decaimento das coerências quânticas essenciais para a dinâmica ultrarrápida.

Metodologia
Os autores desenvolvem um arcabouço de rt-TDDFT dissipativo que trata o subsistema eletrônico como um sistema quântico aberto fracamente acoplado a um banho de fônons térmico. A metodologia integra três componentes principais:

1. Formulação da Matriz de Densidade de Sistema Aberto: Partindo de elementos de matriz e-ph de primeiros princípios, os autores derivam um integral de colisão para a matriz de densidade de um corpo reduzida dentro da aproximação de Born–Markov. A equação de movimento para a matriz de densidade reduzida $\hat{\rho}^{(1)}$ combina a evolução de Liouville–von Neumann coerente, impulsionada pelo Hamiltoniano de Kohn-Sham dependente do tempo ($\hat{H}_{KS}$), com um termo dissipativo $I_{e-ph}$:
   $$\frac{\partial \hat{\rho}^{(1)}(t)}{\partial t} = -i [\hat{H}_{KS}(t), \hat{\rho}^{(1)}(t)] + I_{e-ph}[\hat{\rho}^{(1)}(t)]$$
   O integral de colisão $I_{e-ph}$ é construído para descrever simultaneamente a redistribuição de população assistida por fônons (termos de ganho-perda) e o decaimento das coerências interbandas (amortecimento off-diagonal), garantindo que a evolução permaneça completamente positiva e preservadora de traço (CPTP).

2. Propagação Híbrida Espaço Real/Espaço Ativo: Para manter a viabilidade computacional enquanto preserva todas as capacidades da rt-TDDFT em espaço real, os autores empregam uma representação híbrida.

   • Dinâmica Coerente: A propagação unitária dos orbitais de Kohn-Sham é realizada em espaço real usando o código Octopus, capturando interações luz-matéria não lineares e efeitos de continuum.
   • Atualizações Dissipativas: O integral de colisão é aplicado em um "espaço ativo" selecionado definido por uma base adiabática dependente do tempo (autovetores instantâneos de $\hat{H}_{KS}$). A matriz de densidade é projetada nesta base, atualizada via integral de colisão e, em seguida, mapeada de volta para a representação de orbitais em espaço real.
   • Construção da Base: Para evitar perda de norma e ambiguidades de calibre (gauge), o espaço de Hilbert é particionado em um espaço ativo ($P$) e um espaço residual ($Q$). Um procedimento de fixação de calibre baseado no algoritmo de Procrustes ortogonal garante a continuidade dos orbitais entre os passos de tempo.

3. Taxas de Transição de Granularidade Grossa: Como o espalhamento e-ph de primeiros princípios requer grades de momento ultra-finas para resolver energias de fônons na escala de meV (proibitivo para a propagação em tempo real), os autores implementam um esquema de granulação grossa (coarse-graining). Taxas de transição microscópicas calculadas em uma grade densa (usando o código EPW) são médias sobre patches de momento para definir taxas efetivas para a grade mais grosseira usada na simulação de rt-TDDFT.

Resultados Principais
O arcabouço foi validado através de duas aplicações de referência:

• Relaxação de Portadores Ultrarrápida em Silício Bulk:

  • Ao simular um pulso de laser de 2,5 eV, o método capturou a redistribuição dos portadores fotoexcitados.
  • Relaxação de Energia vs. Momento: Os resultados revelaram uma separação de escalas de tempo: a redistribuição de momento (espalhamento para fora dos vales) ocorreu em uma escala de tempo de 8–14 fs, enquanto o relaxamento de energia (resfriamento) foi mais lento (40–106 fs). Isso apoia a formação de um ensemble de elétrons quentes de quase-equilíbrio antes do resfriamento macroscópico.
  • Concordância Quantitativa: Os tempos de vida de energia extraídos (40 fs a 2,6 eV; 106 fs a 1,6 eV) mostraram concordância quantitativa com medições experimentais, validando a abordagem dissipativa.

• Dinâmica de Vales em WS$_2$ Monocamada:

  • O estudo modelou o espalhamento entre vales após excitação por luz circularmente polarizada (CP) e linearmente polarizada (LP).
  • Despolarização de Vale: Sob luz CP, a polarização de vale inicial (vale K) decaiu conforme os portadores se espalhavam para os vales K', Q e Q'. O tempo de espalhamento entre vales extraído foi de 20 fs (17 fs para LP), consistente com valores experimentais recentes (16 fs).
  • Conexão com tr-ARPES: Ao acoplar a dinâmica dissipativa com o framework de fotoemissão tSURFF, os autores geraram espectros de tr-ARPES totalmente ab initio. Essas simulações reproduziram com sucesso a excitação seletiva de vale e a subsequente redistribuição, ligando taxas de espalhamento microscópicas a sinais de fotoemissão mensuráveis.

Significância e Alegações
O artigo afirma fornecer uma rota prática de primeiros princípios para simular dinâmica quântica dissipativa em materiais cristalinos realistas sem recorrer a parâmetros fenomenológicos. Suas principais contribuições são:

1. Descrição Unificada de Relaxação e Decoerência: Diferente das abordagens de Boltzmann semiclássicas que rastreiam apenas populações, este arcabouço retém a matriz de densidade de partícula única completa, permitindo a descrição simultânea do relaxamento de população e do decaimento de coerência impulsionados pelos mesmos canais de espalhamento microscópicos.
2. Integração com Espectroscopia: Ao incorporar a atualização dissipativa dentro do arcabouço de rt-TDDFT em espaço real do código Octopus, o método preserva o acesso a todo o conjunto de observáveis, incluindo respostas ópticas não lineares e assinaturas espectroscópicas resolvidas no tempo (ex: tr-ARPES).
3. Escalabilidade: A formulação de espaço ativo limita o custo computacional do integral de colisão ao subconjunto relevante de estados, tornando a abordagem viável para materiais realistas enquanto mantém a conexão rigorosa com a teoria e-ph de primeiros princípios.

Os autores posicionam este trabalho como uma realização controlada de um arcabouço geral de rt-TDDFT de sistema aberto, capaz de abordar fenômenos ultrarrápidos onde o acionamento coerente e o espalhamento dissipativo competem, como em geração de harmônicos de alta ordem e espectroscopia de bomba-sonda. Eles observam que, embora a implementação atual assuma um banho de fônons térmico e aproximações de Born–Markov, a estratégia é extensível a outros canais ambientais (ex: fótons, magnons) e regimes envolvendo efeitos não-markovianos fortes.