Spin and Orbital Rashba response in ferroelectric polarized PtSe$_2$/MoSe$_2$/LiNbO$_3$ heterostructures

Resumo Técnico: Resposta de Rashba de Spin e Orbital em Heteroestruturas Ferroelétricas Polarizadas de PtSe2/MoSe2/LiNbO3

Enunciado do Problema
Avanços recentes em emissores de terahertz (THz) espintrônicos demonstraram que pulsos ópticos ultrarrápidos podem gerar radiação THz via conversão de spin-para-carga (SCC) em heteroestruturas ferromagnéticas/não-magnéticas (FM/NM). Embora mecanismos como o Efeito Hall de Spin Inverso (ISHE) e o Efeito Rashba-Edelstein Inverso (IREE) sejam bem estabelecidos, o papel do grau de liberdade orbital nesses processos permanece subexplorado na literatura, particularmente dentro de heteroestruturas de van der Waals (vdW). Além disso, embora substratos ferroelétricos como o LiNbO3 tenham sido mostrados como moduladores da emissão THz via campos elétricos dependentes da polarização, carece-se de uma quantificação por primeiros princípios de como esses campos influenciam tanto as texturas de Rashba de spin quanto de orbital — e a consequente eficiência de SCC. Este trabalho aborda a necessidade de quantificar o efeito Rashba em interfaces de PtSe2/MoSe2, incorporando a dinâmica orbital para explicar assinaturas experimentais de THz e a modulação da emissão pela polarização ferroelétrica.

Metodologia
Os autores empregaram uma abordagem de múltiplos estágios baseada em primeiros princípios, combinando Teoria do Funcional da Densidade (DFT) e teoria de resposta linear:

1. Relaxação Estrutural e Extração de Campo: Utilizando DFT-PBE dentro do código VASP (base de ondas planas, cutoff de 500 eV, correção vdW D2), os autores modelaram uma heteroestrutura de PtSe2/MoSe2/LiNbO3. Eles simularam duas configurações: camadas 2D posicionadas acima e abaixo de uma placa de LiNbO3 com corte em Z para representar direções opostas de polarização ($\hat{P}_{\uparrow}$ e $\hat{P}_{\downarrow}$). Ao comparar o potencial eletrostático da heteroestrutura completa com um bilayer de PtSe2/MoSe2 isolado, eles extraíram o campo elétrico efetivo fora do plano ($E_{eff}$) gerado pelo substrato ferroelétrico na interface 2D.
2. Cálculos de Transporte e Resposta: Utilizando o código SIESTA (conjuntos de base localizados DZP) e o formalismo de Kubo, os autores realizaram cálculos de transporte fora do equilíbrio. Um campo elétrico artificial, derivado dos valores extraídos de $E_{eff}$, foi aplicado a um sistema trilayer puro de PtSe2/MoSe2 para emular o ambiente ferroelétrico.
3. Análise de Resposta Linear: Os autores calcularam as acumulações de spin e orbital (densidades) induzidas por um campo elétrico no plano ($E_x$) usando a teoria de resposta linear. Isso envolveu a integração da fórmula de Kubo sobre uma malha de k-pontos de $200 \times 200$ para determinar as densidades de spin ($\langle \delta \hat{s}_y \rangle$) e orbital ($\langle \delta \hat{l}_y \rangle$) fora do equilíbrio. O estudo focou na janela de energia relevante para a excitação óptica experimental (aproximadamente -2 eV a +2 eV em relação ao nível de Fermi).

Principais Contribuições e Resultados

• Quantificação de Campos Elétricos Efetivos: O estudo quantificou com sucesso os campos elétricos efetivos atuando nas camadas 2D devido ao substrato de LiNbO3. Para a configuração com MoSe2 mais próximo da superfície de LiNbO3, o campo foi calculado em aproximadamente 0,115 V/Å, enquanto para a configuração com PtSe2 mais próximo, obteve-se ~0,047 V/Å. Esses valores foram usados para simular a modulação da polarização do sistema 2D.
• Texturas de Rashba Orbital: Os autores mapearam as texturas orbitais do sistema PtSe2/MoSe2. Observaram que, embora a quiralidade perto do ponto $\Gamma$ permanecesse relativamente estável, ocorreram distorções significativas perto dos pontos $K$ e $K'$ dependendo da direção da polarização. Notavelmente, a textura orbital exibiu uma forma mais hexagonal comparada à textura de spin, com quiralidade anti-horária em $K/K'$ e horária em $\Gamma$.
• Acumulação de Spin vs. Orbital:
  • Spin: As simulações revelaram uma grande janela de energia (aprox. -1,75 a -1,0 eV) onde a acumulação de spin é aumentada pela direção da polarização apontando para a interface ($\hat{P}_{\uparrow}$). Esse aumento alinha-se com observações experimentais de um aumento de ~30% na amplitude do campo elétrico THz.
  • Orbital: A acumulação orbital exibiu um sinal oposto à acumulação de spin dentro da janela [-1, 0] eV. Essa oposição é atribuída a uma correlação spin-órbita negativa ($\langle \hat{l} \cdot \hat{s} \rangle$) através da janela de energia.
• Modulação e Inversão de Sinal: O estudo identificou níveis de energia específicos onde os sinais de acumulação se invertem. Para o spin, uma mudança de sinal ocorre em torno de 0,9 eV, enquanto para a densidade orbital, ocorre em torno de -1,5 eV. Isso sugere que o sinal de THz poderia, teoricamente, ser ajustado alterando a energia do fóton do pulso de excitação ou o potencial químico via gating.
• Força de Acoplamento de Rashba: A constante de acoplamento de Rashba ($\alpha_R$) estimada para a interface PtSe2/MoSe2 é de aproximadamente -300 meV·Å, impulsionada principalmente pelos estados Se-p na interface. Este valor é superior ao encontrado em sistemas PtSe2/Grafeno (-195 meV·Å).
• Contribuição Orbital para a Emissão THz: Os autores calcularam a porcentagem de ganho no comprimento de Rashba-Edelstein inverso ($\Lambda_{REE}$) entre os estados de polarização. Enquanto a contribuição de spin mostrou um ganho de <30%, a contribuição orbital demonstrou uma sensibilidade significativamente maior, atingindo um aumento de ~100% em janelas de energia específicas. Isso indica que o grau de liberdade orbital é altamente suscetível à polarização ferroelétrica.

Significância e Alegações
O artigo afirma que a incorporação do grau de liberdade orbital é crucial para uma compreensão mais profunda da conversão de momento angular-para-carga e da emissão THz em heteroestruturas de vdW. Os autores asseveram que suas simulações de primeiros princípios fornecem uma explicação quantitativa para os resultados experimentais de THz-TDS, especificamente a modulação da emissão pela polarização ferroelétrica.

As principais alegações incluem:

1. Mecanismos Complementares: As respostas de spin e orbital, embora possuam sinais opostos em certos intervalos de energia devido à correlação spin-órbita negativa, equilibram-se para produzir uma acumulação líquida positiva que se alinha com os sinais experimentais de THz.
2. Sensibilidade Orbital: A contribuição orbital para o efeito Rashba-Edelstein mostrou-se mais sensível aos campos de polarização do que a contribuição de spin, oferecendo um caminho potencial para aumentar ou ajustar a eficiência da emissão THz.
3. Estrutura Metodológica: O trabalho estabelece uma metodologia para calcular densidades de spin e orbital fora do equilíbrio na presença de campos ferroelétricos, que pode ser adaptada para projetar propriedades de materiais em heteroestruturas de vdW.
4. Direções Futuras: Os autores observam modestamente que, embora a resposta impulsionada eletricamente seja quantificada, a quantidade exata de polarização orbital injetada do lado ferromagnético permanece uma quantidade desconhecida, dependente dos detalhes da interface. Eles sugerem que a engenharia de materiais futura (por exemplo, usando elementos como Ni) poderia maximizar a contribuição orbital, potencialmente levando a amostras espintrônicas onde a resposta orbital é predominante.

O estudo conclui que compreender a dinâmica de ambos os graus de liberdade, spin e orbital, é essencial para o desenvolvimento de tecnologias THz de próxima geração e dispositivos de torque de spin-órbita, particularmente em materiais centrosimétricos como o PtSe2, onde o levantamento da degenerescência de spin pode gerar eficiências significativas.