Piezoaxial coupling for strain-selected ferroaxial domain control

Resumen Técnico: Acoplamiento Piezoaxial para el Control de Dominios Ferroaxiales Seleccionados por Deformación

Planteamiento del Problema
El orden ferroaxial es un estado ferroico caracterizado por un parámetro de orden que se transforma como un vector axial y es par tanto bajo la inversión espacial ($P$) como bajo la inversión temporal ($T$). Esta simetría lo distingue de los órdenes ferroeléctricos ($P$-impar) y ferromagnéticos ($T$-impar). En consecuencia, los campos eléctricos y magnéticos estándar no pueden servir como campos conjugados lineales universales para controlar el signo de un dominio ferroaxial, ya que no comparten la simetría de vector axial $P$-par y $T$-par requerida. Si bien se han establecido rutas dinámicas utilizando pulsos de terahercios con polarización circular para inducir la conmutación mediante campos axiales fuera del equilibrio, ha permanecido esquivo un campo conjugado estático y homogéneo capaz de sesgar los dominios ferroaxiales en equilibrio. Aunque han surgido firmas experimentales del acoplamiento deformación-ferroaxial (por ejemplo, en 1T-TiSe$_2$), no se había establecido un marco de simetría general que identificara los campos específicos derivados de la deformación para todos los grupos de puntos, junto con la verificación mediante primeros principios del resultante desdoblamiento de energía.

Metodología
Los autores desarrollan una jerarquía basada en la simetría de polinomios de deformación homogéneos estáticos que se transforman en la misma representación irreducible que un parámetro de orden ferroaxial seleccionado ($A$).

1. Formulación de Simetría: Los autores definen un término de acoplamiento piezoaxial en la energía libre, $F_A = -h(\boldsymbol{\varepsilon})A$, donde $h(\boldsymbol{\varepsilon})$ es un campo axial derivado de la deformación. Derivan sistemáticamente los polinomios de deformación de orden más bajo permitidos para grupos de puntos parentales no piroaxiales, clasificándolos por familia cristalina (ortorrómbica, tetragonal, trigonal, hexagonal y cúbica).
2. Verificación mediante Primeros Principios: Para validar el marco teórico, los autores realizan cálculos de la teoría del funcional de la densidad (DFT) en el compuesto ferroaxial trigonal Na$_2$BaMg(PO$_4$)$_2$. En este material, el parámetro de orden corresponde a la rotación colectiva de los tetraedros de PO$_4$ alrededor del eje trigonal.
3. Protocolos Computacionales: Se emplean tres cálculos complementarios:
   • Desdoblamiento de dominios con coordenadas bloqueadas: Cálculo de la diferencia de energía entre estados ferroaxiales pre-construidos con signos de rotación opuestos bajo deformación fija.
   • Escaneos de potencial de doble pozo: Mapeo del perfil de energía a lo largo de la coordenada de rotación para observar la inclinación del potencial.
   • Relajación con deformación fija: Relajación de las posiciones atómicas desde la estructura para-axial ($\phi=0$) bajo deformación fija para determinar si el campo de deformación selecciona por sí solo un pozo de dominio específico.

Contribuciones Clave y Resultados

• Jerarquía de Simetría de los Campos de Deformación: El artículo establece que el orden principal del campo piezoaxial depende del grupo de punto parental:
  • Ortorrómbico: Lineal en la deformación (campo de cizalla).
  • Tetragonal: Cuadrático en la deformación.
  • Trigonal y Hexagonal: Cúbico en la deformación.
  • Cúbico: Dividido en dos clases. Los grupos "Cúbico-I" (23, $\bar{m}3$) permiten campos de deformación total lineales para ejes seleccionados, mientras que los grupos "Cúbico-II" (432, $\bar{4}3m$, $m\bar{3}m$) prohíben los campos lineales, requiriendo combinaciones cuadráticas o cúbicas dependiendo del eje seleccionado ([001], [110] o [111]).
• Especificidades del Caso Trigonal: Para el sistema trigonal (por ejemplo, Na$_2$BaMg(PO$_4$)$_2$), el campo basal de orden principal es cúbico: $h(\boldsymbol{\varepsilon}) \propto Y(3X^2 - Y^2)$, donde $X$ e $Y$ son componentes de deformación desviatoria. En términos de amplitud de deformación con signo $\varepsilon_u$ y ángulo del eje principal $\theta$, esto resulta en $h \propto \varepsilon_u^3 \sin(6\theta)$.
• Verificación mediante Primeros Principios:
  • Dependencia Angular y de Amplitud: Los desdoblamientos de dominio calculados ($\Delta E_{dom}$) siguen estrictamente la dependencia angular $\sin(6\theta)$ predicha y la escala cúbica $\Delta E_{dom} \propto \varepsilon_u^3$.
  • Inclinación del Potencial: Los escaneos de energía confirman que el campo de deformación introduce un término impar en $\phi$ que sesga el potencial de doble pozo, eliminando la degeneración entre los dos dominios.
  • Relajación Seleccionada por Deformación: Partiendo de la estructura para-axial, las relajaciones con deformación fija evolucionan hacia mínimos ferroaxiales opuestos basándose únicamente en el signo de la deformación aplicada ($\varepsilon_u$) o la orientación del eje de deformación ($\theta$).
• Viabilidad Experimental: Los autores estiman que para un núcleo correlacionado de $\sim 100$ celdas unitarias, el sesgo inducido por la deformación con una deformación de $|\varepsilon_u| = 0.020$ puede volverse termodinámicamente significativo cerca de la temperatura de transición ($T_c \approx 540$ K), lo que sugiere la viabilidad para experimentos de enfriamiento mediante campo de deformación.

Significado
El artículo establece la deformación homogénea estática como un campo conjugado permitido por simetría para el orden ferroaxial. Proporciona un marco de simetría general para identificar los polinomios de deformación específicos requeridos para controlar los dominios ferroaxiales en diversos sistemas cristalinos. Al demostrar que la deformación puede seleccionar energéticamente los dominios ferroaxiales en equilibrio, el trabajo propone un protocolo de "enfriamiento por campo de deformación" análogo al enfriamiento por campo magnético. Esto ofrece una ruta para la ingeniería y el control de dominios ferroaxiales sin depender de campos dinámicos fuera del equilibrio. El marco se extiende más allá de las simples transiciones estructurales para incluir estados axiales electrónicos y de vector de onda finito (por ejemplo, estados CDW), siempre que el campo de deformación se entienda como un acoplamiento a un compuesto axial uniforme inducido. Los autores señalan que, si bien el marco es general, la realización experimental más directa se espera en transiciones ferroaxiales estructurales propias, como los fosfatos de tipo glaserita estudiados aquí.