Piezoaxial coupling for strain-selected ferroaxial domain control

Resumo Técnico: Acoplamento Piezoaxial para Controle de Domínios Ferroaxiais Selecionados por Deformação

Definição do Problema
A ordem ferroaxial é um estado ferroico caracterizado por um parâmetro de ordem que se transforma como um vetor axial e é par sob tanto a inversão espacial ($P$) quanto a reversão temporal ($T$). Esta simetria a distingue das ordens ferroelétricas ($P$-ímpar) e ferromagnéticas ($T$-ímpar). Consequentemente, campos elétricos e magnéticos padrão não podem servir como campos conjugados lineares universais para controlar o sinal de um domínio ferroaxial, pois eles não compartilham a simetria de vetor axial $P$-par e $T$-par exigida. Embora rotas dinâmicas usando pulsos de terahertz com polarização circular tenham sido estabelecidas para induzir a comutação via campos axiais fora do equilíbrio, um campo conjugado estático e homogêneo capaz de enviesar domínios ferroaxiais em equilíbrio permaneceu elusivo. Embora assinaturas experimentais do acoplamento deformação-ferroaxial (ex: em 1T-TiSe$_2$) tenham emergido, um arcabouço de simetria geral identificando os campos específicos derivados da deformação para todos os grupos pontuais, juntamente com a verificação por primeiros princípios do resultante desdobramento de energia, não havia sido estabelecido.

Metodologia
Os autores desenvolvem uma hierarquia baseada em simetria de polinômios de deformação homogênea estática que se transformam na mesma representação irredutível de um parâmetro de ordem ferroaxial selecionado ($A$).

1. Formulação de Simetria: Os autores definem um termo de acoplamento piezoaxial na energia livre, $F_A = -h(\boldsymbol{\varepsilon})A$, onde $h(\boldsymbol{\varepsilon})$ é um campo axial derivado da deformação. Eles derivam sistematicamente os polinômios de deformação de ordem mais baixa permitidos para grupos pontuais parentais não-piroaxiais, classificando-os por família cristalina (ortorrômbica, tetragonal, trigonal, hexagonal e cúbica).
2. Verificação por Primeiros Princípios: Para validar o arcabouço teórico, os autores realizam cálculos de teoria do funcional da densidade (DFT) no composto ferroaxial trigonal Na$_2$BaMg(PO$_4$)$_2$. Neste material, o parâmetro de ordem corresponde à rotação coletiva dos tetraedros de PO$_4$ em torno do eixo trigonal.
3. Protocolos Computacionais: Três cálculos complementares são empregados:
   • Desdobramento de domínio por coordenadas presas: Cálculo da diferença de energia entre estados ferroaxiais pré-construídos com sinais de rotação opostos sob deformação fixa.
   • Varreduras de potencial de poço duplo: Mapeamento do perfil de energia ao longo da coordenada de rotação para observar a inclinação do potencial.
   • Relaxação de deformação fixa: Relaxação das posições atômicas a partir da estrutura para-axial ($\phi=0$) sob deformação fixa para determinar se o campo de deformação sozinho seleciona um poço de domínio específico.

Principais Contribuições e Resultados

• Hierarquia de Simetria dos Campos de Deformação: O artigo estabelece que a ordem principal do campo piezoaxial depende do grupo pontual parental:
  • Ortorrômbica: Linear na deformação (campo de cisalhamento).
  • Tetragonal: Quadrática na deformação.
  • Trigonal e Hexagonal: Cúbica na deformação.
  • Cúbica: Dividida em duas classes. Grupos "Cúbicos-I" (23, $\bar{m}3$) permitem campos de deformação total lineares para eixos selecionados, enquanto grupos "Cúbicos-II" (432, $\bar{4}3m$, $m\bar{3}m$) proíbem campos lineares, exigindo combinações quadráticas ou cúbicas dependendo do eixo selecionado ([001], [110] ou [111]).
• Especificidades do Caso Trigonal: Para o sistema trigonal (ex: Na$_2$BaMg(PO$_4$)$_2$), o campo basal principal é cúbico: $h(\boldsymbol{\varepsilon}) \propto Y(3X^2 - Y^2)$, onde $X$ e $Y$ são componentes de deformação desviadora. Em termos de amplitude de deformação assinada $\varepsilon_u$ e ângulo do eixo principal $\theta$, isso resulta em $h \propto \varepsilon_u^3 \sin(6\theta)$.
• Verificação por Primeiros Princípios:
  • Dependência Angular e de Amplitude: Os desdobramentos de domínio calculados ($\Delta E_{dom}$) seguem estritamente a dependência angular $\sin(6\theta)$ prevista e a escala cúbica $\Delta E_{dom} \propto \varepsilon_u^3$.
  • Inclinação do Potencial: Varreduras de energia confirmam que o campo de deformação introduz um termo ímpar em-$\phi$ que enviesa o potencial de poço duplo, levantando a degenerescência entre os dois domínios.
  • Relaxação Selecionada por Deformação: Partindo da estrutura para-axial, as relaxações de deformação fixa evoluem para mínimos ferroaxiais opostos baseando-se unicamente no sinal da deformação aplicada ($\varepsilon_u$) ou na orientação do eixo de deformação ($\theta$).
• Viabilidade Experimental: Os autores estimam que, para um núcleo correlacionado de $\sim 100$ células unitárias, o viés induzido pela deformação em uma deformação de $|\varepsilon_u| = 0.020$ pode tornar-se termodinamicamente significativo próximo à temperatura de transição ($T_c \approx 540$ K), sugerindo viabilidade para experimentos de resfriamento por campo de deformação.

Significância
O artigo estabelece a deformação homogênea estática como um campo conjugado permitido por simetria para a ordem ferroaxial. Ele fornece um arcabouço de simetria geral para identificar os polinômios de deformação específicos necessários para controlar domínios ferroaxiais através de vários sistemas cristalinos. Ao demonstrar que a deformação pode selecionar energeticamente domínios ferroaxiais em equilíbrio, o trabalho propõe um protocolo de "resfriamento por campo de deformação" análogo ao resfriamento por campo magnético. Isso oferece uma rota para a engenharia e controle de domínios ferroaxiais sem depender de campos dinâmicos fora do equilíbrio. O arcabouço estende-se além de transições estruturais simples para incluir estados axiais eletrônicos e de vetor de onda finito (ex: estados CDW), desde que o campo de deformação seja entendido como um acoplamento a um composto axial uniforme induzido. Os autores observam que, embora o arcabouço seja geral, a realização experimental mais direta é esperada em transições ferroaxiais estruturais próprias, como os fosfatos do tipo glaserite estudados aqui.