Analytic first order nonadiabatic coupling matrix elements of spin-adapted open-shell time-dependent density functional theory

Este artigo apresenta a derivação, implementação e avaliação de desempenho dos elementos de matriz de acoplamento não adiabático de primeira ordem analíticos para o método X-TDDFT adaptado ao spin, demonstrando que ele reduz significativamente os erros em comparação com o U-TDDFT padrão e fornece insights qualitativamente corretos sobre a fotofísica de sistemas de camada aberta, como a porfirina de cobre(II).

O essencial

Imagine que você está tentando prever como uma máquina complexa, como um pião giratório feito de ímãs, reagirá quando você o agitar. No mundo da química, essa "máquina" é uma molécula com elétrons desemparelhados (como radicais orgânicos ou complexos de metais de transição), e o "agitar" é a luz incidindo sobre ela, fazendo-a saltar para um estado excitado.

Os cientistas utilizam uma ferramenta chamada Teoria do Funcional da Densidade Dependente do Tempo (TDDFT) para simular essas reações. Pense na TDDFT como uma previsão do tempo sofisticada para moléculas. Ela prevê como a molécula se move e como sua energia muda.

No entanto, há um problema. A TDDFT padrão (vamos chamá-la de U-TDDFT) é como uma previsão do tempo que assume que o vento sempre sopra em linha reta. Funciona razoavelmente bem para moléculas simples, mas para moléculas complexas com elétrons "desemparelhados" (como nosso pião de ímã giratório), ela fica confusa. Ela trata os dois "spins" dos elétrons (vamos chamá-los de Spin A e Spin B) como se fossem independentes, o que leva a erros. É como tentar descrever uma dança onde dois parceiros estão de mãos dadas, mas a previsão assume que eles estão dançando sozinhos.

A Nova Solução: X-TDDFT

Os autores deste artigo desenvolveram uma atualização chamada X-TDDFT. Isso é como um novo modelo meteorológico que entende que os parceiros estão de mãos dadas. Ele força a matemática a respeitar as regras de "spin" da mecânica quântica. Eles já haviam usado isso para prever melhor a energia e a forma dessas moléculas, mas faltava uma peça crucial: os Elementos de Matriz de Acoplamento Não Adiabático (NACMEs).

O que é um NACME?
Imagine que a molécula é um carro dirigindo em uma estrada acidentada.

• A Energia diz o quão rápido o carro está indo.
• Os Gradientes dizem para onde a estrada está inclinada.
• Os NACMEs dizem quão provável é o carro trocar de faixa ou bater em um estado diferente.

Na química, essa "troca de faixa" é chamada de Conversão Interna (CI). É o processo em que uma molécula absorve energia, fica excitada e depois rapidamente libera essa energia de volta para o estado fundamental, frequentemente liberando calor em vez de luz. Se o seu cálculo de NACME estiver errado, você pode achar que o carro permanecerá em sua faixa, quando, na realidade, ele está prestes a desviar violentamente para uma vala.

O que eles fizeram?

A equipe derivou as fórmulas matemáticas para calcular essas probabilidades de "troca de faixa" (NACMEs) usando seu novo método X-TDDFT, consciente do spin. Eles então o testaram de duas maneiras:

1. O Teste Pequeno (Radical de Formaldeído): Eles compararam seu novo método contra um cálculo superpreciso "padrão-ouro" (como verificar um novo GPS contra um mapa de satélite). Eles descobriram que o método antigo (U-TDDFT) frequentemente errava por uma margem enorme — às vezes fora em um terço a dois terços. O novo método (X-TDDFT) corrigiu a maioria desses erros, tornando a previsão de quão rápido a molécula "resfria" (taxa de Conversão Interna) muito mais precisa. Em alguns casos, o novo método previu que a velocidade de resfriamento era 100 vezes mais lenta do que o método antigo previu.

2. O Teste Grande (Porfirinas de Cobre): Eles analisaram moléculas complexas à base de cobre (semelhantes à heme no sangue, mas com cobre).

   • A Visão Antiga (U-TDDFT): Previa que, quando a molécula fica excitada, ela tem chances iguais de resfriar diretamente ou fazer um desvio através de estados intermediários.
   • A Nova Visão (X-TDDFT): Previa que a molécula quase nunca resfria diretamente. Ela quase sempre faz o desvio.
   • O Resultado: Isso mudou completamente a história de como essas moléculas se comportam. O método antigo não apenas errou ligeiramente os números; ele errou a história. Também atrapalhou a comparação entre diferentes versões da molécula (com diferentes decorações químicas), fazendo parecer que uma versão era mais rápida que a outra quando o oposto era verdadeiro.

A Conclusão

O artigo conclui que, para moléculas com elétrons desemparelhados (como radicais ou metais de transição), você não pode confiar na antiga matemática de "linha reta" (U-TDDFT) para prever como elas mudam de estados de energia.

Assim como você não usaria um mapa plano para navegar em uma cadeia de montanhas, você não deve usar o antigo método TDDFT para essas moléculas complexas. O novo método X-TDDFT age como um mapa topográfico 3D, revelando que as "estradas" (caminhos de energia) são muito diferentes do que se pensava anteriormente. Isso é crucial para cientistas tentando projetar melhores células solares, LEDs ou catalisadores, porque, se você não sabe em qual "faixa" a molécula vai saltar, você não pode controlar seu comportamento.

Em resumo: Os autores construíram uma régua melhor para medir como as moléculas "saltam" entre estados de energia. Eles provaram que a régua antiga era tão imprecisa que estava contando histórias completamente diferentes sobre como essas moléculas funcionam, especialmente para aquelas envolvendo cobre e outros metais de transição.

Resumo técnico

Resumo Técnico: Elementos de Matriz de Acoplamento Não Adiabático de Primeira Ordem Analíticos de Teoria do Funcional da Densidade Dependente do Tempo Adaptada ao Spin para Sistemas de Camada Aberta

Enunciado do Problema
A teoria do funcional da densidade dependente do tempo adaptada ao spin (X-TDDFT) foi recentemente estabelecida como um método robusto para o cálculo de energias de excitação e gradientes analíticos de moléculas de camada aberta (como radicais e complexos de metais de transição), oferecendo melhorias significativas em relação à TDDFT não restrita (U-TDDFT) ao lidar adequadamente com a contaminação de spin e incluir implicitamente excitações duplas específicas. No entanto, o efeito da adaptação ao spin sobre os elementos de matriz de acoplamento não adiabático (NACMEs) — críticos para a simulação de taxas de conversão interna (CI) e dinâmicas não adiabáticas — permanecia desconhecido. Embora existam teorias rigorosas de NACMEs em U-TDDFT, elas sofrem de grandes erros devido à contaminação de spin, particularmente para estados dominados por excitações de camada fechada para camada vazia (CV). Teorias existentes de NACMEs adaptadas ao spin são limitadas a Hamiltonianos semiempíricos, excitações específicas de inversão de spin, ou carecem do formalismo necessário para sistemas gerais de camada aberta. Consequentemente, a simulação precisa da fotofísica e fotoquímica em sistemas de camada aberta usando frameworks conservadores de spin foi impedida pela falta de uma teoria consistente de NACMEs em X-TDDFT.

Metodologia
Os autores derivam e implementam NACMEs analíticos de primeira ordem para o método X-TDDFT, especificamente para transições do estado fundamental para o estado excitado (ge) e de estado excitado para estado excitado (ee). A derivação segue uma abordagem de equação de movimento (EOM), substituindo os operadores padrão de excitação e desexcitação da aproximação de fase aleatória (RPA) por seus equivalentes adaptados ao spin.

Componentes teóricos-chave incluem:

1. Operadores Adaptados ao Spin: O formalismo utiliza um estado de referência de Kohn-Sham de camada aberta restrito (ROKS). O operador de excitação $O^\dagger_I$ é construído usando excitações CV(0) e CV(1) adaptadas ao spin, onde a última é misturada com excitações duplas específicas ($\Psi^{\bar{t}a}_{\bar{i}t}$) para garantir a simetria de spin adequada.
2. Formulação Lagrangiana: Para lidar com a dependência implícita dos vetores de excitação e dos coeficientes de orbitais moleculares (OM) nas coordenadas nucleares, emprega-se uma abordagem lagrangiana. Isso permite o cálculo de derivadas analíticas sem diferenciar explicitamente os coeficientes dos OM, incorporando termos de Pulay por meio de multiplicadores de Lagrange (vetores $Z$ e $W$).
3. Termos de Correção: Os NACMEs em X-TDDFT são expressos como os NACMEs padrão em U-TDDFT mais termos de correção específicos decorrentes de:
   • O uso de uma referência adaptada ao spin (ROKS) em vez de uma não restrita (UKS).
   • A modificação do operador de excitação CV(1) para incluir excitações duplas.
   • As mudanças resultantes nas matrizes de densidade de transição (TDMs) e na matriz de Casida (blocos $A$ e $B$).
4. Estratégia de Implementação: A implementação foi projetada para ser construída sobre códigos existentes de NACMEs em U-TDDFT. Os autores demonstram que o custo computacional adicional é aceitável, comparável ao de cálculos de gradientes analíticos, pois envolve principalmente a resolução de equações de vetores Z modificadas e o cálculo de matrizes de densidade intermediárias.

Contribuições Principais

• Derivação Teórica: O artigo fornece a primeira derivação rigorosa de NACMEs analíticos para X-TDDFT conservador de spin, cobrindo tanto transições ge quanto ee.
• Implementação Algorítmica: O método é implementado no pacote de software BDF. Ele aproveita a infraestrutura existente de U-TDDFT, exigindo apenas a adição de termos de correção específicos às matrizes de densidade de transição e aos multiplicadores de Lagrange.
• Validação: Os autores validam o método contra cálculos multirreferenciais de alto nível (XMS-CASPT2) para o cátion radical do formaldeído ($CH_2O^+$).
• Aplicação a Complexos de Metais de Transição: O método é aplicado a porfirinas de cobre(II) (CuP, CuOEP, CuTPP) para simular caminhos de relaxação do estado excitado e taxas de conversão interna.

Resultados

• Precisão da Validação ($CH_2O^+$): Comparados aos dados de referência XMS-CASPT2, os NACMEs em X-TDDFT reduzem o desvio absoluto médio (MAD) das normas dos NACMEs em U-TDDFT em aproximadamente 32% a 65% (dependendo do funcional). As melhorias mais significativas são observadas para estados com caráter CV(1) significativo, onde a U-TDDFT frequentemente falha devido à contaminação de spin. No entanto, os autores observam que, para certos estados mistos (ex.: $1^2A_2$), a X-TDDFT pode ocasionalmente produzir erros maiores do que a U-TDDFT, embora isso seja atribuído a uma cancelamento fortuito de erros nos resultados da U-TDDFT e não a uma falha da teoria adaptada ao spin.
• Porfirinas de Cobre(II):
  • Magnitudes dos NACMEs: A X-TDDFT prevê magnitudes de NACMEs significativamente menores (até uma ordem de grandeza) para a transição $2T_1 \to 2S_0$ em comparação com a U-TDDFT. Essa redução é atribuída ao cancelamento quase perfeito das contribuições $\alpha$ e $\beta$ em estados CV(1) puros dentro do framework adaptado ao spin, um cancelamento que é incompleto na U-TDDFT devido a diferenças na forma dos orbitais.
  • Taxas de Conversão Interna: A redução nas normas dos NACMEs leva a reduções drásticas nas taxas de CI calculadas (até duas ordens de grandeza).
  • Caminhos de Relaxação: Para CuP não substituído, a U-TDDFT prevê que o caminho de relaxação dominante envolve CI direta aproximadamente igual ao estado fundamental e CI indireta via estados de campo de ligante. Em contraste, a X-TDDFT prevê que o caminho indireto via estados de campo de ligante é dominante.
  • Efeitos de Substituintes: As tendências relativas das taxas de CI para porfirinas substituídas (CuP, CuOEP, CuTPP) diferem entre os dois métodos. A U-TDDFT prevê uma tendência de CuOEP < CuTPP < CuP, enquanto a X-TDDFT prevê CuTPP < CuOEP < CuP. Isso indica que os erros da U-TDDFT não são apenas grandes, mas também não sistemáticos, falhando em capturar tendências relativas mesmo para moléculas estreitamente relacionadas.

Significado e Afirmativas
O artigo afirma que o desenvolvimento de NACMEs em X-TDDFT completa o conjunto de ferramentas essencial para o estudo de processos fotofísicos e fotoquímicos de moléculas de camada aberta usando TDDFT adaptada ao spin. Os autores enfatizam que a adaptação ao spin não é meramente uma correção para energias e gradientes, mas é crucial para NACMEs e, consequentemente, para a precisão das dinâmicas não adiabáticas simuladas e das taxas de CI.

O estudo conclui que confiar em NACMEs de U-TDDFT para sistemas de camada aberta pode levar a previsões qualitativamente incorretas de mecanismos de relaxação e erros quantitativos em constantes de taxa que são tanto grandes quanto não sistemáticos. Portanto, os autores recomendam o uso de NACMEs em X-TDDFT para sistemas de camada aberta em geral, mesmo em casos onde as energias de excitação em U-TDDFT e X-TDDFT parecem semelhantes. O trabalho estabelece as bases para futuras aplicações em dinâmica molecular não adiabática e no estudo da fotofísica ultrarrápida em radicais e complexos de metais de transição.