Piezoaxial coupling for strain-selected ferroaxial domain control

Résumé technique : Couplage piézoaxial pour le contrôle des domaines ferroaxiaux par sélection de déformation

Énoncé du problème
L'ordre ferroaxial est un état ferroïque caractérisé par un paramètre d'ordre qui se transforme comme un vecteur axial et qui est pair sous l'inversion spatiale ($P$) et la renversement du temps ($T$). Cette symétrie le distingue des ordres ferroélectriques ($P$-impair) et ferromagnétiques ($T$-impair). Par conséquent, les champs électriques et magnétiques standards ne peuvent pas servir de champs conjugués linéaires universels pour contrôler le signe d'un domaine ferroaxial, car ils ne partagent pas la symétrie de vecteur axial $P$-pair et $T$-pair requise. Bien que des voies dynamiques utilisant des impulsions térahertz à polarisation circulaire aient été établies pour induire une commutation via des champs axiaux hors équilibre, un champ conjugué statique et homogène capable de biaiser les domaines ferroaxiaux en équilibre est resté insaisissable. Bien que des signatures expérimentales du couplage déformation-ferroaxial (par exemple, dans 1T-TiSe$_2$) soient apparues, un cadre de symétrie général identifiant les champs spécifiques dérivés de la déformation pour tous les groupes ponctuels, couplé à une vérification par premiers principes de l'éclatement d'énergie résultant, n'avait pas encore été établi.

Méthodologie
Les auteurs développent une hiérarchie de polynômes de déformation homogènes basés sur la symétrie qui se transforment selon la même représentation irréductible qu'un paramètre d'ordre ferroaxial sélectionné ($A$).

1. Formulation de la symétrie : Les auteurs définissent un terme de couplage piézoaxial dans l'énergie libre, $F_A = -h(\boldsymbol{\varepsilon})A$, où $h(\boldsymbol{\varepsilon})$ est un champ axial dérivé de la déformation. Ils dérivent systématiquement les polynômes de déformation de plus bas ordre autorisés pour les groupes ponctuels parents non pyroaxiaux, en les classant par famille cristalline (orthorhombique, tétragonale, trigonale, hexagonale et cubique).
2. Vérification par premiers principes : Pour valider le cadre théorique, les auteurs effectuent des calculs de la théorie de la fonctionnelle de la densité (DFT) sur le composé ferroaxial trigonal Na$_2$BaMg(PO$_4$)$_2$. Dans ce matériau, le paramètre d'ordre correspond à la rotation collective des tétraèdres PO$_4$ autour de l'axe trigonal.
3. Protocoles de calcul : Trois calculs complémentaires sont employés :
   • Éclatement de domaine en coordonnées contraintes : Calcul de la différence d'énergie entre des états ferroaxiaux préconstruits avec des signes de rotation opposés sous déformation fixe.
   • Scans de potentiel à double puits : Cartographie du profil d'énergie le long de la coordonnée de rotation pour observer l'inclinaison du potentiel.
   • Relaxation à déformation fixe : Relaxation des positions atomiques à partir de la structure para-axiale ($\phi=0$) sous déformation fixe pour déterminer si le champ de déformation sélectionne à lui seul un bassin de domaine spécifique.

Contributions clés et résultats

• Hiérarchie de symétrie des champs de déformation : L'article établit que l'ordre principal du champ piézoaxial dépend du groupe ponctuel parent :
  • Orthorhombique : Linéaire en déformation (champ de cisaillement).
  • Tétragonal : Quadratique en déformation.
  • Trigonal et Hexagonal : Cubique en déformation.
  • Cubique : Divisé en deux classes. Les groupes « Cubique-I » (23, $\bar{m}3$) permettent des champs de déformation totale linéaires pour des axes sélectionnés, tandis que les groupes « Cubique-II » (432, $\bar{4}3m$, $m\bar{3}m$) interdisent les champs linéaires, nécessitant des combinaisons quadratiques ou cubiques selon l'axe sélectionné ([001], [110] ou [111]).
• Spécificités du cas trigonal : Pour le système trigonal (ex. Na$_2$BaMg(PO$_4$)$_2$), le champ basal principal est cubique : $h(\boldsymbol{\varepsilon}) \propto Y(3X^2 - Y^2)$, où $X$ et $Y$ sont des composantes de déformation déviatorique. En termes d'amplitude de déformation signée $\varepsilon_u$ et d'angle $\theta$ de l'axe principal, cela donne $h \propto \varepsilon_u^3 \sin(6\theta)$.
• Vérification par premiers principes :
  • Dépendance angulaire et d'amplitude : Les éclatements de domaine calculés ($\Delta E_{dom}$) suivent strictement la dépendance angulaire $\sin(6\theta)$ prédite et l'échelle cubique $\Delta E_{dom} \propto \varepsilon_u^3$.
  • Inclinaison du potentiel : Les scans d'énergie confirment que le champ de déformation introduit un terme impair en $\phi$ qui biaise le potentiel à double puits, levant la dégénérescence entre les deux domaines.
  • Relaxation sélectionnée par la déformation : En partant de la structure para-axiale, les relaxations sous déformation fixe évoluent vers des minima ferroaxiaux opposés uniquement en fonction du signe de la déformation appliquée ($\varepsilon_u$) ou de l'orientation de l'axe de déformation ($\theta$).
• Faisabilité expérimentale : Les auteurs estiment que pour un noyau corrélé d'environ 100 cellules unitaires, le biais induit par la déformation à une déformation de $|\varepsilon_u| = 0,020$ peut devenir thermodynamiquement significatif près de la température de transition ($T_c \approx 540$ K), suggérant la viabilité d'expériences de refroidissement par champ de déformation.

Signification
Cet article établit la déformation homogène statique comme un champ conjugué permis par la symétrie pour l'ordre ferroaxial. Il fournit un cadre de symétrie général pour identifier les polynômes de déformation spécifiques requis pour contrôler les domaines ferroaxiaux à travers divers systèmes cristallins. En démontrant que la déformation peut sélectionner énergétiquement les domaines ferroaxiaux en équilibre, ce travail propose un protocole de « refroidissement par champ de déformation » analogue au refroidissement par champ magnétique. Cela offre une voie pour l'ingénierie et le contrôle des domaines ferroaxiaux sans dépendre de champs dynamiques hors équilibre. Le cadre s'étend au-delà des simples transitions structurelles pour inclure des états axiaux électroniques et à vecteur d'onde fini (ex. états CDW), à condition que le champ de déformation soit compris comme un couplage à un composite axial uniforme induit. Les auteurs notent que bien que le cadre soit général, la réalisation expérimentale la plus directe est attendue dans les transitions ferroaxiales structurelles propres, telles que les phosphates de type glaserite étudiés ici.