Symmetry-Based Design Rules for Second-Harmonic Generation in Stacked and Twisted MoS2 Bilayers

Resumo Técnico: Regras de Design Baseadas em Simetria para Geração de Segundo Harmônico em Bilayer de MoS$_2$ Empilhados e Torcidos

Problema
A resposta óptica não linear de dicalcogenetos de metais de transição (TMDs) bidimensionais, particularmente o MoS$_2$, é altamente sensível às configurações de empilhamento e à rotação interlaminar. Embora cálculos de primeiros princípios possam modelar sistemas idealizados, eles se tornam computacionalmente proibitivos para configurações de baixa simetria realistas e superredes de moiré de grande escala. Há uma necessidade de um arcabouço rigoroso, baseado em simetria, que possa mapear configurações estruturais para seus tensores de suscetibilidade de segunda ordem ($\chi^{(2)}$) sem depender de cálculos de força bruta para cada ângulo de torção ou variante de empilhamento específica.

Metodologia
Os autores estabelecem um arcabouço abrangente combinando análise de teoria de grupos com validações direcionadas de primeiros princípios:

1. Análise de Simetria: Usando o princípio de Neumann, os autores derivam os componentes permitidos e proibidos do tensor $\chi^{(2)}$ para quatro grupos de pontos distintos relevantes para bilayers de MoS$_2$: $D_{3h}$ (monocamada e empilhamento AA), $D_{3d}$ (empilhamento AB(2H)), $C_{3v}$ (empilhamento AB(3R)) e $C_3$ (bilayers torcidos). Eles analisam como elementos de simetria específicos (centro de inversão $I$, espelho horizontal $\sigma_h$, espelhos verticais $\sigma_v$) restringem a estrutura do tensor.
2. Cálculos de Primeiros Princípios: As previsões teóricas são validadas contra cálculos de teoria do funcional da densidade (DFT) usando o código VASP com o funcional PBE e correções DFT-D3 para forças de van der Waals. A resposta óptica não linear é computada usando funções de Wannier maximamente localizadas (via Wannier90 e postw90) dentro da aproximação de partícula independente. Os cálculos incluem acoplamento spin-órbita e decompõem a resposta em contribuições interbanda e intrabanda.
3. Modelagem de Resposta Angular: Os autores derivam expressões analíticas para a intensidade de geração de segundo harmônico (SHG) copolarizada $I_{\parallel}(\theta)$ como uma função do ângulo de polarização incidente $\theta$ para cada classe de simetria, investigando especificamente o surgimento de deslocamentos de fase em configurações torcidas.

Principais Contribuições

• Mapeamento Simetria-Resposta: O artigo fornece um mapeamento definitivo de como as configurações de empilhamento (AA, AB(2H), AB(3R) e torcido) ditam a estrutura do tensor $\chi^{(2)}$ através dos grupos de pontos $D_{3h}$, $D_{3d}$, $C_{3v}$ e $C_3$.
• Regras de Design para Supressão e Ativação de SHG:
  • Simetria de Inversão ($D_{3d}$): A presença de um centro de inversão no empilhamento AB(2H) impõe estritamente $\chi^{(2)} = 0$, levando à supressão completa de SHG via interferência destrutiva interlaminar.
  • Quebra de Espelho Horizontal ($C_{3v}$ e $C_3$): A quebra da simetria $\sigma_h$ em bilayers AB(3R) e torcidos ativa componentes de tensor fora do plano ($\chi_{zzz}^{(2)}, \chi_{zxx}^{(2)}, \chi_{xxz}^{(2)}$), que são proibidas em sistemas $D_{3h}$.
  • Quebra de Espelho Vertical ($C_3$): A remoção de planos de espelho vertical em bilayers torcidos ativa o componente no plano $\chi_{xxx}^{(2)}$, que é estritamente zero em configurações de maior simetria.
• Bloqueio Geométrico dos Lobos de SHG: Os autores demonstram que a ativação de $\chi_{xxx}^{(2)}$ em bilayers torcidos induz uma rotação azimutal rígida dos lobos polares de SHG. Crucialmente, eles descobrem que o ângulo de rotação $\phi_0$ é independente da frequência em todo o espectro óptico e está geometricamente travado ao ângulo de torção estrutural $\theta_{tw}$ via a relação $\theta_{tw} = 2\phi_0$.

Resultados

• Validação do Tensor: Cálculos de primeiros princípios confirmam as previsões da teoria de grupos. O bilayer AB(2H) mostra sinais de SHG evanescentes. O bilayer AB(3R) recupera a resposta no plano semelhante à monocamada, mas ganha componentes significativas fora do plano. O bilayer torcido ($C_3$) exibe simultaneamente $\chi_{yyy}^{(2)}$ e $\chi_{xxx}^{(2)}$ não nulos.
• Características Espectrais: Os espectros de SHG são dominados por ressonâncias de dois fótons correspondentes às transições excitônicas A e B (bandas de valência com divisão por spin-órbita) e um pico C de maior energia associado a transições de band-nesting.
• Rotação do Padrão Polar: Para um bilayer torcido com $\theta_{tw} = 21,8^\circ$, o padrão polar de SHG rotaciona em $\phi_0 \approx 10,9^\circ$. Este deslocamento permanece constante em toda a faixa de frequência computada, confirmando que a relação de fase entre $\chi_{xxx}^{(2)}$ e $\chi_{yyy}^{(2)}$ é ditada pela geometria estrutural, e não por variações espectrais dinâmicas.
• Componentes Quirais: Embora o componente quiral fora do plano $\chi_{xyz}^{(2)}$ seja permitido pela simetria em $C_3$, sua magnitude é encontrada como excepcionalmente pequena devido ao fraco acoplamento interlaminar e à preservação local da simetria de espelho dentro das monocamadas individuais.

Significância
O artigo estabelece uma rota óptica direta e independente do comprimento de onda para determinar o ângulo de torção em materiais 2D. Ao vincular o ângulo de torção estrutural diretamente à rotação dos lobos polares de SHG ($\theta_{tw} = 2\phi_0$), os autores fornecem uma ferramenta de caracterização não destrutiva e de alto rendimento para superredes de moiré e domínios de empilhamento. Além disso, o trabalho oferece um roteiro claro para a engenharia de respostas ópticas não lineares em materiais 2D puramente através do controle de empilhamento, permitindo o design de heteroestruturas com funcionalidades customizadas, como supressão completa de SHG ou ativação de componentes de tensor específicos. O arcabouço é apresentado como universal para bilayers de 2H-TMD e cristais de van der Waals análogos, fornecendo uma base fundamental para futuras triagens de alto rendimento de arquiteturas ópticas não lineares.