Thesis 2
Thesis 2
A N D M U LT I F E R R O I C M AT E R I A L S B Y
F I R S T - P R I N C I P L E S C A L C U L AT I O N S
—-
Faculté des Sciences
Université de Liège
&
Centro de Investigación y de Estudios Avanzados IPN
Abril ����
S T U D Y A N D C H A R A C T E R I Z AT I O N O F M A G N E T I C
A N D M U LT I F E R R O I C M AT E R I A L S B Y
F I R S T - P R I N C I P L E S C A L C U L AT I O N S
Supervisors:
Prof. Aldo Humberto Romero Castro, CINVESTAV-Querétaro
Dr. Eric Bousquet, Université de Liège (Secretary)
Jury Members:
Prof. Nicola Spaldin, ETH-Zürich
Prof. Philippe Ghosez, Université de Liège (President)
Prof. Sergio Jiménez Sandoval, CINVESTAV-Querétaro
Prof. Rebeca Castanedo Pérez, CINVESTAV-Querétaro
Prof. Arturo Mendoza Galván, CINVESTAV-Querétaro
—-
Ph.D-Candidate
Physics Department
Université de Liège
&
Centro de Investigación y de Estudios Avanzados IPN
Abril ����
Andrés Camilo García Castro: Study and Characterization of Magnetic and Mul-
tiferroic Materials by First-principles Calculations, —-, c Abril ����
Dedicated to my father and my grandfather memory that always encouraged me
to follow my dreams and be consistent with it.
In the last fifteen years, multifunctional materials, and more specifically, multi-
functional oxides have been widely studied due to its wide range of properties.
Properties that go from superconductivity to ferroelectricity passing through mag-
netism and multiferroism have been reported. Nonetheless, the fluoride family was
left aside and little information is known about its possible ferroelectricity or mul-
tiferroism. In this Ph.D thesis, we explored the electronic, vibrational, structural
and magnetic properties of fluoride perovskite-based compounds. To such pur-
poses, We performed ab-initio calculations based in the density-functional theory
(DFT) as implemented in VASP and CRYSTAL codes.
Our first step was to perform vibrational analyses in a large set of fluoroper-
ovskites ABF3 . Based on the results, we proposed a model that establishes an
A-site geometrically driven ferroelectric vibrational instability in fluorides. Our
studies reveal a different behavior as a function of isotropic pressure for NaBF3
with respect to oxides (e.g. BaTiO3 ) with B = Ca, V, Mn, and Zn. For these
compounds we found an increase of the ferroelectric instability as a function of
hydrostatic pressure. This probably due to the “transformation” of eigendisplace-
ments responsible for the mode that creates the corresponding instability. In
particular, an increase of ionic A-site radii present a strong influence in FE-polar
instability.
vii
Magnetic ordering reveals a non-collinear ground state with the Gz Ax Fy repre-
sentation. Even more interesting, and non-linear magnetoelectric coupling was
found under the strained Pna�1 ground state becoming the first known multifer-
roic/magnetoelectric perovskite fluoride.
viii
RÉSUMÉ
Durant les quinze dernières années, les matériaux multifonctionnels et plus par-
ticulièrement les matériaux multifonctionnels à base d’oxyde ont été extrême-
ment étudiés du fait de leurs propriétés très diverses. Ces propriétés vont de la
supraconductivité à la ferroélectricité en passant par le magnétisme ou le multi-
ferroïsme. Dans cet engouement pour les oxydes, la famille des fluorites a cepen-
dant été délaissée de telle manière que peu d’informations sont connues quant
à leur propension à être ferroélectrique ou multiférroïque. Dans cette thèse de
doctorat, nous nous sommes proposés d’explorer les propriétés électroniques, vi-
brationnelles, structurales et magnétiques des matériaux perovskites à base de
fluor. Cette étude repose principalement sur des calculs ab-initio basés sur la
théorie de la fonctionnelle de la densité. Les codes de calcul utilisés sont princi-
palement VASP et CRYSTAL.
Dans un premier temps nous avons analysé les propriétés vibrationnelles d’une
large gamme de fluoro-perovskites ABF3 . Cette première étude nous a permis
de montrer que ces cristaux présentent une instabilité ferroélectrique dans leur
phase cubique qui est d’origine géométrique et qui est dominée par le cation A.
Cette origine géométrique rend la dépendance de l’instabilité ferroélectrique à
une pression hydrostatique très différente de celle observée dans les oxydes. Dans
les fluorites la pression a tendance à augmenter l’instabilité ferroélectrique tandis
que dans les oxydes elle a tendance à la supprimer.
En nous basant sur une analyse des symétries et sur des calculs ab-initio, nous
avons aussi montré que les phases perovskites Pnma, observées à basse pression,
et les phases post-pérovskites, observées à haute pression, contenant un cation
B magnétique ont toutes une structure magnétique non-colinéaire. Cette étude a
permis de clarifier les propriétés magnétiques observées expérimentalement et de
prédire une transition Pnma/post-perovskite dans le NaMnF3 sous haute pres-
sion avec un ferromagnétisme léger bien plus important que dans tous les autres
systèmes reportés jusqu’à présent.
ix
réponse inhabituelle à une modification du vecteur de déplacement ferroélectrique
des atomes où sous tension la polarisation est dominée par un déplacement
géométrique du cation A et sous compression par le déplacement des cations
A et Mn. L’analyse du magnétisme montre un ordre antiferromagnétique domi-
nant de type G associé à des ordres non-colinéaires de type A et F (Gz Ax Fy ).
L’association de l’ordre ferroélectrique et du ferromagnétisme léger (Fy ) fait que
NaMnF3 sous contrainte d’épitaxie est un nouveau multiferroïque. L’étude de
la réponse magnétoélectrique sous champ magnétique montre qu’elle s’effectue
également à l’ordre non-linéaire.
Enfin, dans un dernier temps nous avons aussi étudié la structure électronique
et les propriétés structurales d’heterostructures oxyfluorites (KTaO3 )n /(KBF3 )l
où B = Zn et Ni. Dans ces interfaces entre un oxyde et une fluorite nous avons
observé que la levée de dégénérescence des orbitales atomiques est fortement mod-
ifiée par la liaison O–B–F. Un phénomène de catastrophe polaire est aussi observé
où l”epaisseur critique pour l’apparition d’un gaz bidimensionnel est plus petite
que celle reportée dans les super-réseaux SrTiO3 /LaAlO3 . Dans le cas B = Ni, le
magnétisme dans la couche de KNiF3 reste globalement du type G mais présente
une modulation en amplitude des moments magnétiques près de l’interface. De
manière surprenante, nous avons aussi observé une levée de dégénérescence de
type k3 -Rashba à l’interface de ces heterostructures oxyfluorites et d’amplitude
au moins quatre fois supérieure à celle rapportée dans l’interface SrTiO3 /LaAlO3
et deux fois supérieur à celle des transistors à base de KTaO3 .
En résumé, notre étude a permis de conclure que les cristaux perovskites à base
de fluorites sont de bons candidats à considérer dans la recherche de propriétés
multifonctionnelles. Nous espérons que nos résultats motiveront des expérimenta-
teurs et d’autres théoriciens dans l’étude et la caractérisation des fluorites dans
la recherche de nouveaux matériaux aux propriétés inédites et à haut potentiel
pour des applications technologiques et industrielles.
x
RESÚMEN
También se ha demostrado, que todos las post-perovskitas ABX3 con una sitio
B magnéticamente activo exhiben la posibilidad de presentar un comportamiento
magnético no colineal permitido por la simetría. Lo anterior con base en los cálcu-
los de primeros principios y el análisis de teoría de grupos. Con estos resultados
fueron aclaradas las propiedades magnéticas encontrados experimentalmente en
esta fase perovksita de alta presión encontrada principalmente en el manto de la
Tierra. Además, fue posible predecir la existencia de esta transición a fase post-
perovskita en el compuesto NaMnF3 . Con esto, este fluoruro se convertirse en un
nuevo miembro de esta familia en el cual, las propiedades magnéticas, y sobre
todo, su ferromagnetismo no colineal muestra un alta componente de momento
magnético.
xi
perovskita. Además, se observó una polarización fuera de plano en la tendencia
opuesta a lo encontrado en óxidos. Se encontró ademas un ordenamiento ferroeléc-
trico cooperativo entre los sitios de Na y Mn para el esfuerzo de compresión. Lo
anterior en contraste con una ferroelectricidad geométricamente dominada por el
sitio A a valores de tensión en el plano ortorrómbico xz. El ordenamiento mag-
nético revela un estado fundamental no colineal bajo la representación Gz Ax Fy .
Aún más interesante, se observó un acople magnetoeléctrico no lineal en este ma-
terial bajo la simetría Pna�1 de estado fundamental, convirtiéndose con esto en
el primer multiferroico/magnetoeléctrico conocido en estructura perovskita base
fluoruro.
xii
P U B L I C AT I O N S
Some ideas and figures have appeared previously in the following publications:
xiii
[�] Sobhit Singh, A. C. Garcia-Castro, Irais Valencia-Jaime, Francisco Munoz,
and Aldo H. Romero, Prediction of a controllable Weyl semi-metallic phase in
inversion-asymmetryc BiSb, Advanced Materials, (���� - To be submitted).
xiv
“Because success is stumbling from failure to failure with no loss of enthusiasm”
—Wiston S. Churchill.
ACKNOWLEDGEMENTS
Sometimes, you face the difficult decision to leave your country and your home in
the search for better opportunities and to explore new horizons. This path that
I traced in my life � years ago now is rewarded with this amazing achievement
not just for my professional life, but also for my personal formation as well. Thus,
the culmination of the last � year of this journey, full of nice and not that good
experiences, is condensed in the present manuscript. Therefore, I’m taking this
chance, and a few words, to acknowledge all this people that were and have been
there for me and support me along this journey, because, I have always believed
that when you achieve something worthy in your life, this is the sum of all the
efforts of the people around you.
First, I would like to thank to prof. Aldo H. Romero for his great supervision
and friendship that always extended to me kindly, and for bring me the oppor-
tunity to join his research group. I also want to thank him for encouraging and
support my scientific travels used to interact with the scientific community.
I would also equally thank to Dr. Eric Bousquet for all his help, guidance and
friendly patience to explain me a lot of scientific aspects, but more important
than that, for being there to help me even in small things such as finding an
apartment to rent.
I am very grateful to my fiancée, Sara Orozco, who was always there for me
listening my scientific explanations and ideas, even knowing that she is not a
scientist. Who has been there for me and had understood when I have to travel
for a long time and leave her alone.
To the Euraxess fellowship and its stuff for their friendly help and support that
allowed me to be in Liège for more than a year and enjoy the delicious chocolates,
xv
nice beers and tasty food.
I thank to prof. Nicola A. Spaldin and her research group at the ETH-Zürich
for the chance to work with them and have the opportunity to be in the beautiful
Zürich for several months.
To my friend and colleague Irais Valencia for her support and patience to hear
my ideas even when they did not take me to a successful results.
It is worth to mention that all the computational work presented here was pos-
sible thanks to support —for more than � millions of computational hours— from
the Extreme Science and Engineering Discovery Environment, XSEDE at United
States with the supercomputers Stampede and Ranger (project TG-DMR������).
Thanks also to the computational time provided by the Belgian CECI projects
and specially the Zenobe stuff.
At last, but not least, I would like to thanks to my parents, Alberto and Martha,
for teaching me the value of the honesty, hard work, and respect. Thanks to my
sister Cindy Johana, that is constantly in touch with me giving me the courage
to go forward. I thank to my friends Luis Alberto, Sebastian, and Catalina that
were always worried about my thesis progress and were continuously asking: “how
is the thesis going?”. Finally and specially, thanks to Wilfredo and his wife Rosa
Elisa, that opened their home for me and were a friendly face every time that I
traveled from México to Belgium.
xvi
CONTENTS
Acronyms xx
Introduction xxi
xvii
xviii contents
�.� Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ��
� non-collinear magnetism in c a r h o 3 and n a n i f 3 post-
perovskites ��
�.� Computational Details . . . . . . . . . . . . . . . . . . . . . . . . . ��
�.� Structural Characterization of CaRhO 3 and NaNiF 3 . . . . . . ��
�.� Non-collinear Ordering . . . . . . . . . . . . . . . . . . . . . . . . ��
�.� Exploring the Magnetic Interactions . . . . . . . . . . . . . . . . . ��
�.� Summary and Conclusions . . . . . . . . . . . . . . . . . . . . . . ��
� induced multiferroic and magnetoelectric behav-
ior in n a m n f 3 ��
�.� Computational Details . . . . . . . . . . . . . . . . . . . . . . . . . ��
�.� Epitaxial strain in Pnma NaBF 3 Fluorides . . . . . . . . . . . . . ��
�.� Strain-Induced Ferroelectricity in NaMnF 3 . . . . . . . . . . . . ��
�.�.� Second-order piezoelectricity in Pna�1 NaMnF3 . . . . . . ��
�.� Non-Collinear magnetism and ME coupling in Pna�1 NaMnF 3 . ��
�.�.� Field-induced distortions as the source of the ME coupling ��
�.�.� Effect of the U and J into the ME coupling in Pna�1 NaMnF3 ��
�.�.� Special remark: some ferroelectrics under strain . . . . . . ��
�.� Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ��
v appendix ���
a overview of the magnetic interactions in fluoroper-
ovskites ���
b group theory analysis to the non-collinear mag-
netism in ppv ���
c b a c o f 4 thin films ���
c.� Computational Details . . . . . . . . . . . . . . . . . . . . . . . . . ���
c.� Strain Effect in B a C o F4 Fluoride: . . . . . . . . . . . . . . . . . . ���
c.� Noncollinear Magnetism and Multiferroic Behavior . . . . . . . . . ���
contents xix
references ���
ACRONYMS
AFD Antiferrodistortive
AFM Antiferromagnetism
A-AFM A-type Antiferromagnetism
BZ Brillouin Zone
C-AFM C-type Antiferromagnetism
DFPT Density-Functional Perturbation Theory
DFT Density-Functional Theory
FE Ferroelectric
FM Ferromagnetism
G-AFM G-type Antiferromagnetism
GGA Generalized Gradient Approximation
HSE�� Heyd-Scuseria-Ernzerhof
IMT Insulator-to-Metal Transition
LDA Local Density Approximation
LDOS Local-Density Of States
LSDA Local Spin Density Approximation
ME Magnetoelectric
MF Multiferroic
PAW Projected Augmented-Wave
PBE Perdew Burke Erznzerhof
PBEsol Perdew Burke Erznzerhof Revised for Solids
Pv Perovskite
pPv Post-perovskite
SOC Spin-Orbit Coupling
VASP Vienna Ab-initio Simulation Package
xx
INTRODUCTION
In the field of multiferroics and magnetoelectrics, there are still some open ques-
tions and challenges to be addressed that concerns the achievement of a room
temperature multiferroic or a magnetic-field switchable polarization. The break-
through results we present in this thesis stand for the realization of a strain engi-
neering multiferroism/magnetoelectric in a class of materials that nobody has an-
ticipated. Moreover, in the community of ferroelectrics no one would have looked
at ionic compounds with very low polarizability. Besides, due to the low polariz-
ability and large electronegativity of the fluorine atom, novel and unexpected non-
linear strain/polarization properties can be found, which to our knowledge, has
never been reported before in any ferroelectric perovskite-based systems. Thus,
the unexplored behavior in perovskite fluorides is very appealing at the funda-
xxi
xxii introduction
mental level —that it has never been considered— and at the experimental level
—that it has never been observed— and it opens totally new and unexpected
perspectives in the hot field of ferroelectrics and magnetoelectrics, which is in a
constant search for new materials that can be addressed experimentally and drive
novel responses and applications.
From the industrial point of view, fluorides have proven to be of high inter-
est for numerous long term applications, such as fluorides-based glasses with a
large thermal expansion, low refractive and non-linear index [��, ��], strong mag-
nets with an optically transparency in the visible light [��, ��], electrochemical
devices, solid-state batteries, gas sensors, and electrochromic systems [��] or cata-
lyst surfaces based on metal fluorides such as AlF3 [��]. Thus, adding multiferroic
properties in fluoroperovskites would open an exciting opportunity to their use
in novel and extended applications.
At the end of this document several appendixes are presented. The magnetic
properties —exchange constants, group theoretical details, etc— of fluoride per-
ovskites and post-perovskites is explained in Appendixes A and B. We also present
in Appendices C, D, and E additional side studies we performed during our thesis
such as strained BaCoF4 thin films, the role of water adsorptions in STO/LAO
interfaces and the spin texture at the SrBO3 (B = Ti and Hf) surfaces in the
presence of an oxygen vacancy respectively.
Part I
D E N S I T Y F U N C T I O N A L T H E O RY A N D
ELECTRONIC STRUCTURE
Ĥ | i = E | i , (�.�)
where the Hamiltonian operator (Ĥ) acting over the wavefunction ket (| i) has
as a result, the energy eigenvalues (E) times the wavefunction. The Ĥ operator
for the ions composed by the nuclei and electrons is described by the kinetic and
potential energy interactions taking the expression:
Ĥ(R, r) = T̂N (R) + T̂e (r) + V̂NN (R) + V̂ee (r) + V̂eN (r, R), (�.�)
� The time-dependent Schrödinger’s equation present the form: Ĥ(t) | i=ih( @t
@
)| i
�
� theoretical background
here, the operators for the kinetic energy of the nuclei (T̂N ) and the electrons
(T̂e ) are as follows:
h2 X 1 2
T̂N (R) = - r , (�.�)
2 MI I
I
h2 X 2
T̂e (r) = - ri , (�.�)
2me
i
by notation, the coordinates related to the ions and electrons are given in upper
and lower cases subscripts respectively (i.e. I, J, etc. and i, j, etc.). MI and me
are the masses for the ions and the electrons respectively.
The electrostatic potentials for the ion-ion, electron-electron, and electron-ion
interactions are given by:
e2 X ZI ZJ
V̂NN (R, R) = , (�.�)
2 kRI - RJ k
I6=J
e2 X 1
V̂ee (r, r) = , (�.�)
2 kri - rj k
i6=j
X ZI
V̂eN (r, R) = -e2 , (�.�)
kri - RI k
i,I
Ĥe (R, r) = T̂e (r) + V̂NN (R) + V̂ee (r) + V̂eN (r, R). (�.�)
The ion-ion potential, V̂NN is recognized such as the Madelung energy, which
can be easily computed in periodic systems. The exact calculation of this term
can be reviewed in the appendix F by Kaxiras [��]. The V̂eN potential is now
independent of the ion RI coordinates and can be represented only as a function of
the electron positions (ri ). After these considerations, the electronic Hamiltonian,
Ĥe , assuming a fixed nuclei, is described by the sum shown in Eq �.�:
X
Ĥe = T̂e + V̂ee + V̂ext = T̂e + V̂ee + v(ri ). (�.�)
i
�.� density-functional theory �
Z
N = n(r)dr. (�.��)
Z Z
n(r) = N . . . | (r1 , r2 , . . . , rN )|2 dr1 , dr2 , . . . , drN , (�.��)
here, if we think into invert the dependence of terms in the expression in Eq.
�.��, we will see that the wavefunction will become = [n(r)], which is indeed
� theoretical background
a functional of the density� . Even if the later reasoning looks quite simple, it is
indeed the foundation of the density-functional theory as shown next.
Theorem �: The electron density that minimizes the energy of the overall
functional is the true electron density corresponding to the full solution of the
Schrödinger’s equation [��].
The proofs for these theorems are beyond the aim of this chapter and it can be
find in several electronic structure textbooks such as Martin [��]. The previous
theorems establish that there is a one-to-one relationship between the ground
state energy and the electronic density n(r) determined by a unique external
potential. Besides, if the “true” functional is known, then, we can perform a min-
imization of the energy in order to find the related electron density.
Thus, just after considering these theorems, we can redefine the energy (pre-
sented in Eq. �.��) as:
being n(r) the electronic density for the ground state and F[n(r)] a functional of
the density that regroups the kinetic and potential energies added to the external
potential energy. Nevertheless, in all the development shown by Hohenberg and
Kohn, the F[n(r)] is unknown and to find a solution to the Schrödinger’s equation
it is still a problem at this stage of the framework. Then, it was only after the
seminal work by Kohn and Sham [��] that some light was shed onto this problem.
� As a simple example of a functional lets take first the function f(x) = x2 + 1, then a func-
R1
tional F[f(x)] of the function f(x) can be written as: F[f] = f(x)dx. After evaluating the later
-1
expression we have that F[f] = 38 .
�.� density-functional theory �
The main idea of Kohn and Sham was to map the many-electron system into
a single-particle description [��]. For that purpose, the F[n(r)] functional was
redefined as:
here, Tks is the kinetic energy of a system with density n but, excluding the
electron-electron interactions. VH is known as the Hartree potential derived from
the classic Coulomb potential for the electrons. Finally, Exc defines the exchange-
correlation energy missing in Tks and VH . It is important to note at this point
that Tks is not the real kinetic energy, nonetheless, its magnitude is comparable
and can be computed without any approximation. Then, all the terms in the Eq.
�.�� can be evaluated exactly, except for the exchange-correlation energy, which
as we will see, it is the main concern in the theory. As a next step, by considering
the Eqs. �.�, �.�, �.��, and �.�� we have the Schrödinger’s equation with the
form:
" Z #
h2 2 n(r 0 ) 0
- r + dr + vxc (r) + v(r) | ii = ✏i | ii , (�.��)
2me |r - r 0 |
Exc [n(r)]
vxc = . (�.��)
n(r)
As we can see in the Eq. �.��, if we regroup the terms with a potential character,
we have the non-interacting one-electron Schrödinger’s equation expressed in the
Eq. �.��:
h2 2
- r + veff | ii = ✏i | ii . (�.��)
2me
The Eq. �.�� will be renamed from now on as the Kohn-Sham equation. Based
in the new Konh-Sham orbitals, | i i, the electronic density can be computed as:
X
nKS (r) = | 2
i (r)| . (�.��)
i
The most basic form, but at the same time, the most used methodology to ap-
proximate the exchange–correlation energy is by means of the local-density ap-
�.� exchange–correlation energy and functionals �
where ✏hom
xc is the exchange–correlation energy per particle in a homogeneous
�
electron gas . Here, it is assumed that the energy per particle at the point r is
only a function of this point, besides, the exchange energy of a homogeneous
electron is known having the form as follows:
3 1
✏hom
x [n] = - (3⇡2 n) 3 . (�.��)
4⇡
The missing correlation part can be computed using the quantum Monte Carlo
approach as shown in Ref. [��]. Surprisingly, even when this approximation looks
quite rudimentary, nowadays, it is one of the most used and gives a good agree-
ment with experimental results (e.g. ferroelectric properties in materials [��]).
Nonetheless, one of its disadvantages is the over estimation of the cohesive en-
ergy that it can be appreciated in smaller lattice cell parameters in comparison
to experimental data. The wrong computation of the electronic band gap with an
underestimation of its value is also a known disadvantage of the approximation.
The updated representation of the exchange–correlation energy to overcome
some misdescribed issues is based onto the generalized gradient approximation
(GGA). In this approximation, the variation of the electron density is taken into
account in the ✏xc term as:
Z
EGGA
xc [n(r)] = n(r)✏xc [n(r), |rn(r)|]dr. (�.��)
�.�.� DFT+(U,J)
here, the correction term EHub is known as Hubbard Hamiltonian for corre-
lated states. Edc is the correction term for the double-counting of the states. This
problem arises because the Kohn-Sham energies already include the interactions
effects within the Hartree and exchange correlations terms. Then, adding the
Hubbard Hamiltonian leads to a double-counting of some interactions that needs
to be corrected. Note that in the EHub and Edc the atomic site index, the state,
and the spin are taken into account with the l, m, and lower and upper indices
respectively. Finally, the ⇢lmm 0 are the occupation numbers of the localized or-
bitals. In order to solve the Eq. �.��, the EHub and Edc terms can be computed
as follows:
X Ul X 0 Ul l l
EHub
[⇢lmm 0 ] - Edc [⇢l ]= ⇢lm ⇢lm 0 - ⇢ (⇢ - 1) , (�.��)
2 2
l j
the KS states.
One of the problems within this approach is that there is no exact determination
of the U and J values, and then, they must be adjusted based on experiment
data such as magnetic moment, magnetic anisotropy, and electronic band gap.
An alternative method to overcome this issue was proposed by Cococcioni and
de Gironcoli [��]. In this method, a linear response approach is used in order
to calculate the interaction parameters entering to the LDA+U functional. This
�.� exchange–correlation energy and functionals ��
The hybrid functionals own their name to the mixture character of the exchange–
correlation energy. This energy is defined in terms of an amount of an orbital-
dependent Hartree-Fock (HF) energy plus and explicit DFT part. As far as we
know, those are among the most accurate methods up to day —in addition to
Dynamical Mean-Field Theory, DMFT [��] and GW [��] methods for example—.
The exchange-correlation functionals can be expressed in their most general form
as:
EB3LYP
xc = ELDA
xc + a0 (EHF DFT
x - Ex ) + ax EB88
x + ac EPW91
c , (�.��)
here the a0 , ax , and ac coefficients are fixed to �.��, �.�� and �.�� values
respectively adjusted to fit atomic and molecular data. EB88 x is the exchange
energy term taken from the Becke’s ���� gradient correlation within LSDA func-
tional [��]. Finally, the correlation term EPW91
c is taken from the Perdew-Wang
�� functional [��].
Another example is the functional developed by J. Perdew, M. Ernzerhof, and
K. Burke (PBE�) [��]:
1
Exc = ELDA
xc + (EHF - EDFT
xc ), (�.��)
4 x
in this functional, �/� of mixing is employed based in several test in which the
EDFT
xc is set to the PBE functional.
�� theoretical background
The B�WC [��] hybrid functional was developed especially to be used in per-
ovskite oxides. In this functional the GGA approach is taken within the Wu-Cohen
[��] representation and it is employing an ��% of mixing as shown in the Eq. �.��:
EB1 GGA
xc = Exc + 0.16(EHF GGA
x - Exc ). (�.��)
EPBE0
x = a[EHF,SR
x (!) + EHF,LR
x (!)[+(1 - a)[EPBE,SR
x (!) + EPBE,LR
x (!)], (�.��)
E!PBEh
xc = aEHF,SR
x (!) + (1 - a)EPBE,SR
x (!) + EPBE,LR
x (!) + EPBE
c . (�.��)
In the Eq. �.�� the !PBEh term converges to PBE0 and PBE for ! ! � and
! ! 1 respectively. For more detailed information about the derivation of the
SR and LR terms, please refer to the reference of Heyd et al. [��].
It is important to note that some of these hybrid functionals will be used in
the analysis of oxide/fluoride interfaces in the chapter �.
strain ⌘ depending the case. Thus, taking the expansion up to the first-order we
can access to the forces (F ⌘ @E/@⌧), polarization (P) and strain constants ( ) as
shown in the Table �.�. When the expansion is carried out up to second-order,
information related to phonons (i.e Interatomic force constants, C), Born effective
charges (Z⇤ ⌘ @2 E/@⌧@E ⌘ @P/@⌧) and the internal strain coupling parameters
( ) can be obtained (see the Table �.�). Other quantities like the clamped-ion
elastic constant (c0 ), the clamped-ion piezoelectric tensor (e0 ), and the optical
dielectric tensor (✏1 ) can be also computed as presented in the Table �.�. The
third order derivatives give information about nonlinear response such as Raman
tensors, optical susceptibilities and electro-optic properties.
X @E 1 X @2 E
E( ) = E(0) + i + i j+
@ i 0 2 @ i@ j 0
i i,j
(�.��)
1X @3 E
+ i j k +...
6 @ i@ j@ k 0
i,j,k
The calculations of the derivatives shown in Eq. �.�� are based on a variational
principle involving a minimization procedure giving also access to non-linear re-
sponses based in the "�n+�" theorem [��]. In this theorem is established that the
energy of the �n+� state can be obtained with the knowledge of the states from
0 through n.
Table �.�: Available physical quantities through the first- and second-order energy
derivatives with = ⌧, ⌘, and E, (see the Eq. �.��). Table adapted from
Ref. [��].
�st -order �nd -order
@/@⌧ @/@⌘ @/@E
@/@⌧ F C Z⇤
@/@⌘ c0 e0
@/@E P Z⇤ e0 ✏1
1
n,k (r)= p un (k, r)eikr (�.��)
⌦
where k is a lattice wave-vector into the reciprocal space, ⌦ represents the cell
volume and n stands for the band number. un is the periodic function that holds
all the information related to the crystal periodicity. Later, for a more practical
view of the Bloch wavefunctions, those are decomposed into a Fourier expansion
of the form:
1 X
n,k (r) =p Cn,k (G)ei(k+G)r (�.��)
⌦ G
Here, Cn,k are the Fourier coefficients and G is a vector in the crystal recip-
rocal space. Therefore, for any reciprocal lattice vector we have that n,k (r) =
n,k+G (r). Nonetheless, from the computational point of view this expansion
must be limited up to a finite number of plane-waves. For this purposes, a cut-off
energy Ecut condition is established as:
h2
|k + G|2 6 Ecut (�.��)
2m
The Ecut value is independent for each crystal, symmetry and atomic species
and then, a correct selection of an energy value should be set in each particular
case which is commonly known as the energy cutoff convergence.
One of the disadvantages of this methodology is that a large number of plane-
waves is required to describe the electronic wavefunctions close to the nucleus.
This due to the faster change of the valence electrons wavefunctions in the region
close to the atomic nucleus. Thus, a very large Ecut —and thus, a large number
of plane-waves— is necessary to describe these wavefunctions correctly in com-
parison to the energy needed for the electrons close the valence states. In order
to overcome this issue, another approximation has been developed known as the
pseudopotential approach as it will be described later.
�.� from the theory to the implementation ��
The Eq. �.�� take advantage of the Born-von Karman periodic boundary condi-
tions, which allow us to represent an infinite periodic crystal into a representative
building block of construction known as unit cell. At the same time, thanks to
the periodic boundary conditions, the wavefunctions can be presented in the form
of the Bloch functions as commented previously. In order to calculate the wave-
funtions from the Eq. �.�� —and consequently, the electronic density as shown
previously—, an integration over the reciprocal space involving all the bands and
the Bloch functions over an infinite k-points number is required. In practice, this
is impossible and needs to be truncated into a finite number of k-points in a
similar way as done with the Ecut . Therefore, the convergence of the calculations
require enough k-points in order to reach a sufficient precision within a reasonable
computation time. The later process is usually known as the k-points sampling
and determine the accuracy level in the calculations (i.e. bigger is the number
of k-points, higher is the calculation’s precision but larger is the computational
time).
Several methods to sample the reciprocal space and then, to select the k-points
mesh have been developed. This is the case for example, of the Monkhorst-Pack
[��] mesh in which essentially the k-points are distributed homogeneously in the
Brillouin zone. In this technique, rows and/or columns of k-points are running
parallel with respect to the reciprocal lattice vectors in the specific Brillouin zone.
The previous considerations within the periodic boundary conditions are very
convenient when bulk materials and/or superlattices are represented. Neverthe-
less, the modeling of molecules, isolated surfaces, thin films and/or material’s
defects require the definition of special computational unit cells. For more details
in these special systems, please check the reference of Sholl and Steckel [��].
�.�.� Pseudopotentials
The pseudopotential approach [��] is based on the idea of replacing the strong
Coulomb potential of the nucleus by an effective ionic potential (i.e. shielding
potential) acting on the valence electrons only. This approximation can be done
because most of the electronic properties of materials are mainly determined by
the valence electrons while the core electrons can be considered as frozen. The
later statement has a few exceptions in which, deeper electronic levels in compar-
ison to the valence states have an influence on the properties (e.g. ferroelectric
materials under pressure [��, ��]). Thus, the pseudo wavefunction have no radial
nodes in the core region and the pseudopotential is smoother than the all-electron
potential near to the nucleus but identical beyond the cut-off radius, rc . The key
feature of a pseudopotential lies in its transferable and accurate character, then,
even if the pseudopotential was constructed in one environment (the atom itself),
it should faithfully reproduces the properties of a material when it is placed in
different environments such as molecules, solids and ions.
�� theoretical background
The PAW method was proposed by Blöchl [��] in ����. This method deal with
the main shortcoming of the pseudopotential approach, in which, most of the
information regarding to the wavefuntion is lost in the region close to the nuclei.
Here is important to recall that in the bonding region the wavefunction is smooth,
but close to the nuclei it displays rapid oscillations that are very demanding on
the numerical representation (i.e. large number of plane-waves in the pseudopo-
tential approach). Then, to address this problem the PAW method proposed the
following mathematical formalism:
X
T̂ = 1 + T̂ a , (�.��)
a
↵
| a
ni = (1 + Tˆa ) ˜ a ˆa ˜ a
i ,T | ii=|
a ˜a
i i-| i i. (�.��)
Here is important↵ to remark that the T̂ operator is only acting inside the sphere,
then, | a
n i = ˜ a when r > ra . Under the assumption that the smooth partial
n c
waves form a complete set inside the augmentation sphere, we can expand the
smooth all-electron wavefunctions as:
↵ X a ↵
˜n = Pni ˜ a
i , (�.��)
i
�.� from the theory to the implementation ��
where Pni
a are the expansion coefficients. If we multiply by T̂ in both sides of the
Eq. �.��Pwe find that the same is also true for the non-smooth wavefunctions (i.e.
| in = i Pni a | a i). Therefore, the P a expansion coefficients are the same for
n ni
both, smooth and non-smooth wavefunctions. ↵In order to find such coefficients,
we choose a smooth projector functions, p̃a i that satisfies the orthogonality
condition and must be complete inside the augmentation spheres, then:
hp̃a ˜a
i| ji = ij . (�.��)
X
|˜a a
i i hp̃i | = 1. (�.��)
i
⌦
Then, multiplying the Eq. �.�� by p̃a
i and using the Eq. �.�� condition we
have that:
X X
hp̃a ˜
i | ni =
a
Pnj hp̃a ˜
i | ji =
a
Pnj ij
a ,
=Pni (�.��)
j j
X X
T̂ a = T̂ a · 1 = T̂ a |˜a a
j i hp̃i | = (T̂ a | ˜ a a
j i) hp̃i |
i i
X (�.��)
= (| a ˜ a a
i i - | i i) hp̃i | ,
i
Then, the full operator (based in the Eq. �.��) can be expressed as:
X
T̂ = 1 + (| a ˜a a
i i - | i i) hp̃i | , (�.��)
i,a
↵ X X ↵ a
| ni = ˜n + | a a
i i Pni - ˜a
i Pni . (�.��)
i,a i,a
The relation �.�� can be schematically represented as in the Fig. �.�. In the
figure it can be noted that the full wavefunction is expressed as the sum of the
pseudo plane waves (first term in Eq. �.�) plus the all-electron “on-site” interaction
(close to the nuclei and inside de augmentation spheres, second term in Eq. �.�)
minus the “on-site” part from the pseudo wavefunctions (third term in Eq. �.�).
�� theoretical background
= + -
Figure �.�: Schematic representation of the main PAW method that exemplifies the Eq
�.��.
In this way, we are keeping the information about the wavefunction around the
core region without expending large computational time. The later because the
number of plane waves is considerably reduced in comparison to the regular pseu-
dopotential approach discussed in the previous section. With this, the accuracy
is significantly increased and the computational efforts are reduced.
where n 0 (r) is the spin-density matrix, which takes into account a non-
collinear framework in the �⇥� matrix allowing to take into account the mag-
netization spin-density variations (i.e. its magnitude and direction). The EHK
term is the zero-field energy functional as shown in the previous sections. H is
the applied magnetic field, µ0 is the magnetic permeability of free space, and
mtot is the total spin magnetic moment of the system obtained from the spatial
integration of the m(r). In the same way as shown previously, the Eq. �.�� can be
solved by the Kohn-Sham approach into the single-particle spinor orbitals schema.
Then, in addition to the usual terms, an external magnetic field acting on the
spinor orbitals is obtained and the related potential shift have the form:
!
g Hz Hx + iHy
V 0 = - µB µ0 , (�.��)
2 Hx - iHy -Hz
�.� from the theory to the implementation ��
here, the g term is the Landé factor, Hj (with j= x, y, and z) are the magnetic
field vector components and µB is the Bohr magneton. This implementation allow
us the computation of several properties of interest such as the magnetoelectric
coupling as shown in the Chapter � for the NaMnF3 compound. This, opening the
possibility to unveil the properties in multifunctional compounds with a related
magnetic behavior.
VASP— The Vienna ab-initio simulation package (VASP) is the main code used
in the development of this thesis [��, ��]. Projector Augmented-Waves (PAW)
[��] pseudo-potentials were used to represent the valence and core electrons. The
exchange correlation was represented within the GGA - PBEsol parametrization
[��]. The magnetic character of these systems was include and d -electron localiza-
tion corrected by means of the GGA+U within the Liechtenstein formalism [��],
for which, the U value was found by fitting the gap and the magnetic moment as
a function of this parameter. The periodic solutions of these crystalline structures
were represented by using Bloch states with fully converged k-points and energy
cut-off values. The spin-orbit coupling was include into non-collinear calculations
[��]. Born effective charges and phonon calculations were performed within the
density functional perturbation theory (DFPT) [��]. The dynamical matrix was
then unfolded using the Phonopy code [��]. The ferroelectric spontaneous polar-
ization computed through the Berry phase approach [��].
CRYSTAL��— For big cells and large number of atoms, hybrid calculations
have been performed using the linear combination of atomic orbitals (LCAO) as
implemented in the CRYSTAL code [��]. This code is based in a Gaussian-type
basis sets including polarization orbitals with the advantage of decrease consid-
erably the computation time comparing to plane-wave basis set as in VASP.
S T R U C T U R A L P R O P E RT I E S O F
PEROVSKITES ABF3
��
�� vibrational behavior in fluoroperovskites abf 3
lar ferroelectric ground state [��]. This was first attributed by a subtle balance
between the short-range (SR) forces and the long-range (LR) Coulomb interac-
tions in which the later destabilizes the crystal toward a ferroelectric state [��].
Later on, first-principles calculations evidenced that the microscopic key ingredi-
ent for ferroelectricity is a dynamical charge transfer between the oxygen anions
and the B-cation due to hybridization between O-�p and empty B-site d-orbitals
[��]. This effect is highlighted by the presence of anomalously large Born effective
charges (Z⇤ ) on these atoms [��] and by an unstable polar phonon mode at the
center of the Brillouin zone of the cubic high temperature phase [��]. Thus, to
our knowledge, almost all experimental and theoretical researches in the field of
magnetoelectrics and multiferroics have been dedicated to oxide systems, with
special focus to perovskite structures and only a few works and reports —in com-
parison to oxides— have been focused in analogous systems such as fluorides.
Then, the understanding of physical phenomena is particular to each material
and cannot be generalized. Here, as an alternative, emerge some structures based
on fluoride systems with a similar crystal structure than the already discussed
oxides, such as KBF3 [��], and NaBF3 [��, ��]. In these systems, the physical
phenomena such as superexchange interaction, stabilization of octahedral-tilting,
and possible multiferroic or magnetoelectric behavior is still unclear. According
to our previous discussion, we carried out a research on ABF3 fluoroperovskites
in order to explore the possibility to include this material as a new multiferroic
in this selected group.
In terms of the tolerance factor (t), perovskite fluorides follows pretty well the
trend of the phase as a function of the t factor. For example, for values below t =
�.�� the ABF3 adopts the R�c structure, the Pnma symmetry can be found for t
between �.�� to �.��. The Pm3̄m in �.�� to �.�� range and P�3 /mcm from �.��
to �.�� respectively. There are a few exceptions due to the Jahn-Teller distortion
in KCuF3 , KCrF3 with a I�/mcm structure [��]. Besides, taking into account the
large change on the ionic radii of the A-site going from Li to Cs, it is expected
to observe a large change in the expected structure and then, the vibrational
landscape that suppress the ferroelectricity. A systematic study of the structural
instabilities in the high symmetry cubic halide perovskites was carried out in the
��’s using interionic potentials [��]. Phonon frequencies were calculated at the ,
X, M and R high symmetry points of the cubic Brillouin zone for ABX3 , with A
= Li, Na, K, Rb or Cs, B = Be, Mg or Ca, and X a halide (F, Cl, Br or I). While
none of these perovskites has a ferroelectric ground state —although some, mostly
Li members of the series have the non-perovskite ferroelectric LiNbO3 structure
as their ground state—, �� of the �� cases were found to show a ferroelectric
(FE) instability. The instability was usually absent in compounds with a larger
A cations, which anyway usually form in a non-perovskite hexagonal structure,
and so it was concluded that the Na compounds should be the most promising
candidates for ferroelectricity. It was found, however, that in all cases a competing
antiferrodistortive (AFD) zone boundary instability corresponding to rotations
�.� computational details ��
and tilts of the X octahedra around the B cation, dominated over the weaker
polar instability, resulting in a centrosymmetric AFD ground state [��]. Such
competition between AFD and FE instabilities has been widely discussed in the
perovskite oxides literature [��, ��, ��, ��] and recently shown, in the case of
Pnma perovskites to result from a mutual coupling to an X+ 5 mode involving
antipolar displacements of the A-cations and X-sites along the [���] axis. Indeed
artificially removing the X+5 mode in the calculations for perovskite oxides was
shown to induce a ferroelectric ground state [��]. Thus, in the next sections we
are going to go further in this vibrational analysis in order to obtain a deeper
understanding of the structural properties in fluoroperovskites to be later used
in the “tuning” and study of their properties.
period IV. This with the aim to understand the possible influence of ionic radii
size and magnetic ordering in the structure and the vibrational behavior. The
later will help us to understand the ground state at the fluoroperovskite familie
by studying and knowing the nature of the displacive transitions posible in the
perovskite fluorides. In this type of transition, the modes that drives the trans-
formation are the same phonon modes driving the lattice dynamics. Therefore,
it is of high importance to understand the vibrational landscape of the materials
investigated on this thesis. For a better understanding of the nature of the dis-
placive transitions and the theory beneath it, the reader can check the references
from Dove [��] and Lines and Glass [��].
First at all, we start by defining the phonon frequencies and eigenvectors within
the harmonic approximation as:
X
0 0
D̃k↵,k 0 (q) mq (k ) = !2mq mq (k ↵) (�.�)
k 0,
where ↵ and are the cartesian coordinates of displacement, and ’ are the
atoms of interest, !mq is the phonon frequency of a m mode at a wave vector q,
mq is the associated
p phonon eigenvector corresponding to the eigendisplacement
defined as: ⌘mq = M mq and finally, D̃k↵,k 0 is the dynamical matrix, formed
by the Fourier transform of the interatomic force constants matrix (IFC, C̃k↵,k 0 )
as:
X p
D̃k↵,k 0 (q) = C̃k↵,k 0 (q)/ Mk Mk 0 (�.�)
k 0,
In the real space the IFC matrix is defined in terms of the energy of the system
as follows:
@2 E
C↵, 0 (l, l 0 ) = 0)
, (�.�)
@⌧↵ (l)@⌧ 0 (l
�.� ionic A - and B -site radi i size effects ��
where E is the total energy, l and l’ are the cells of respective atoms and ⌧ is
the displacement. In Fig. �.� are presented the fully relaxed phonon-dispersion
bands for NaMnF3 computed in the high temperature symmetric cubic Pm3̄m
perovskite [��]. Structural instabilities (i.e. imaginary modes represented as neg-
ative frequency values) can be seen from the phonon-dispersion curve through
the entire Brillouin zone. A polar ferroelectric mode is observed at zone-center
-point that has been also found in all ferroelectric systems based on oxides such
as BaTiO3 , PbTiO3 , and antiferroelectric PbZrO3 [��].
Additionally, two instabilities with high intensity in frequency value were also
observed at R-point showing a flat branch along the R–M path characteristic of
the AFD modes reported in some oxides too. It is believed that this geometrical
instability competes with FE mode and, in some cases can completely suppress
�� vibrational behavior in fluoroperovskites abf 3
it in the ground state. Thus, all these instabilities combined leads most of the
NaBF3 compounds to an orthorhombic Pnma (��) ground state, in which, fer-
roelectricity is totally suppressed by symmetry. As an exception to the rule, the
NaCaF3 case has been reported in a rhombohedral R�c structure in ground state
symmetry with possible ferroelectric ordering [��]. Nevertheless, up to now there
is no experimental results corroborating the FE behavior in this material. The
phonon-dispersion bands for B = Ca, V, Ni, and Zn are similar to B = Mn and
because of that, are not presented here.
According to these vibrational results, we have found for all structures a clear
dependence of FE-instability as a function of the ionic B-site radii for B2+ state,
of course, taken into account the number of coordination (i.e. B-F6 bonding) [��].
Then, as the ionic B-site radii decrease, the FE-mode decrease as well as a function
of this parameter (see Table �.�). Additionally, we observed the same behavior
for one of two instabilities at R-point, meanwhile the second one remains almost
constant. In the same table we also presented the computed values for BaTiO3
for comparison with oxide systems.
Figure �.�: Polar 4- mode, into the Bethe’s notation [��], as a function of the A- and
B-site in ABF3 perovskites, the imaginary modes are presented as negative
values.
Now, we are going to talk about the A-site effect, for this proposes the B-site
was fixed to Ni case, this in order to keep small lattice parameters. From Table
�.� it can be appreciated that FE polar instability present a more critical depen-
dence with respect to the ionic A-site radii size than B- one. The later findings
are presented graphically in the Fig. �.�. Therefore, an increase of this parameter
can totally suppress the FE imaginary mode and AFD instabilities, which leads
the structure to a totally different crystal ground state symmetries as listed in
Table �.�. This behavior can be explained in terms of eigendisplacements, ⌘ (see
Table �.�) for this mode, in which the dominant contribution in the FE-mode is
precisely the A-site displacements. Increasing the A-site radii, the displacements
�.� ionic A - and B -site radi i size effects ��
decreases as well and with that a strong reduction of the frequency it is appreci-
ated, meanwhile, the imaginary mode is suppressed for A = K, Rb, and Cs (see
Table �.�).
��
Table �.�: Summary of most unstable high symmetry point phonon eigenvalues, relaxed lattice constants for cubic Pm 3̄ m symmetry,
ionic A and B sites radii and reported ground state symmetry for all compounds [��]. HS: High Spin, LS: Low Spin. All modes
+ + + -1 .
in ( 4- , R+
4 , R 5 , X 5 and M 3 ) in cm
-
Compound ao [Å] B-radiia [pm] Godschmidt, t FE - 4 R+
4 R+
5 X+
5 M+
3 Ground State [��]
BaTiO3 �.���� ��.� �.�� ���i Non Non ���i ���i P�mm [���]
NaCaF3 �.���� ��� �.�� ���i ���i ���i ���i ���i R�c [��]
NaMnF3 �.���� �� (HS) �.�� ���i ���i ��i ��i ���i Pnma
NaVF3 �.���� �� �.�� ��i ���i ��i ��i ���i Pnma
NaZnF3 �.���� �� �.�� ��i ���i ��i ��i ���i Pnma
NaNiF3 �.���� �� (LS) �.�� ��i ���i �� ��i ���i Pnma
A-radii [pm]
CsNiF3 �.���� ��� �.�� ��� ��� ��� ��� ��� P�3 /mmc, z =�
RbNiF3 �.���� ��� �.�� ��� ��� ��� ��� ��� P�3 /mmc, z =�
KNiF3 �.���� ��� �.�� ��� �� ��� ��� �� Pm3̄m
NaNiF3 �.���� ��� �.�� ��i ���i ��i ��i ���i Pnma
LiNiF3 �.���� �� �.�� ���i ���i ���i ���i ���i R�c b
a Ions are �-coordinate [��]
vibrational behavior in fluoroperovskites abf 3
structure perovskites as a function of the ionic radii ratio. We start from the
respectively tolerance factor (i.e. t, Table �.�) defined by Goldschmidt as follows:
rA + rX
t= p , (�.�)
2(rB + rX )
where rA , rB , and rX are the ionic radii for the A-, B- and X-sites respectively.
It can be noted that the Eq. �.� shows a close relationship between the tolerance
factor and ionic sites radii, which at the same time, it can be related to the ground
state and the FE/AFD instabilities ratio. Then, for smaller t values (i.e. LiNiF3
and NaCaF3 ) we observed that the FE frequency is higher in comparison to
AFD one. Here, this partial conclusion suggest that the FE freezing is sufficient
enough to lead these systems from the cubic high temperature structure to a
rhombohedral R�c ground state as previously reported for NaCaF3 but, to out
knowledge, found for a first time in LiNiF3 . This opens the possibility for a novel
G-type AFM multiferroic fluoride-based perovskite. For other compounds, an
orthorhombic Pnma structure is expected at tolerance factors that lies between
�.�� to �.��. From �.�� to �.�� a cubic Pm3̄m ground state is expected. At higher
values, above t = �.�� are anticipated different hexagonal symmetries. A few
exceptions such as KCuF3 and KCrF3 have been reported due to Jahn-Teller
distortions [��]. All these results are in full agreement with the our computed
vibrational properties listed in Table �.�.
In Fig. �.� are presented the results for the ferroelectric instability as a function
of the tolerance factor. Here, the LiBF3 and NaBF3 series are shown, in which, it
can be observed that the values of the frequency for the Li family are well above
��� icm-1 . These values are close to the ones observed in oxides perovskites. For
the Na family we found frequency values between �� to ��� icm-1 showing and
stronger dependence on the B-site radii.
Figure �.�: Ferroelectric instability in LiBF3 and NaBF3 (B = Ca, Mn, V, Zn, and Ni)
as a function of the tolerance factor, t.
�.� ionic A - and B -site radi i size effects ��
@P
Z⇤,↵ = ⌦ | , (�.�)
@⌧↵, E=0
where ↵ and are the cartesian coordinates, P is the polarization and E is the
macroscopic electric field. Our computed Z⇤ values for all fluorides in question
and BaTiO3 compounds are presented in Table �.�. For both, A- and B-sites
we found that Z⇤ s are close to nominal charge (i.e. A1+ , B2+ ). However, for
F-ions the Z⇤ are �.� e- away from its valence, which could be possibly due to
a weak covalent interaction in the F-(B-site)-F bonding. It is clear that these
values are lower than in the oxides case, but are strong enough to contribute
to the superexchange magnetic interaction. The later interaction responsible for
the antiferromagnetic ordering observed experimentally in most of perovskite
fluorides —with a magnetically active B-site— at their ground state structure
[��]. In the same table, the Z⇤ components for barium titanate are also presented.
Anomalous Z⇤ for Ti and O sites parallel to the bonding are observed as expected
due to the strong covalent contribution that allows the stabilization of the FE
ground state in this oxide. These Z⇤ close to nominal charges stand for another
substantial difference of the fluorides in comparison to oxides perovskites that
can be exploited.
In order to achieve a better understanding of the vibrational properties in
fluorides, we analyzed the interatomic force constants (IFC) calculated from Eq.
�.�. The analysis of IFCs reveals information about interatomic interactions that
are predominant in the vibrational instabilities and lead to a stable structure.
We present our calculated on-site force constants (on-site IFC) in Table �.�. The
on-site IFC are positive for all the elements in the ABX3 compounds, showing
that they are all stable under an isolated atomic displacement. Then, the phonon
instabilities observed in phonon dispersions and discussed above, only can be
removed or compensated by a cooperative movement of all atoms in order to
gain energy and lower the total enthalpy of the system. One exception inside the
analyzed compounds is the negative on-site IFC for Li atoms in LiNiF3 perovskite
suggesting a necessary collective movement of Li ions in lattice to reach a stable
state. This fact can explain the R�c crystalline group expected as ground state for
LiNiF3 in which a huge displacement and distortion of cubic structure is required
to achieve the structural transition. From the Table �.� can be also noted that B
= V and Ni presents the highest values of on-site IFC values for A, B, and F sites,
this fact suggest a low movement of self-site ions from its position in comparison
with the other compounds. The remarkable low value in the on-site IFC for Na in
NaCaF3 suggest a possible A-site distortion needed to stabilize the structure and
then, suppress the unstable phonon branch exhibited at zone-center. Additionally,
�� vibrational behavior in fluoroperovskites abf 3
Table �.�: Born Effective Charges (Z⇤ , e- units) in for all compounds computed at
relaxed cubic structure computed according to the Eq. �.�. The X? and Xk
are the components perpendicular and parallel to the B–X bonding.
Z⇤ (ABX3 ) [e- ]
Compound A-site B-site X? Xk
BaTiO3 �.��� �.��� -�.��� -�.���
NaCaF3 �.��� �.��� -�.��� -�.���
NaMnF3 �.��� �.��� -�.��� -�.���
NaVF3 �.��� �.��� -�.��� -�.���
NaZnF3 �.��� �.��� -�.��� -�.���
NaNiF3 �.��� �.��� -�.��� -�.���
CsNiF3 �.��� �.��� -�.��� -�.���
RbNiF3 �.��� �.��� -�.��� -�.���
KNiF3 �.��� �.��� -�.��� -�.���
NaNiF3 �.��� �.��� -�.��� -�.���
LiNiF3 �.��� �.��� -�.��� -�.���
Table �.�: On-site interatomic constants (on-site IFC, in eV/Å2 units) for the ABX3
flouroperovskites. The on-site IFC for the BaTiO3 and PbTiO3 are also pre-
sented for comparison. The Fx=y sites lies in A-F plane meanwhile the Fz lies
along B–F bonding.
On-site IFC for ABX3 [eV/Å2 ]
Compound A-site B-site X X
x=y=z x=y=z x=y z
BaTiO3 �.�� ��.�� +�.�� +��.��
PbTiO3 [��] �.�� ��.�� +�.�� +��.��
NaCaF3 �.�� ��.�� �.�� ��.��
NaMnF3 �.�� ��.�� �.�� ��.��
NaVF3 �.�� ��.�� �.�� ��.��
NaZnF3 �.�� ��.�� �.�� ��.��
NaNiF3 �.�� ��.�� �.�� ��.��
CsNiF3 �.�� �.�� �.�� �.��
RbNiF3 �.�� ��.�� �.�� �.��
KNiF3 �.�� ��.�� �.�� ��.��
NaNiF3 �.�� ��.�� �.�� ��.��
LiNiF3 -�.�� ��.�� �.�� ��.��
Figure �.�: Selected IFC interactions in the � ⇥ � ⇥ � supercell. A - B - and F-sites are
depicted in green, violet and grey respectively. Additionally, the atoms are
labeled to correlate the IFCs shown in the Table �.�.
��
Table �.�: Selected IFC according to interactions illustrated in the Fig.�.� for the A -, B -, and F-sites respectively.
IFC for ABF3 compounds one quarter of unit cell [eV/Å2 ]
Compound A1 -A3 A1 -A2 A1 -B1 B1 -B2 A1 -F1 B1 -F1
xx=yy zz k ? xy xx k ? xx=yy xy zz k ?
NaCaF3 -�.�� +�.�� -�.�� +�.�� -�.�� -�.�� -�.�� +�.�� +�.�� +�.�� -�.�� -�.�� -�.��
NaMnF3 -�.�� +�.�� -�.�� +�.�� -�.�� -�.�� -�.�� -�.�� +�.�� +�.�� -�.�� -�.�� -�.��
NaVF3 -�.�� +�.�� -�.�� +�.�� -�.�� -�.�� -�.�� +�.�� +�.�� +�.�� -�.�� -�.�� -�.��
NaZnF3 -�.�� +�.�� -�.�� +�.�� -�.�� -�.�� -�.�� +�.�� +�.�� +�.�� -�.�� -�.�� -�.��
NaNiF3 -�.�� +�.�� -�.�� +�.�� -�.�� -�.�� -�.�� +�.�� +�.�� +�.�� -�.�� -�.�� -�.��
CsNiF3 -�.�� +�.�� -�.�� +�.�� -�.�� -�.�� -�.�� +�.�� -�.�� -�.�� -�.�� -�.�� -�.��
RbNiF3 -�.�� +�.�� -�.�� +�.�� -�.�� -�.�� -�.�� +�.�� -�.�� -�.�� -�.�� -�.�� -�.��
KNiF3 -�.�� -�.�� -�.�� +�.�� -�.�� -�.�� -�.�� +�.�� -�.�� -�.�� -�.�� -�.�� -�.��
NaNiF3 -�.�� +�.�� -�.�� +�.�� -�.�� -�.�� -�.�� +�.�� +�.�� +�.�� -�.�� -�.�� -�.��
LiNiF3 -�.�� +�.�� -�.�� +�.�� -�.�� -�.�� -�.�� +�.�� +�.�� +�.�� -�.�� -�.�� -�.��
vibrational behavior in fluoroperovskites abf 3
�.� ionic A - and B -site radi i size effects ��
Figure �.�: Double-well energy obtained at the FE-mode freezing in NaBF3 compounds.
The gain in energy is in agreement with previous results that shows that the
bigger is the B-site, bigger is the FE instability and with this, the possibility
to be ferroelectric. ± denotes the full amplitude of the freezing.
Additionally to on-site IFC, some selected IFC for interactions were computed
within a quarter of a �⇥�⇥� supercell as labeled in Fig.�.� and presented in
Table �.�. The analysis were taken into account for components parallel and per-
pendicular to the direction of line that connects the related atoms. This analysis
was performed with the aim to study some selected interactions and analyze the
relationship between the ions that present the highest contribution in the total
energy. First, it is noted the significantly high repulsion force between (B1 –F1 )k
sites, which induce the movement of all B-sites with the fluorine ions instead of
against them, like in barium and lead titanates. The IFCs values as a function of
A-site radii shows that B–F repulsion interaction decrease when A-radii increase
as well, and then, a freezing of octahedra tilting and stabilization of cubic struc-
ture is achieved in the KNiF3 case. Besides, the B1 –B2 k IFCs values and the
high on-site IFC obtained for the center-sites suggest that B-ions cannot move
easily in comparison with the center-site displacement observed in some oxides
(e.g. barium titanate). Nevertheless, these B-sites are still presenting a small con-
tribution in the FE-mode that cannot be neglected as discussed above based on
eigendisplacements values. The contributions from another interactions such as
A–A and A–F listed in the Table �.� are smaller than those discussed previously
and do not present a considerable contribution.
Moving forward, we analyze the energy landscape of the ferroelectric mode. In
Fig.�.� are presented the double-well energy curves obtained through the freezing
of ferroelectric instabilities along z -axis. From this figure it is observed a clear
dependence of energy minimum as a function B-site, in which the higher value
of FE frequency corresponds to a lowest energy value. These results are in full
agreement with our previous analysis and are summarized in the next section.
In NaCaF3 case, the double-well analysis suggests a prominent predisposition of
�� vibrational behavior in fluoroperovskites abf 3
these materials to freeze the FE ground state over the AFD one due to the high
energy gain.
BaTiO3
Geometric%FE%ins
(a) BaTiO3 because%no%ch
Geometric%FE%inst
(a) BaTiO3 BaTiO3
(a) BaTiO
(a) 3
transfer%and%smal
because%no%cha
in%a%big%cell%(4
(a) BaTiO3
transfer%and%small
in%a%big%cell%(4.1
NaMnF3
(b) NaMnF3
(b) NaMnF3
Figure 3.4: Double-Well energy curves(b)
Figure 3.4:Andrés
Double-Well energy
Camilo García curves
Castro,
and
NaMnF
and
(b)
NaMnF
schematic representation of atomic
3 3
schematic representation of atomic
NaMnF
[email protected]
displace-
displace-
ments in different models of ferroelectricity for oxides 3 a) BaTiO3 and b) NaMnF3 WVU-Physics
ments3.4:
Figure in different models
Double-Well of ferroelectricity
energy for oxides
curves and schematic a) BaTiO3 and
representation b) NaMnF
of atomic displace-
3
Figure 3.4:Andrés
Double-Well energy
Camilo García curves
Castro, and (b)
schematic representation of atomic displace-
[email protected]
ments in different models of ferroelectricity for oxides a) BaTiO3 and b) NaMnF3 WVU-Physics D
mentsFerroelectric
3.3 in different models of ferroelectricity
Instability for oxides a) BaTiO
Enhancement 3 and b)3NaMnF
in NaMF Fluorides
3
3.3 �.�:Ferroelectric
Figure Double-Well energy Instability Enhancement
curves and in NaMF3ofFluorides
schematic representation atomic displace-
An
3.3 interesting behavior
Ferroelectric of ferroelectricity
Instability and then,
Enhancement of the FE-instability
in NaMF as a function
Fluorides
mentsbehavior
An interesting in different models of
of ferroelectricity ferroelectricity
and then, of thein forNaMF
oxides 3a)Fluorides
3
FE-instability asBaTiO 3 and b)
a function
of 3.3 Ferroelectric
isotropic external strainInstability
has been found Enhancement
in some oxide based compounds like BaTiO
Anof interesting
isotropic external
NaMnF 3 . strain
behavior has been foundand
of ferroelectricity in some
then, oxide
of thebased compounds
FE-instability aslike BaTiO33
a function
and
AnPbTiO 3 [23,behavior
interesting 26], in which the FE ordering can be totally suppressedasbya function pressure.
ofand PbTiO
isotropic 3 [23, 26],
external strainin of ferroelectricity
which
has been the found
FE ordering and then,
in somein can of
oxidebethe FE-instability
totally
based suppressedlike
compounds by pressure.
BaTiO3is
Then,
Then,
this
of isotropic
this
physical
physical
phenomena
externalphenomena
strain has been widely
widely founddiscussedin someinoxide
discussed
oxides
oxides
but
based
but
less
compounds
less
known
known
in fluorides
like
in BaTiO3is
fluorides
As
and
and allPbTiO
PbTiO
proposed the
for previous
[23,
3 study
[23,
26],
in
26],
in
the discussion
in
which
framework
which
the
the and
FE
of
FE this analysis
ordering
research.
ordering
can
can suggest,
be totally
First,
be we
totally the
have mechanism
suppressed
studied
suppressed
by
the
by responsible
pressure.
effects
pressure. of
proposed
Then, this for studyphenomena
3
physical in the framework widely ofdiscussed
this research. in First,but
oxides weless
haveknown studied in the effectsisof
fluorides
for the
pressure ferroelectric
Then, this
pressure ininvibrational
physical
vibrational
mode
properties
phenomena is
properties
based
widely on
ofoffluorides-likethe
discussed in
fluorides-like
A -site movement,
perovskites.
oxides but less
perovskites. In
In figure
which
known3.5
figure in is
3.5
is strongly
presented
fluorides
is presented is de-
proposed for study in the framework of this research. First, we have studied the effects of
pendent
the square
proposed
the of of
square vacuum
foroffrequency
study inspace
frequency thefor or
FE-mode
framework
for free space
FE-mode instability
of for A
this research.
instability as -site
as displacements
aafunction
First,
function we ofof strain.
have studied
strain. related
From
Fromtheexpansionthen,
effects
expansion of to a
pressure
strain zone
pressure ininvibrational
can be
vibrational properties
observed that
properties ofthe
offluorides-like
square
fluorides-like of perovskites.
frequency
perovskites. for Inall
In figure 3.5
compounds
figure 3.5 is
is presented
converge
presented
steric
the effect.
strain
square zone Thus,
ofoffrequency a
can be observed decrease
for of
that the
FE-mode B -site radii,
square of
instability as reduce
frequency
aafunction the
for
of lattice
all compounds
strain. volume
From convergethen,
expansionand
almost
the to
square the same value.
frequency This
for could
FE-mode be a consequence
instability as of
function structural
of strain. distortions
From induced
expansion
reduce
almost
strain the
zoneto free
the
can space
same
be value.
observed for This
A -site
that could
the movement.
be a consequence
square of Also,ofthis
frequency factcompounds
structural
for all explains the
distortionsconverge strongly
induced
bystrain
bymixed
mixedzone can
phenomena be
phenomena observed
between
between that
ionic the
ionic square
charge
charge of
effectsfrequency
effects and
and for
ionic
ionic all
B-sitecompounds
B-site radii
radii converge
that
that needs
needs
correlated
almost
almostto
further to dependence
the
studies.thesame
same
For value.
value. ofThis
FE-mode
This could
could be beisas a function
aaobserved
consequence
consequence ofstructural
of
of A-site
structural radii, thus,induced
distortions
distortions a complete
induced
byfurther
mixed studies.
phenomena Forcompression
compression
between
strain
strain
ionic chargeis observed
effects
aanon-linear
non-linear
and ionic
increase
increase
B-site
of
of instability
radii instability
that needs
in
in
by mixed
opposite
opposite phenomena
direction
direction as as between
observed in
observed instrain ionic
BaTiO
BaTiO charge
3 [23] ineffects
which
[23] in which and
a ionic
suppresion B-site by
a suppresion by pressure radii
pressure that has needs
been
has been
further
furtherstudies.
studies. For compression
Forwas compression isis3observed aanon-linear
non-linear increase of
of instability in
demostrated.
demostrated. Also,
Also, was observed
observed astrain
a“transition”
“transition” observed ofofeigendisplacements
eigendisplacements increase instability
responsible
responsible of
of in
this
this
opposite
opposite direction
direction asas observed
observed in BaTiO
inNaNiF
BaTiO 3 [23]
3 [23]
in which
in which a suppresion
a suppresion by pressure
by pressure has been
has been
mode
mode (figure
(figure 3.6).
3.6).Inwas
Influorides,
fluorides, NaNiF 3 3case
casepresent
present aatotal
total suppression
suppression of
of FE-instability
FE-instability
demostrated.
demostrated. Also,
Also, was observed
observed aa“transition”
“transition” ofofeigendisplacements
eigendisplacements responsible
responsible of
of this
this
atat
external
external strain above
strainIn 2% as
above 2% asNaNiF observed for
observed for present BaTiO in the
BaTiO33 in the same same figure 3.5. This
This different
figure 3.5.FE-instability different
mode
mode (figure
(figure3.6).
3.6). Influorides,
fluorides, NaNiF3 3case case presentaatotal totalsuppression
suppression of of FE-instability
atat
external
externalstrain
strainabove
above2% 2%asasobserved
observedfor forBaTiO
BaTiO3 in inthe
thesame
same figure
figure 3.5. 3.5. This
This different
different
�.� ferroelectric instability enhancement in NaBF 3 fluorides ��
suppression of the FE mode is observed for ionic radii beyond Potassium. Addi-
tionally, it is noted that F-ions present a strong participation in the FE-instability,
and also exhibit the same trend as a function of ionic B-site radii. Besides that,
the ferroelectric polar mode is less dependent on the covalent interaction present
along the B–F bonding, making this phenomena more stable than in oxides. In
order to corroborate this idea, double-well energy curves were computed for the
NaBF3 compounds (As discussed before and shown in Fig.�.�).
To compute this curves, the polar FE phonon-branches were frozen in order to
obtain the energy barrier for this mode. Additionally, we computed the curves
for the independent freezing of each A and B-sites eigendisplacements in fluo-
rides and BaTiO3 oxide. In Fig.�.�a it can be appreciated the energy curves for
which in the oxide case, the lowering in energy is induced by the B-site movement
as widely studied and expected. Any contribution of the A-site (i.e. Ba) was
observed according to the curves. In contrast, for the fluoride case —here choos-
ing the NaMnF3 as example— the largest contribution comes from the A-site
displacements (see Fig.�.�b), and no contribution of B-site was observed. The
later is in agreement with all the previous results that explains the A-site ge-
ometrically driven FE mode. In the same figures are also schematically shown
—with arrows— the atomic displacements for FE-modes in each case. As dis-
cussed above, we found different contributions from the A- and B- sites in the
fluoride and oxide cases respectively. This substantial difference is making this
phenomena completely characteristic of the fluoroperovskites in comparison to
oxides. In conclusion, all the results demonstrates an A-site geometrically-driven
FE-instability mechanism in fluorides rather than a B-site charge-transfer one as
in most of the oxides.
influence of the A-site for pressure regions located close to the relaxed lattice and
expansion strain. In the region of constrain values, below -�% with respect to the
relaxed volume in cubic Pm3̄m structure, the influence of the A-site decreases
and the mode is dominated by B-site displacement. We found that the eigendis-
placements behave similarly in all of the NaBF3 materials, with the exception
of B = Ni. To reinforce this change in A- and B-sites displacements, the F-sites
also present a ⌘z transformation. For ⌘z parallel to B–F bonding, a “transition” of
eigendisplacements from negative, at expansion strain values, to positive values
are observed beyond strain pressures of -�%. The eigendisplacements for F-site
perpendicular to B–F bonding remains negative but goes from ⌘z = -� to ⌘z =
-�� a.u. indicating an increase on their movement amplitude. The eigendisplace-
ments are schematically represented in the Fig.�.� for two regions of strain. The
transformation on the nature of the eigendisplacements is possibly the cause for
the frequency variation as a function of pressure observed in BaTiO3 (Fig.�.�).
In this oxide, the same change in ⌘ is observed for the ferroelectric mode at the
FE-mode suppression when the hydrostatic pressure is applied.
Strain < 0%
Strain > 0%
Figure �.�: Eigendisplacements for NaCaF3 compound showing the ⌘ transition also
found in B -site: V, Mn, and Zn (left). Schematic representation of eigendis-
placement in the fluorides compounds for strain values above and below of
the relaxed structure (right).
In Table �.� we presented the on-site IFC values for all NaBF3 compounds
computed at a fixed volume that is equivalent to an expansion pressure in all
of them. The later analysis in order to understand the B-site radii effect on the
IFC. It can be noted that the A-site IFC constants are negative and close to zero,
this indicating a possible uniform displacement of these sites required to lower
the total energy. The on-site IFC values for F-sites increase as the B-site radii
decrease, going from Ca to Ni. Therefore, the cubic lattice volume reduction leads
to a lower F free space. The Fz values decrease with B-site radii, which implies a
reduction of the attractive interaction force allowing a larger Fz -site displacement.
The close value observed for on-site IFC site of A-ions could be the reason of the
�� vibrational behavior in fluoroperovskites abf 3
Figure �.�: Eigendisplacements for BaTiO3 showing the change of the mode as cause
of suppression of ferroelectric instability as a function of isotropic pressure.
Table �.�: On-site IFC for all NaMF3 compounds at same lattice parameter equivalent
to expansion.
On-site IFC [eV/Å2 ], a = �.�� Å
Compound A-site B-site F F
x=y=z x=y=z x=y z
NaCaF3 -�.�� +��.�� +�.�� +��.��
NaMnF3 -�.�� +�.�� +�.�� +�.��
NaVF3 -�.�� +��.�� +�.�� +�.��
NaZnF3 -�.�� +�.�� +�.�� +�.��
NaNiF3 -�.�� +�.�� +�.�� +�.��
On-site IFCs as a function of pressure are presented in Fig. �.�� for each site.
It can be seen that for a certain value in each case, between -�% to -�%, each
of the Fx=y constants exhibit a “transition” from positive to negative valued,
which implies a change of interaction from attraction to repulsion respectively for
these atoms. This transformation in the nature of interaction possibly induces
some instabilities responsible for the increment of the FE-mode frequency as a
function of pressure observed in the Fig.�.�. It can be also noted that on-site
IFC values for NaNiF3 remain positive, this could explain the different behavior
of this compound in comparison with the rest of the members of this family
discussed above. Additionally, A-, B- and Fz constants are quite larger for B =
Ni in comparison with the Ca, Mn, V, and Zn.
�.� ferroelectric instability enhancement in NaBF 3 fluorides ��
Figure �.��: On-site IFC for A -, B - and F-sites (x=y, and z ) as a function of strain for
each of sites in Na B F3 with B = Ca, Mn, V, Zn, and Ni.
The analysis and study of vibrational properties allow us to understand the nature
of the ferroelectric behavior in perovskite compounds. Based on that, a deeper
analysis was carried out in order to understand the coexistence of ferroelectric,
antipolar, rotational and tilting distortions in the NaBF3 family of fluorides. First
of all, it can be noted that in fluorides the contribution of the antipolar X+ 5
mode is higher than in oxides (compare Fig.�.�� with the values reported in
Ref. [��]). This mode is related to antipolar distortions dominated by the A-site
eigendisplacements. The R+ 5 mode contribution can be appreciated as well which
is almost missing in oxides [��]. R+
5 mode is also related to the A-site motion. M2
+
mode present a smaller effect and, in the range of scale it cannot be observed.
The R+ 4 mode and M3 are also present as in oxides. This fact clearly modifies
+
the balance picture between the contribution of modes to distortion when going
from Pm3̄m to Pnma.
Here, the picture of the competition of modes is more interesting because the
polar ferroelectric instability —present at -point—, the antipolar mode —at
X-point— and the antiferroelectric mode —at R-point— depend as well to the A-
site displacements. Despite the picture in oxides, where the ferroelectric instability
is leaded by B-site displacement meanwhile the A-site dominated antipolar X+ 5
mode suppresses the ferroelectric one in the Pnma symmetry as demonstrated
by Benedek and Fennie [��].
�� vibrational behavior in fluoroperovskites abf 3
Figure �.��: Contribution of modes in the phase transition from Pm 3̄ m to Pnma in the
Na B F3 fluoroperovskites.
Finally, to conclude the ABF3 analysis, we stress that the ionic A and B site
radii size have a direct influence on FE-instability. From the theoretical analysis
carried out and discussed above, the LiNiF3 compound was chosen due to its
higher FE-instability. This compound presents a relaxed cubic lattice parameter
of �.�� Šand its ferroelectric instability is higher than AFD-modes. This fact
could lead to a stable ferroelectric ordering at the relaxed ground state struc-
ture, in which a R�c crystal symmetry type —like in LiNbO3 — is expected as
suggested by Claeyssens et al. [���]. In Fig.�.�� are presented the analysis of in-
stabilities as a function of strain at the Brillouin zone-center and boundary zone.
For strain values around �% of the relaxed volume, the FE-instability shows the
highest frequency value but below -�% the AFD-instability becomes dominant
in the vibrational landscape for this compound. The eigendisplacements for the
�� vibrational behavior in fluoroperovskites abf 3
Figure �.��: Squared frequency at the vibrational instabilities. We can observe the FE-
mode at and AFD-modes at R -point as a function of hydrostatic strain
for LiNiF3 (left). Eigendisplacements of FE-mode instability in LiNiF3 are
also shown (right).
�.� Conclusions
Vibrational analysis were carried out in fluoroperovskites ABF3 . Different behav-
ior for in NaBF3 in comparison with oxides (e.g. BaTiO3 ) was found for B =
Ca, V, Mn, and Zn. For these compounds we found an increase of ferroelectric
instability as a function of hydrostatic strain, this possible due to the “transition”
of eigendisplacements responsible for polar mode. The increase of ionic A-site
radii present a strong influence in FE-polar instability, generating totally differ-
ent ground states symmetries for all compounds going from Li to K. The com-
pound NaNiF3 present a quite different trend in comparison with fluoride family,
for this case, we observed a suppression of FE instability. However, it was not
found change in mode eigendisplacements. The compound LiNiF3 reveal a appar-
ently predisposition to ferroelectric behavior with a strong FE-mode instability
in comparison with the others instabilities observed along M–R branch in the
phonon-dispersion curve. A model for an A-site geometrically driven ferroelec-
tricity in fluorides is proposed according to results and analysis. Some questions
�.� conclusions ��
remain unclear such as: what is the mechanism for cancellation of FE mode at
the ground state? How can be tuned the polar mode in a fluoride structure to en-
gineer a ferroelectric fluoroperovskite? Some of these questions will be addressed
in the next chapters.
Part III
Perovskite materials with stoichiometry ABX3 are well known in material sci-
ences due to their numerous properties and their high interest for technological
applications. These properties are very diverse such as ferroelectricity [��], ferro-
magnetism [���], piezoelectricity [���], pyroelectricity [���] or superconductivity
[���]. Going further, recent researches focused on the achievement of multifunc-
tionalities in a single material that would open the door to new types of applica-
tions such as computer memories or spintronics [���, ���].
One of the main features of these compounds is that their properties can be con-
veniently tuned by external control parameters such as electric and/or magnetic
fields, temperature, strain or pressure [���, ���]. Hydrostatic pressure appeared to
be an interesting external parameter to be considered to enhance our understand-
ing of the perovskite physics, since they exhibit profuse structural, electronic and
magnetic phase transitions [���] under this parameter. A particularly appealing
transition is the perovskite to post-perovskite (pPv ) structural phase transition
(see Fig. �.�a) that has been observed in MgSiO3 crystal in ���� [���, ���]. The
discovery of this pPv phase transition have important consequences in the physics
of the earth mantle where MgSiO3 is abundant, under very high pressure con-
ditions [���, ���]. The pPv transition has then attracted great interest and it
has been observed in other oxide and fluoride perovskites such as MgGeO3 [���],
CaSnO3 [���], NaMgF3 [���] or NaCoF3 [���, ���] (see Table �.� for more details).
Interestingly, the fluoroperovskite crystals have a lower pressure and tempera-
ture pPv phase transition than the oxides and thus making them attractive for
the study of the pPv phase. Furthermore, it has been reported the possibility
of keeping the pPv phase at atmospheric pressure after releasing the pressure
in fluroroperovskites [���], while all the oxides return to their Pnma distorted
perovskite ground state [���].
Besides the interest for the pPv transition in earth mantle studies, we have
recently highlighted the existence of a ferroelectric instability in the cubic phase
of ABF3 perovskites that originates from a geometric effect rather than a dynam-
ical charge transfer, as is it the case in oxides —for more details see Chapter �
and Ref. [�]—. It was demonstrated that ferroelectric polarization in NaMnF3 is
induced by applying an epitaxial strain, an observation that leads to new possibil-
ities of engineering multiferroic materials within the family of fluoroperovskites.
This makes NaMnF3 a good candidate for multifunctional properties based on its
particular behavior as a function of external constraints. Additionally, in the fam-
ily of �d metal transition pPv fluorides, NaMnF3 remains unexplored along the
��
�� prediction of the p p v transition in n a m n f 3
pPv transition. The study of NaMnF3 under pressure is then of high interest for
scrutinizing the origins of its unique properties and its potential multifunctional
character.
Table �.�: Lattice parameters of the pPv phase (a, b, and c) and related temperature
(T) and pressure (P) transition parameters from the Pnma phase to the pPv
phase of fluoride and oxide perovskites. In the pPv phase, the A and B cations
occupy the 4a and 4c Wyckoff positions and the O/F anions occupy the 4c
and 8d positions. Most of the oxides have higher pressure and temperature
transitions than the fluorides.
Crystal a/b/c [Å] P [GPa] T [K] Ref.
Oxides ABO3
MgSiO3 �.���/�.���/�.��� ��� ��� [���]
MgGeO3 �.���/�.���/�.��� �� ��� [���]
MnGeO3 �.���/�.���/�.��� �� - �� ���� - ���� [���]
CaSnO3 �.���/�.���/�.��� �� ��� [���]
CaIrO3 �.���/�.���/�.��� �-� ���� - ���� [���, ���]
CaRuO3 �.���/�.���/�.��� �� - �� ���� - ���� [���, ���, ���]
CaPtO3 �.���/�.���/�.��� � ���� [���, ���]
CaRhO3 �.���/�.���/�.��� �� - �� ���� - ���� [���, ���]
NaIrO3 �.���/��.���/�.��� �-� ���� [���]
Fluorides ABF3
NaMgF3 �.���/�.���/�.��� �� - �� ��� - [���]
NaZnF3 �.���/��.���/�.��� �� - �� ��� - [���]
NaFeF3 �.���/��.���/�.��� � ��� - [���]
NaCoF3 �.���/��.���/�.��� �� - �� ��� - [���, ���]
NaNiF3 �.���/��.���/�.��� �� - �� ���� - ���� [���, ���]
NaMnF3 �.���/�.���/�.��� � � [�]
In this chapter we present the results on the structural, electronic and magnetic
properties of NaMnF3 as function of applied hydrostatic pressure. The results are
presented as follows. We first determine the structural ground state structures
over a wide range of pressure conditions and we predict a transition from the
Pnma ground state to a pPv phase that has not been previously explored exper-
imentally. We analyse the pPv phase transition through a vibrational character-
ization and a structural study. We then focus on the analysis of the magnetic
behaviour with a special interest on the non-collinear magnetic ground states
allowed in the pPv. We found a high-spin ground state associated with a large
ferromagnetic canting that has not been reported in any of the previous pPv
�.� computational details ��
transition observed in other compounds. The validity of the results are discussed
with respect to the values of the Coulomb, U and exchange, J parameters within
the DFT+U approach.
In Tables �.� and �.� we report the calculated silent, Raman and infra-red
(IR) phonon modes as well as the mode Grüneisen parameters ( i ) of the relaxed
Pnma ground state. The irreducible representation of the Pnma perovskite phase
is =�Au ��B1u �B2u ��B3u �Ag �B1g �B2g �B3g . Both the
Raman (Ag , B1g , B2g and B3g labels) and IR modes (B1u , B2u and B3u labels)
are in good agreement with the available experimental reports [���, ���]. We
note that a very soft polar B2u mode is present with a frequency of �� cm-1 .
This indicates that the Pnma phase is close to a ferroelectric phase transition
and this mode is responsible for the epitaxial strain induced ferroelectricity of
NaMnF3 as we reported in the Refs. [�, �] and explain in detail in Chapter
�. From the Grüneisen parameters reported in Table �.�, we observe that two
phonon modes are particularly sensitive to a change in volume: the soft B2u
mode at �� cm-1 with = �.��� and the Ag mode at �� cm-1 with =
�.���, both can be easily detected from Raman experiments. Phonon modes,
corresponding to a Au symmetry (at �� cm-1 and ��� cm-1 ), have negative
and correspond to the modes that become soft as the volume decreases, which
indicates that probably these two phonon modes are close to a nearest neighbour
force instability. Besides, these modes can contribute to a competition between
positive and negative thermal expansion of the material, depending on how the
phonon modes are populated [���].
In regards to the magnetic ordering, the total energy differences between the
ferromagnetic (FM), C- and A-type antiferromagnetic (AFM) with respect to
G-type AFM ordering in the Pnma structure are respectively ��.��, ��.�� and
��.�� meV/u.c. The lowest energy is thus given by the G-type AFM ordering,
in agreement with the experimental measurements that establish a G-type AFM
ordering for the Pnma phase of NaMnF3 [���]. All of that confirms that our
choice of parameters for the DFT calculations (GGA PBEsol and U = �.� eV) of
NaMnF3 gives good results when comparing to the experimental measurements.
Table �.�: Calculated Raman, Infra-red (IR) and silent modes frequencies (cm-1 ) of
the Pnma phase of NaMnF3 (at zero pressure). In brackets are presented the
experimental values of the Raman [���] and the IR [���] modes.
Raman modes
Ag B1g B2g B3g
�� (��) �� (—–) ��� (���) �� (��)
��� (���) ��� (���) ��� (���) ��� (—–)
��� (���) ��� (���) ��� (���) ��� (—–)
��� (���) ��� (���) ��� (���) ��� (���)
��� (���) ��� (���) ��� (—–) ��� (—–)
��� (���) — ��� (���) —
��� (���) — ��� (—–) —
Silent IR modes
Au B1u B2u B3u
�� (—–) ��� (���) �� (—–) �� (���)
��� (—–) ��� (���) ���(���) ��� (���)
��� (—–) ��� (���) ��� (���) ��� (���)
��� (—–) ��� (���) ��� (—–) ��� (���)
��� (—–) ��� (���) ��� (—–) ��� (���)
��� (—–) ��� (���) ��� (—–) ��� (���)
��� (—–) ��� (���) ��� (—–) ��� (���)
��� (—–) ��� (—–) — ��� (—–)
— ��� (—–) — ��� (—–)
Going from the lowest to the highest pressure, we found the following low en-
thalpy phases. Between -�� GPa and -�.� GPa the hexagonal P�3 /mmc phase is
the most stable structure. Between -�.� GPa and -�.� GPa we found the P�1 /c
structure to be the most stable. Then between -�.� GPa and �.� GPa the or-
thorhombic Pnma structure is the ground state and thus passing through the
atmospheric pressure as expected from the experimental results. At higher pres-
sures (above � GPa, up to �� GPa) the orthorhombic Cmcm pPv structure is the
lowest energy phase. Then, we found that NaMnF3 presents a phase transition
at relatively low pressure (� GPa) from the Pnma perovskite phase to the Cmcm
pPv phase.
In Fig. �.�a we represent a schematic view of the Pnma and the pPv phases.
It is clear that the structural transition between the Pnma and the pPv phase is
reconstructive since a breaking of Mn–F–Mn bonds is needed to connect the two
�� prediction of the p p v transition in n a m n f 3
Table �.�: Calculated mode Grüneisen parameters ( i ) for the NaMnF3 Pnma phase
at and around atmospheric pressure.
Ag B1g B2g B3g Au B1u B2u B3u
�.��� �.��� �.��� �.��� -�.��� �.��� �.��� �.���
�.��� �.��� �.��� �.��� -�.��� �.��� �.��� �.���
�.��� �.��� �.��� �.��� �.��� �.��� -�.��� �.���
�.��� �.��� �.��� �.��� �.��� �.��� �.��� �.���
�.��� �.��� �.��� �.��� �.��� �.��� �.��� �.���
�.��� — �.��� — �.��� �.��� �.��� �.���
�.��� — �.��� — �.��� �.��� �.��� �.���
— — — — �.��� �.��� — �.���
— — — — — �.��� — �.���
(a)
(b)
Figure �.�: (a) Total energy of NaMnF3 as a function of the volume per atom for all
different crystal structures. Cmcm and Pnma structures show the lowest
energy values for lower volumes. (b) H/atom as a function of external
pressure in NaMnF3 structures. Pnma to Cmcm structural transition can
be seen at the critical pressure of � GPa.
(a)
P
(b)
Intensity [a.u.]
Pnma
Cmcm
10 20 30 40 50 60 70 80
!-2! [deg]
(c)
Figure �.�: (a) Schematic view of the phase transition from Pnma to Cmcm. Na, Mn and
F atoms are represented in green, violet and grey colors respectively. The
breaking of bondings for the transition from corner-shared to edge-shared
octahedra can be appreciated. (b) Simulated XRD patterns for Pnma and
Cmcm structures. The room temperature Pnma phase is in good agreement
with the experimental reports (open triangles) [���, ���, ���]. (c) Octahe-
dra tilt angle as a function of pressure in NaBF3 . For B = Zn, Ni, Co and
Mg experimental results are taken from Yusa et al. [���]. All of these com-
pounds present the transition from the Pnma to the Cmcm structures at a
tilting angle close to �� . Experimental results for NaMnF3 at low pressure
were taken from Katrusiak et al. [��] Full agreement between experimen-
tal (filled squares) and our theoretical calculations (empty squares) can be
appreciated.
�.� study of NaMnF 3 under hydrostatic pressure ��
should occur according to the O’Keeffe rule, which is indeed in agreement with
our BM-EOS prediction of � GPa as the critical pressure.
Our theoretical findings for the transition between the Pnma and the Cmcm
pPv phase do not follow the recent experimental report of Akaogi et al. [���]. In
their pressure measurements at high temperature (T = ���� K), Akaogi et al. do
not observe a transition from the Pnma to the pPv structure up to �� GPa but
found, instead, a MnF2 + Na3 Mn2 F7 phase decomposition around � GPa. We
note that Katrusiak et al. also report a transition from the Pnma to the cubic
phase Pm3̄m at Tc = ��� K and ambient pressure [��] for this compound, which is
not observed by Akaogi et al. at very high temperature (higher than Tc = ���K).
We suggest that at the very high temperature where Akaogi et al. performed their
measurements, it is indeed possible that a phase decomposition appears as they
do not observe the pPv transition. We also remark that at such high temperature,
the O’Keeffe rule is not fulfilled since the crystal is cubic or with titling angle
smaller than the required critical value of �� to observe the Pnma to pPv phase
transition under pressure, therefore, our calculations are still valid at low temper-
atures. Following the suggested experimental analysis and in order to understand
the decomposition phases, we performed calculations for Na3 Mn2 F7 and MnF2
compounds. The first stoichiometry is the one of the Ruddlesen-Popper (RP)
phases [���] with An-1 A20 Bn X3n+1 formula, which in our particular case corre-
sponds to n = � and A 0 = A = Na, B = Mn and X = F. To our knowledge, there
is no previous report of the sodium-manganese based compound in the RP phase.
Therefore, all reported phases for similar compounds (oxides [���] and fluorides
[���]) were tested and the energy minimum per structure was analysed. Neverthe-
less, none of our simulated XRD patterns of the considered RP phases match the
one reported by Akaogi et al. [���], while for the pristine MnF2 compound, we
obtain good agreement with the experiments [���]. We are thus confident that
the crystal phase transition of NaMnF3 from the perovskite to the pPv phase
occurs at the predicted � GPa pressure. In any case, more experimental measure-
ments would be required in order to clarify the pressure phase transition that
takes place in NaMnF3 at lower (room) temperature.
We now focus on the vibrational characterization of the pPv phase at � GPa.
In Table �.� we report the computed phonon mode frequencies and the mode
Grüneisen parameters ( i ) in the Cmcm pPv structure with the lattice parame-
ters relaxed at � GPa (also reported in Table �.�). The irreducible representation
at the zone center point is: �Au �B1u �B2u �B3u �Ag �B1g B2g
�B3g . Two silent modes (Au label), sixteen IR active modes (B1u , B2u and
B3u ) and twelve Raman active modes (B1g , B2g and B3g ) are thus expected. As
we can see in Table �.�, all the calculated modes at the point of the BZ are
positive, indicating a vibrational local stability of the pPv phase at �GPa. The
lowest frequency mode is a silent mode Au at �� cm-1 and the lowest polar mode
is a B1u mode with a frequency of �� cm-1 . The obtained Grüneisen parameters
are important in the vibrational identification of the transition related to Raman
experiments under pressure. In general, they are smaller than the one obtained
�� prediction of the p p v transition in n a m n f 3
Table �.�: Computed Raman, Infra-red and silent modes of NaMnF3 in post-perovskite
Cmcm phase as well as the mode Grüneisen parameters ( i ) at around ��
GPa.
Raman, IR and silent modes [cm-1 ]
Ag B1g B2g B3g Au B1u B2u B3u
��� ��� ��� ��� �� �� ��� ���
��� ��� — ��� ��� ��� ��� ���
��� ��� — ��� — ��� ��� ���
��� — — ��� — ��� ��� —
— — — — — ��� ��� —
Mode Grüneisen parameters ( i )
�.��� �.��� �.��� �.��� �.��� -�.��� -�.��� �.���
�.��� �.��� — �.��� �.��� -�.��� �.��� �.���
�.��� �.��� — �.��� — �.��� �.��� �.���
�.��� — — �.��� — �.��� �.��� —
— — — — — �.��� �.��� —
in the Pnma ground state structure. This can be explained from an interatomic
force stiffness larger in the pPv phase than in the Pnma phase. Even though
there are some phonon modes with negative Grüneisen parameter (B1u modes
at �� cm-1 and ��� cm-1 and the B2u mode at ��� cm-1 ), they are not large
enough to have an effect on the thermal expansion such as the one observed in
↵-ZrW2 O8 [���].
We also report in Fig. �.� the evolution of the Raman and IR modes as a
function of pressure from � GPa to �� GPa. After the transition, the system goes
to a reduced crystallographic unit cell such as the number of modes is two times
smaller in the Cmcm pPv phase than it is in the Pnma phase. We note that
for both Pnma and pPv structure, all the phonon mode frequencies harden when
increasing the pressure. This is also valid for the soft B2u polar mode of the Pnma
structure that hardens as the pressure increases and thus in agreement with the
usual behaviour of perovskites in the presence of the AFD distortions, which are
inimical to the ferroelectric distortions and are enhanced with pressure.
Figure �.�: Vibrational behavior as a function of isotropic pressure for Pnma and Cmcm
pPv. The open and filled symbols are the Raman and the Infra-red modes
respectively.
group [���, ���] and the Mn atoms are placed at the 4c Wyckoff position with
coordinates (�,�,�) for Mn1 , (�,�, 12 ) for Mn2 , ( 12 , 12 ,�) for Mn3 and ( 12 , 12 , 12 ) for
Mn4 . From these � magnetic cations, we can define � collinear magnetic orderings
depending on the sign of the magnetic moments of (Mn1 , Mn2 , Mn3 , Mn4 ): F
(+, +, +, +) for the FM ordering, A (+, +, -, -) for the A-type AFM ordering,
C (+, -, +, -) for the C-type AFM ordering and G (+, -, -, +) for the G-
type AFM ordering, where + and - correspond to spin up and down respectively.
Of course, the ferrimagnetic ordering (+, +, +, -) could be considered but it
is unlikely to appear in real conditions because of the symmetry of the system
and the superexchange interaction [���, ���]. For the non-collinearity we have
to consider the three Cartesian directions x, y and z for each magnetic moment
which gives �⇥� possible states: Fx , Fy , Fz , Ax , Ay , Az , Cx , Cy , Cz , Gx , Gy and
Gz . Then, taking into account the fact that the spins are axial vectors, one has
to apply the symmetry operations of the D2h point group on each of the non-
collinear magnetic orders defined above and see how they transform. The result
of this approach applied on the pPv phase of NaMnF3 is summarised in Table
�.�.
�� prediction of the p p v transition in n a m n f 3
Table �.�: Magnetic ordering allowed in the Cmcm structure according to the D2h
point group [���, ���]. The transformation of each magnetic ordering under
each symmetry operation is labeled as +� and -�, which indicates when the
ordering is invariant under the related transformation (i.e. +�) or when it is
reversed under the application of the symmetry operation (i.e. -�).
Magnetic Orderinga I �z �y �x -� mz my mx
Fx , Ax � -� -� � � -� -� � B3g
Fy , Ay , Cz , Gz � -� � -� � -� � -� B2g
Fz , Az , Cy , Gy � � -� -� � � -� -� B1g
Cx , Gx , � � � � � � � � Ag
As we can see from Table �.�, the Fx and Ax orders have to the same char-
acter. This means that in an energy expansion with respect to these two order
parameters the second order term Fx · Ax is invariant by symmetry. This implies
that if one of the two order parameters condenses in the structure, the second
one can develop as well but as a secondary or a slave order (pseudo-proper �
[���]). This secondary order parameter usually appears as a spin canting in the
real system and its intensity will depend on the amplitude of the Fx · Ax coupling.
Going through the other magnetic orders, in Table �.� we also observe that Fy ,
Ay , Cz and Gz belong to the same B2g character, Fz , Az , Cy and Gy belong
to the B1g character and Cx and Gx belong to the Ag character. Interestingly,
whichever is the main magnetic order we observe in the NaMnF3 pPv phase, a
canting of the spins is always possible. Additionally, we also remark that three
of the four possible characters to which belong all the possible magnetic orders,
three of them allow for ferromagnetism and so have the potential to exhibit weak
FM.
As discussed above, NaMnF3 compound in the Pnma phase has a G-type AFM
ordering [��, ���, ���]. In Table �.� we report the computed Mn–F–Mn distances
and angles of the Pnma phase at room pressure and of the pPv phase at � GPa.
We note that in the Pnma structure, the Mn–Mn bonding in the three directions
have similar Mn–F–Mn angles (about ��� ) and distances (about � Å) while
their are much more anisotropic in the pPv phase (angles of �� and ��� and
one non-existent angle due to the octahedral edge sharing of the structure with
distances of �.�� Šand �.�� Šin the x - z plane and �.�� Šin the y direction).
Having angles closer to ��� in the Pnma explains the G-type AFM ordering
� A pseudo-proper transition appears when a secondary order parameter (Y ) has the same sym-
metry as the main transition order parameter (X) such as they coupled linearly (first order
coupling) in a free energy expansion: F=aX2 +bY 2 +cXY 2 . In the present case the ferromagnetic
(secondary) order parameter transforms as the main antiferromagnetic order parameter, giving
rise to the so-called weak-ferromagnetism.
�.� magnetic properties of the ppv phase ��
Table �.�: Comparison of the distances and the bond angles along the Mn–F–Mn bonds
between the Pnma and the Cmcm structures of NaMnF3 .
Distance and Angle Pnma Cmcm
[���] �.���� [���] �.����
dMn-Mn [Å] [��-�] �.���� [���] �.����
[���] �.���� [���] �.����
[���] ���.�� [���] ��.��
Mn-F-Mn [deg] [��-�] ���.�� [���] ���.��
[���] ���.�� [���] —
since with such bond angle the superexchange interaction between the Mn-d5
spins favours an AFM alignment [���, ���, ���]. In the pPv case, we can see that
there is no direct bonding between the Mn atoms along the y direction such as
we can expect a weak spin coordination between the spins along the y direction
and one would expect a quasi-�D magnetic behaviour in the pPv structure. In
the plane perpendicular to the y direction we have two types of bonding giving
rise to an angle close to ��� along the z direction and to an angle close to ��
along the x direction. On the top of that, we also note that the �� bonding
along the x direction is given through two bonds since it is the direction where
the octahedra share their edge. The superexchange interaction should thus drives
strong FM ordering along the x axis and AFM ordering along the z axis and a
weak interaction along the y axis. Interestingly, the magnetic ordering that fulfils
these superexchange rules is the C-type AFM and it is indeed the lowest energy
collinear magnetic ordering that we found in our calculations (see Fig. �.�b).
Therefore, we can only consider the set of non-collinear calculations containing
the C-type AFM: (Cz , Fy , Gz , Ay ), (Cy , Fz , Az , Gy ) and (Cx , Gz ) as predicted
by group theory. We can also expect the Cx case to be unlikely observed since
it forces the spins to be antiferromagnetically aligned along the x direction while
a strong FM interaction is expected from the double �� bond angle connecting
the Mn atoms in this direction.
Relaxing the spin structure within the non-collinear regime and within the
C-type AFM ordering, we found that the lowest energy is obtained for the Cy
orientation with a large Fz component along the z direction (no canting along
the x axis is observed). This magnetic ground state is however dependent on
the U and J values of the GGA+U exchange correlation functional. In Fig. �.�a
we report the energy difference between the Cy Fz ordering and the other lowest
energy cases Cz Fy , Cz Ay , and Gy Az versus the J value (fixing J = � eV). Here
we remark that the Cy Fz , Cz Fy , Cz Ay , and Gy Az states are very close in energy
such as we can expect an easy magnetic phase transitions with small perturbations
such as an applied magnetic field. Such a transition under a magnetic field has
been observed experimentally in the pPv phase of NaNiF3 where a field parallel
�� prediction of the p p v transition in n a m n f 3
CyFz
CzAy
(a) (b)
Figure �.�: (a) Energy differences between the different magnetic orderings of the pPv
phase of NaMnF3 versus the U parameter. (b) Sketch of the C y F z and the
C z A y magnetic orderings in the pPv phase of NaMnF3 . For U < �.� eV
the system exhibits a AFM+FM magnetic state with a high ferromagnetic
canting. For U > �.� eV, the system transit toward a canted AFM state
with no weak FM.
is difficult to predict what would be the best J value for NaMnF3 without any
experimental feedback. Previous studies about the DFT+U,J parameter report
that low J values for d5 orbital filling are usually the best values [���] such as
the large values of the Fz canting we obtain at low J (close to � µB /atom) in the
pPv phase of NaMnF3 would be the optimal value.
In order to estimate the error bar of our calculations, we report in Table �.�
a list of known magnetic pPv phase in oxide and fluoride perovskites with their
related magnetic orderings, d-orbital filling and S quantum number. We note
that in all oxide cases as well as in NaNiF3 , the ground state spin configuration
always corresponds to the low spin state. According to our calculations, only
NaMnF3 has the high spin configuration as ground state. In order to verify if this
is due to the DFT parameters used in our calculations or if this has a physical
meaning, we performed calculations of the pPv phase of NaNiF3 and CaRhO3
with U =�.� eV. We found a low spin ground state in both NaNiF3 (S=�) and
CaRhO3 (S= 12 ) and so in good agreement with the experimental results. Going
further, we also estimated the non-collinear spin ground state of NaNiF3 and
CaRhO3 in their pPv phase. As expected from the symmetry analysis, we observe
canted magnetic ground state in the two systems. The ground state magnetic
ordering of NaNiF3 is found to be Az with a Gy canting. Then, we did not find
any weak ferromagnetic spin canting in NaNiF3 in agreement with experimental
data. We note that experimentally a transition towards an AFM state with a
ferromagnetic canting is observed in the pPv phase of NaNiF3 when applying
a small magnetic field. The direction of the field induced canted FM spins in
NaNiF3 are the same than the one we report in the pPv phase of NaMnF3 without
magnetic field. This shows that the high-spin large FM canting we report in
NaMnF3 might be metastable in NaNiF3 such as a small magnetic perturbation
can switch to this state. More theoretical and experimental analysis would be
welcome to understand this transition. In the case of CaRhO3 we found Cy with
a Fz canting as magnetic ground state.
�� prediction of the p p v transition in n a m n f 3
Table �.�: Magnetic orderings allowed in the Cmcm structure for ABX3 systems with U
= �.� eV and J = � eV. The large magnetic canting for fluorides compounds
can be noted. The canting angle (') is measured with respect to the ab-plane
on the pPv phase. For NaNiF3 , the two possible magnetic states Cz Fy and
Cy Fz were computed.
Compound Fi [µB /atom] ' [deg] Magnetic Ground-state
CaRhO3 �.�� �.� Cy Fz
NaNiF3 �.�� (�.��) ��.� (�.�) Cz Fy (Cy Fz )
NaMnF3 �.�� ��.� Cy Fz
We found that when considering the cases with ferromagnetic spin canting, the
canting angle is much smaller in CaRhO3 (�.� ) and in NaNiF3 (�.� ) than it
is in NaMnF3 (��.� , see Table �.�). If we look at the amplitude of the canting
per atom we observe a very large amplitude for NaMnF3 (�.�� µB /atom) due to
the high spin configuration while it is much smaller in the case of NaNiF3 (�.��
µB /atom) and of CaRhO3 (�.�� µB /atom). This large canting can be the result
of a rather large Dzyaloshinsky-Moriya interaction [���, ���] or a large single-ion
anisotropy [���]. A more detailed analysis in order to identify the exact origin of
this large spin canting can be found in Chapter �. These previous calculations
allow us to be more confident about the unique high spin ground state of NaMnF3
with a large ferromagnetic spin canting in its pPv phase.
We also report in Fig. �.� the local density of states (LDOS) of the pPv phase
of NaMnF3 . The LDOS shows an overlap between the Mn-3d and F-2p states,
suggesting a covalent contribution between the Mn and F atoms due to the p - d
hybridization, which is a fundamental requirement for the superexchange interac-
tion. Moreover, we observe that the fluorine states are lower in energy then the
Mn-d orbitals, which corresponds to the expected bonding-antibonding picture
of covalency between magnetic cations and anions. Close to the Fermi level, no
contribution of Na states to the DOS is observed such as mostly Na size effects
are expected within the structure, a key-point parameter that is at the origin
of the geometric ferroelectric instabilities observed in cubic fluoroperovskites [�].
The dxy , dxz , dyz , dz2 and dx2 -y2 orbital occupancy gives the � orbitals filled
below the Fermi level for a spin channel and the � orbitals empty above of the
Fermi level for the other spin channel, which confirms the high spin electronic
configuration of Mn2+ in the pPv phase of NaMnF3 .
�.� Conclusions
We have studied, by means of first-principles calculations, the high pressure phase
diagram of NaMnF3 . We predicted a phase transition from the Pnma distorted
perovskite to a Cmcm post-perovskite phase at a critical pressure of � GPa, which
�.� conclusions ��
Table �.�: Magnetic ordering in fluorides and oxides post-perovskite systems, AFM: an-
tiferromagnetic, DM: diamagnetic, FM: ferromagnetic, CWP: Curie-Weiss
paramagnetic. All pPv compounds exhibit a low-spin magnetic configuration
but NaMnF3 that according to our calculations is found to be in the high-spin
state.
Compound Magnetic ordering Ref.
pPv Oxides ABO3
CaIrO3 quasi-�D AFM S= 12 Ir4+ : �d5 (LS) [���]
CaRuO3 quasi-�D AFM S=� Ru4+ : �d4 (LS) [���, ���]
CaPtO3 DM S=� Pt4+ : �d6 (LS) [���, ���]
CaRhO3 AFM+weak-FM S= 12 Rh4+ : �d5 (LS) [���]
NaIrO3 CWP S=� Ir5+ : �d4 (LS) [���]
pPv Fluorides ABF3
NaFeF3 AFM S=� Fe2+ : �d6 (LS) [���]
NaCoF3 AFM — Co2+ : �d7 (LS) [���]
NaNiF3 AFMa S=� Ni2+ : �d8 (LS) [���, �]
NaMnF3 AFM+FM S= 52 Mn2+ : �d5 (HS) [�]
a NaNiF3 compound present a magnetic transition from AFM to AFM+FM ordering at applied
external magnetic field [���].
Figure �.�: Projected local density of states (LDOS) of the Cmcm phase of NaMnF3 .
A regular F-p bonding Mn-d antibonding is observed with small Mn:3d –
F:2p overlaping DOS is responsible for the small superexchange interaction
present in NaMnF3 .
�� prediction of the p p v transition in n a m n f 3
is relatively low with respect to the similar phase transitions observed in other
compounds. As expected from the tilting criteria[���, ���], this phase transition
fits very well with the octahedra tilting angle requirements where the transition
to the pPv is observed when the octahedra tilting angle reach the critical value
of �� with an initial ambient pressure tilting angle of �� . Unfortunately high
pressure experiments have been performed at very high temperature by Akaogi
et al. et al.[���] where the tilting criteria requirement is not fulfilled and also by
Katrusiak et al. [��] who unfortunately did not go beyond � GPa such as none
of them observed the pPv transition. Because of the partially filled d orbitals of
Mn2+ cation we also determined the magnetic ground state of NaMnF3 in its pPv
phase. Our calculations predicted a high spin ground state of the d5 spins with
a magnetic moment of �µB on each Mn. This high spin structure is particularly
interesting and unique since all the other known magnetic compounds presenting
the pPv phase shows a low spin ground state. We predicted an AFM magnetic
ordering ground state of the C-type with the spins lying along the y-direction
and with a surprisingly large ferromagnetic spin canting along the z-direction.
This makes NaMnF3 of particular interest for fundamental and geophysical stud-
ies since it has this unique property of high spin magnetic ground state with
potentially large ferromagnetic spin canting and because it can be synthesized
at relatively small pressure. The fact that the fluoride post-perovskites can be
maintained when the pressure is released makes NaMnF3 even more interesting
to study.
It is important to remark that very recently, another transition at higher pres-
sures was reported and named as a post-pPv [���]. In this phase, the structure
transforms to a face-shared octahedra with a substantial compression along the
crystallographic a-axis into a ppPv-Pnma phase —see Fig. � in Ref [���]—. More-
over, it has been found that some oxide compounds transforms directly from the
room conditions Pnma to the post-pPv under isotropic pressures above �� GPa.
Interestingly, the magnetic ordering in the ppPv materials, with magnetically
active A- and B-sites, shows magnetic “double” one-dimensional channels in ad-
dition to a triangular-like coordination of the B-site generating a geometrically
frustrated magnetic system.
Therefore, we hope that our results on the pPv phase and the reported ppPv
transition will motivate the experimentalist to carry out more studies of NaMnF3
in order to confirm the pPv phase transition at low pressure and low tempera-
ture and to check its potential for high spin structure with a sizeable canted
ferromagnetism.
NON-COLLINEAR MAGNETISM IN
4
CaRhO3 AND NaNiF3
POST-PEROVSKITES
��
�� non-collinear magnetism in c a r h o 3 and n a n i f 3 post-perovskites
Figure �.�: Simulated XRD patterns with K ↵ :Cu for CaRhO3 (above) and NaNiF3
(below) in the Cmcm phase at atmospheric pressure are observed.
ments. In Table �.� we present the calculated Raman, Infra-Red (IR) and silent
mode frequencies of the pPv phase of CaRhO3 and NaNiF3 at ambient pressure.
Even if the pPv is not the ground state structure at room pressure (The Pnma
phase is the ground state), we do not observe any phonon instabilities confirm-
ing that this phase is a metastable phase (The Pnma phase has a lower energy
than the pPv phase) and as observed experimentally when releasing the pressure
[���, ���]. The stability of the pPv phase at room conditions is of great interest
due to the possibility to explore the properties of the pPv structure and its use
for technological applications.
Table �.�: Vibrational characterization modes of CaRhO3 and NaNiF 3 Cmcm phase at
atmospheric pressure with the representation: = �Au �B1u �B2u
�B3u �Ag �B1g B2g �B3g .
Raman , IR and silent (Au ) modes [cm-1 ]
CaRhO3
Ag B1g B2g B3g Au B1u B2u B3u
��� ��� ��� ��� ��� ��� ��� ���
��� ��� — ��� ��� ��� ��� ���
��� ��� — ��� — ��� ��� ���
��� — — ��� — ��� ��� —
— — — — — ��� ��� —
NaNiF3
��� ��� ��� ��� �� ��� ��� ���
��� ��� — ��� ��� ��� ��� ���
��� ��� — ��� — ��� ��� ���
��� — — ��� — ��� ��� —
— — — — — ��� ��� —
the pPv phase when comparing with the regular distorted perovskite structures
as the Pnma phase.
Figure �.�: The non-collinear magnetic orderings in the pPv Cmcm symmetry labeled as
Cz Fy , Cz Ay , Cy Fz and Gy Az are presented in the figure. Two of them with
pure AFM (Cz Ay and Gy Az ) behavior and two AFM+weak-FM (Cz Fy and
Cy Fz ) can be observed (Fig. elaborated with the VESTA software [���]).
can neglect unlikely states and thus reducing the possible magnetic states to
Cz Fy , Cz Ay , Cy Fz and Gy Az , as summarized in Table �.� and depicted in Fig.
�.�.
Table �.�: Magnetic orderings allowed in the Cmcm pPv structure according to the D2h
point group [�]. Here F is ferromagnetic ordering and A - C - and G - are the
different AFM magnetic orderings where the subindex indicate the direction
along the particular Cartesian axis.
Magnetic Ordering Character
Fx , Ax B3g
Fy , Ay , Cz , Gz B2g
Fz , Az , Cy , Gy B1g
Cx , Gx , Ag
According to our total energy calculations, the ground state corresponds to the
Cz Fy and Gy Az magnetic states for CaRhO3 and NaNiF3 respectively (see Table
�.�). It can be noted that these magnetic ground states as well as their canting
amplitude are in good agreement with available experimental data. In Table �.�,
we also report the experimental values of the magnetic state of NaNiF3 under
an applied magnetic field and we compare it with a calculation made with a
higher energy configuration. The later corresponding to the magnetic ground
state observed experimentally where we get a good agreement. We remark that
the atomic magnetic moment in the NaNiF3 case is almost three times larger
that in the oxide CaRhO3 . This is due to the occupation of the B-site magnetic
cation that drives a S = � state in NaNiF3 and S = 12 in CaRhO3 .
Table �.�: Non-collinear ordering and moments in CaRhO3 and NaNiF3 . Experimental
and our theoretical calculations are in good agreement. The magnetic mo-
ments are in µB /atom for each magnetic B-site cation. The induced magnetic
ground state by applied magnetic field (Hext ) in NaNiF3 is also presented.
Compound Ordering C, G F, A Fi exp.
CaRhO3 Cz Fy �.�� �.�� �.�� [���]
NaNiF3 Gy Az �.�� �.�� — [���]
NaNiF3 (Hext > � kOe) Cz Fy �.�� �.�� ⇡ �.� [���]
the AFeO3 systems). Here, we propose to analyse the origin of the non-collinear
magnetism in the pPv phase of NaNiF3 and CaRhO3 by decomposing the mag-
netic interactions into its isotropic exchange (J ), DM and SIA components. To
that end, it is convenient to start by defining the magnetic Heisenberg Hamilto-
nian (H) as follows:
Da Dc
Jab
Ja
Jc
Figure �.�: Schematic representation of the atomic sites and magnetic interactions in
ABX3 (A, B and X -sites in blue, yellow and red respectively). Exchange
constants labeled as Ja , Jc and Jab and DM interactions represented as Da
and Dc equivalent to D12 and D13 respectively are shown. (Fig. elaborated
with the VESTA software [���]).
Eij ij ij ij
"! + E# - E" - E#!
Dk
i 0j 0 = , (�.�)
4S2
Ej" + Ej# - Ek k
" - E#
Ai = , (�.�)
2S2
The values for the Jij constants (see Table �.�) demonstrate a ferromagnetic
interaction along the x-axis with Ja > �. Additionally, we can observe a strong
AFM coupling along the z-axis and a weak interlayer interaction, in agreement
with the non-collinear magnetism results reported in Table �.�. In the CaRhO3
compound, Jab >�, which corresponds to a FM interlayer interaction as expected
from the Cz Fy magnetic ground state. In NaNiF3 Jab < 0 which confirms the
Gy Az ordering and in agreement with the group theory analysis. It is important
to remark that the magnetic constants in CaRhO3 are at least one order of
magnitude larger that in NaNiF3 . This can be explained by the lower covalence
in the B–F–B bonding, when compared to the B–O–B and related also to the
expected JT -distortion with Rh4+ cation as well as the fact that �d and �d
electrons are known to have higher exchange constants than the �d electrons
[���].
As demonstrated by Kim and Min [���] in d3 , d5 and d8 , B-sites, the DM
vector in perovskites can be related to the inter ionic vectors as Dij / (x̂i ⇥
x̂j ) [���], where x̂i and x̂j are unitary vectors along the B–X–B bonds. Thus,
according to the previous relation and the symmetry of the pPv compounds, the
Da must be zero and the Dc has to be along the x-axis, then: Da x = Dy = Dz =
a a
in Mn3+ (d4 ), Cu2+ (d9 ) or Rh4+ with t52g occupation), there is a change in the
orbital symmetry, breaking the mirror symmetry and creating a deviation from
the Kim and Min rule (as it happens in LaMnO3 [���]). In Table �.� we report
the components of the DM vectors for the two considered compounds. Analyzing
the DM components, we conclude that NaNiF3 , with Ni2+ , perfectly fulfill the
Kim and Min rule with a non-zero Dcx = �.��� meV. However, in the CaRhO3
case, a deviation is observed where Da x = -�.��� meV. This is probably due to the
t52g state that induces a JT -distortion, which in turn changes the orbital symmetry
and then, change the DM vectors.
According to the relation between the exchange constants and the magnitude
of the DM vector, the magnetic canting angle ' can be defined as [���]:
!
1 ~ ij |
|D
' = tan-1 (�.�)
2 Jij
if there is no SIA in the direction of the canting [���, ���]. We can estimate
the canting angle from this simple formulae and compare the results with the
�� non-collinear magnetism in c a r h o 3 and n a n i f 3 post-perovskites
relaxed values obtained through the DFT calculations. For NaNiF3 the canting
angle coming from the relaxed DFT calculations are �.�� and �.�� for Gy Az
and Cz Fy respectively, and we obtain ' = �.�� according to the calculated Dij
and Jij parameters and the use of Eq.�.�. In CaRhO3 the canting angle coming
from the fully relaxed calculations is �.�� while it is �.�� from Eq.�.�. The large
discrepancy between the relaxed and estimated canting angle in CaRhO3 can be
related to the competition between JT -distortions, orbital symmetry [���] and
electronic Coulomb localization [���] through the SIA. The contribution of SIA
to the canting angle has been first reported by Moriya in NiF2 where the DM
interaction is forbidden by symmetry [���, ���].
Table �.�: Computed magnetic exchange constants for pPv CaRhO3 and NaNiF3 . Units
in meV.
Ja Jc Jab
CaRhO3 ��.�� -��.�� �.��
NaNiF3 �.�� -�.�� -�.��
Table �.�: Dzyaloshinsky-Moriya (Dij ) vectors computed for pPv CaRhO3 and NaNiF3 .
Here, D12 ⌘ Da and D13 ⌘ Dc as shown in Fig. �.�. Units in meV.
Da Dc
Da
x Da
y Da
z Dcx Dcy Dcz
CaRhO3 -�.��� �.��� �.��� -�.��� �.��� �.���
NaNiF3 �.��� �.��� �.��� �.��� �.��� �.���
Table �.�: Computed SIA parameters for pPv CaRhO3 and NaNiF3 . Units in meV.
Ax Ay Az
CaRhO3 -�.��� �.��� �.���
NaNiF3 �.��� �.��� -�.���
In the search for new and innovative materials with exciting multifunctional prop-
erties, multiferroics and magnetoelectrics has been a keystone in condensed mat-
ter physics in the last decade [���, ���]. Several systems and compounds based
on oxide perovskites have been reported as ideal candidates. In these systems, the
presence of the ferroic orders can be induced and/or tuned by different physical
constrains such as the case of the biaxial epitaxial strain [���] in thin films and
superlattices form [���]. In this search, one of the main challenges is to induce a
ferroelectric order in magnetic and non-magnetic perovskites. Then, taking the
biaxial epitaxial strain as a control parameter, the condensation of a stable fer-
roelectric state has been successfully reported in CaTiO3 [���], SrTiO3 [��] and
CaMnO3 [���, ���] the first two of them paramagnetic materials at their ground
state and the later antiferromagnetic (AFM) with a weak ferromagnetism (weak-
FM) due to a magnetic canted structure [���]. Here, the strain couples with the
softest mode and thus, induce a polar displacement achieving with this a phase
transition to a polar state in a already AFM structure [���] and thus exceeding the
so-called d0 -ness rule that is expected to prevent the formation of a FE phase in
magnetic perovskites [���]. However, the possibility for new materials with these
multifunctional properties in new stoichiometries is still evasive. Recently, fluo-
rides have been reported as possible candidates for a multiferroic behavior with
promising properties [��]. Nonetheless, none of the reported fluoride candidates
belong to the most claimed perovskite family.
��
�� induced multiferroic and magnetoelectric behavior in n a m n f 3
phase approach [��] as implemented in the VASP code. The ME coupling was
obtained by computing the spontaneous polarization as a function of the applied
Zeeman magnetic field as implemented by Bousquet et al. [��] within the LDA
approach —for more details see Chapter �—. Besides, we have evaluated the ME
response against the U and J parameters of the DFT+U as shown later.
Figure �.�: Ferroelectric instability T1u at the cubic Pm3̄m as well as the frequency
values for the softest B2u mode in the Pnma NaBF3 (B = Mn, V, Zn, and
Ni) compounds.
Based on an experimental
p point of view, we suppose a cubic perovskite sub-
strate where a = c = 2ao . Thus, the �% strain value corresponds to ao =
ai +ci
2 where ai and ci are the unstrained relaxed cell parameters of the Pnma
phase defining " = a-aao 100%. After the analysis and calculations, the only flu-
o
Fig. �.�). Thus, a detailed research in the later was carried out and all the results
are presented in the following sections.
Chapter � and Ref. [�]). Thus, the competition between ferroelectric, antipolar
displacements and AFD distortions at the ground state are present in this fluoride
and the balance between those is really delicate. Despite the presence of the FE
mode, as discussed above, the ground state belongs to a non-polar orthorhombic
Pnma structure. However, as we found, the ferroelectric instability in fluorides
present a totally different origin, which is related to A-site geometrically driven
polar displacements [�].
In FE oxides it is well known that in some particular cases, the epitaxial strain
can induce ferroelectricity in the Pnma phase [���] or enhance it in compounds
with different ground states and symmetries such as PbTiO3 or BiFeO3 [���]. In
NaMnF3 the lowest FE mode has the B2u symmetry and it is polarized perpen-
dicular to the biaxial constraint (the orthorhombic b-axis). Unexpectedly, this
mode becomes unstable whatever the value of the epitaxial strain in compression
or expansion regime as presented in Fig. �.�a. By freezing the unstable mode at
and performing the structural relaxation, the new ferroelectric ground state is the
polar Pna�1 space group as it can be observed from the double well energy curve
presented in Fig. �.�c. Then, the ground state of the epitaxially constrained struc-
ture is always ferroelectric and thus multiferroic due to the G-type AFM NaMnF3 .
This finding is completely different from oxides in which some of them, originally
paraelectric, can only reach the ferroelectric state for strain values above �.�%
such as in CaTiO3 [���], �.�% in SrTiO3 [��] and �.�% in CaMnO3 [���]. Ad-
ditionally, for negative or positive strain, the polarization is enhanced up to a
value of �� µC·cm-2 at +�% and � µC·cm-2 at -�% respectively, and here again,
a strong difference with previous observations and reports in ferroelectric oxides
[���]. We note also that the polarization develops in the direction where the anti-
polar motions of the X+ 5 mode are absent and thus, where there is still possibility
for the A-site to move in the remaining flat energy landscape direction obtaining
a FE stable state along b-axis with an anti-polar character in the ac-plane.
�.� strain-induced ferroelectricity in NaMnF 3 ��
Figure �.�: a) Epitaxially strain induced B2u instability in Pnma NaMnF3 in which
a nonlinear coupling is observed, imaginary frequencies are represented as
negative values by notation. b) Eigendisplacements of the polar B2u along
b-axis (⌘b ) for each site (Na in green squares, Mn in violet circles and
F? , Fk in upper and down triangles respectively) in NaMnF3 in Pna�1 . c)
Ferroelectric mode freezing in which the transition from the paraelectric
Pnma to a polar Pna�1 phase is observed at equivalent �% in the ac-plane.
In order to understand the origin and physical phenomena within the transi-
tion and the stable FE ground state, we analyzed the eigendisplacements of the
unstable mode B2u under epitaxial strain (see Fig. �.�b). From this figure it can
be concluded first that as in oxides, the Na alkali-metal ion (in this case the
A-site) tries to annihilate the ferroelectricity when the strain goes from �%, to
-�% with eigendisplacements of �.��� and �.��� for Na respectively. However, the
antipolar motion of B-site (Mn in this case) decrease and the eigendisplacements
of this site begins to present a polar contribution to the b-axis that compensates
the decrease in contribution of the Na site (see Table �.�). This can explain the
anomalous behavior of FE-instability at compression epitaxial strain. The fluo-
rine displacements remain almost invariant for expansion or compression strain
with no significant effect on the instability behavior. The analysis of the Born
effective charges (Z⇤ in the Table �.�) reveal that values of this parameter are
�� induced multiferroic and magnetoelectric behavior in n a m n f 3
Table �.�: Born effective charges (Z⇤ , in e- ) of selected strain values ("). F? and Fk rep-
resent the Born effective charge of the F atom when it is displaced respectively
against the Mn cation and along the perpendicular direction. Eigendisplace-
ments (⌘b ) of the FE unstable mode along the b-axis in NaMnF3 . F? and
Fk have the same sense as given in Born effective charges.
Z⇤ [e- ] ⌘b [a.u.]
" [%] Na Mn F? Fk Na Mn F? Fk
Nominal � � -� -� — — — —
+�% �.�� �.�� -�.�� -�.�� �.��� �.��� -�.��� -�.���
�% �.�� �.�� -�.�� -�.�� �.��� �.��� -�.��� -�.���
–�% �.�� �.�� -�.�� -�.�� �.��� �.��� -�.��� -�.���
close to nominal ones (+�, +� and -� for Na, Mn, and F respectively) which is
related to a geometric origin of the ferroelectricity in NaMnF3 as expected [�]
instead of a charge-transfer phenomena found in most of the ferroelectric oxide
perovskites [��]. A x
Ax
Gz
Py
Fy
Ax
NaMnF3
G z Ax Fy + Py
a) b)
Figure �.�: a) Polarization along b -axis in NaMnF3 in Pna�1 symmetry. The second-
order dependence of the polarization as a function of ac -plane epitaxial
strain is observed as a result of the non-linear piezoelectricity in this per-
ovskite fluoride. b) NaMnF3 , Pna�1 structure is also presented where Na,
Mn and F ions are depicted in green, violet and grey respectively. Addi-
tionally, the non-collinear magnetic ground state as well as the polarization
vector along the b -axis are schematically depicted in the figure showing the
multiferroic character of this fluoroperovskite.
X X
Pµ = eµj "j + Bµjk "j "k , (�.�)
j j,k
⇣ e ⌘
Py = e21 + e23 - 22 "1 + ...
✓ n ◆
1 B221 + B223 B222 (�.�)
+ B211 + B233 + 2B213 - 2 + 2 "21
2 n n
The computed linear piezoelectric coefficients e21 , e23 , e22 , e16 , and e34 at
�% strain are equal to �.���, -�.���, �.���, -�.���, and �.��� C·m-2 respectively
and are not far from the ones observed in BaTiO3 [���, ���]. Nonetheless, the
computation of quadratic order components is beyond of the reach of this chap-
ter and the description and analysis of those parameters needs further studies.
Nevertheless, from the quadratic fitting observed in Fig. �.� we can extract that
the non-linear contribution related to the sum of the Bµjk component is at least
the ��% of the linear piezoelectric constants which points out the necessity for a
further study of this property in this material.
� In the Pnma structure that a and c parameters are not exactly the same and an average between
them is taken in order to apply the bi-axial strain. Then, an small plus and minus strain is
intrinsically applied to the a and c parameters respectively.
�� induced multiferroic and magnetoelectric behavior in n a m n f 3
1
Pi = ↵ik Hk + ijk Hj Hk + jik Hj Ek , (�.�)
2
�.� non-collinear magnetism and me coupling in pna� 1 NaMnF 3 ��
Table �.�: Allowed magnetic orderings in the D2h point symmetry group [���, ���]
for Pnma [���] and epitaxially-induced phase Pna�1 . Besides, the magneto-
electric couplings [���] for Pna�1 are shown depending to the non-collinear
ordering and the magnetic group.
Magnetic Pna�1
Ordering Character ME-linear ME-second order
Cx , Gy , Fz A2 ↵yz , ↵zy ijk , jik
Ax , Fy , Gz B1 ↵xx , ↵yy , ↵zz ijk , jik
Fx , Ay , Cz B2 ↵xz , ↵zx ijk , jik
Gx , Cy , Az A1 ↵xy , ↵yx ijk , jik
The magnetoelectric coupling was computed for several strain values according
to the implementation of finite magnetic field given in Ref. [��] and explained in
Chapter �. The results are condensed in Fig. �.� and in Table �.�. In Fig. �.� we
report the evolution of the electric polarization with respect to the amplitude of
the applied magnetic field along different directions. In Table �.� we report the
extracted ME coefficients at three different epitaxial strains. Our results reveal
that a sizeable non-linear ME coupling is present when the magnetic field is
applied along the y-axis (parallel to the weak-FM moment and the spontaneous
polarization). When the field is applied along the x-axis, we observe a linear ME
response along the same direction (↵xx ) and a non-linear one along the y-axis
( yxx ).
In comparison to Cr2 O3 that has been widely studied and in which the exper-
imental reports show a ME coupling in the range between �.�� – �.�� ps·m-1
[���, ���, ���, ���], the linear magnetoelectric response in NaMnF3 is in the
range, when the strain is close to +�%. Additionally, the ME response in this
fluoroperovskite is of the same order as in CaMnO3 [���].
Additionally, it can be noted that all the ME coefficients increase when the
strain goes from positive to negative values (for example the ↵yy coefficient
varies from -�.��� to -�.��� ps·m-1 for +�% to -�% respectively), which is
consistent with the fact that the polarization contribution of the magnetic Mn
atoms is larger for compressive strains and thus we can expect a stronger cou-
pling between electric polarization and magnetism. The global form of the ME
responses in this fluoroperovskite is very similar to the one predicted for CaMnO3
[���], which shows that this might be a general rule for strain-induced FE in the
Pnma structure. Finally, we also note that the ME response is observed in spite
�� induced multiferroic and magnetoelectric behavior in n a m n f 3
Table �.�: ME coefficients of strained NaMnF3 Pna�1 expressed as in the Eq. �.�, ↵ik
in [ps·m-1 ] and jik in [⇥��-8 ps·A-1 ].
Strain, " [%] ↵yy ↵xx yyy xxx yxx
+�% -�.��� �.��� -�.��� �.��� -�.���
�% -�.��� �.��� -�.��� �.��� -�.���
–�% -�.��� �.��� -�.��� �.��� -�.���
–�% -�.��� �.��� -�.��� �.��� -�.���
of a non Mn-driven polar distortion which claims for a study of the ME-coupling
source as discussed next.
�.� non-collinear magnetism and me coupling in pna� 1 NaMnF 3 ��
A phonon projection was carried out, then, the effect of the magnetic field and
the related modes activated as a function of it were analyzed. In Fig. �.� we plot
the overlap between the magnetic B-field induced distortions in the Pna�1 phase
(at �� T amplitude of the field along the y-axis, here denoted as | i in Fig. �.�)
and the Pnma phonon eigendisplacements basis set (|⌘i in Fig. �.�) for epitaxial
strain values of -�%, �%, and +�%. The eigendisplacements are normalized as
follows h⌘i |M|⌘i i = 1 where M is the mass matrix of the atoms.
Figure �.�: Projection of the magnetic B-field induced distortions (at �� T in the Pna�1
FE phase) against the phonon basis set of the Pnma phase for -�%, �%, and
+�% values of the epitaxial strain. We note that the B 2u FE unstable mode
(the one with negative frequencies) is strongly dominant at compressive
strain while the A g modes are more dominant at tensile strains. This fulfills
the fact that at compressive strain the FE distortions are much more Mn-
driven (and thus much more sensitive to magnetic field) than at tensile
strains and this also explains why under compressive strain the ME response
is larger than under tensile strains.
Our results suggest that the mechanism driving the magnetoelectric coupling in
the NaMnF3 is related to the activation of the Ag modes (Na+F sites) induced
by applied magnetic field (see Fig. �.�). Thus, the activation of those modes
modify the eigendisplacements of the polar B2u phonon and then, tunes the
FE polarization. This type of ME effect in which, a field-induced distortions
are dominated by different phonon modes others than the responsible for the
�� induced multiferroic and magnetoelectric behavior in n a m n f 3
In order to analyze the effect of the Coulomb repulsion and exchange interaction
on the magnetoelectric coupling, calculations for different U and J constants were
performed. In Table �.� presentes the computed constants for the ↵ik and jik
terms. It can be observed that the linear ME terms decreases by approximately �.�
ps·m-1 for ↵xx and ↵yy . However, the second-order coupling term jik increases
by one order of magnitude. Additionally, it can be appreciated that all the tensor
components decreases when the J term increase from �.� to �.� eV. In conclusion,
even when the values for the magnetoelectric coupling changes, the deviation is
rather small and the effect in the ME coupling is not very large when LDA+U
approach is used. Besides, the presence of the ME-coupling is strongly supported
by the symmetry analysis.
Table �.�: ME coefficients (↵ik and jik in ps·m-1 and ⇥��-8 ps·A-1 respectively) of
strained NaMnF3 as a function of U and J parameters [in eV] for the �%
value.
U, J [eV] ↵yy ↵xx yyy yxx
�, �.� -�.��� �.��� -�.��� -�.���
�, �.� -�.��� �.��� -��.��� -��.���
�, �.� -�.��� �.��� -��.��� -��.���
�, �.� -�.��� �.��� -��.��� -��.���
The B term in the Eq. �.� reflects the non-linear character of the strain-polarization
coupling of a FE crystal. In order to check whether this nonlinearity is really spe-
cific to materials presenting a steric geometric origin of the FE, we have analyzed
the polarization-strain coupling in different cases. We first looked at this response
in BaCoF4 (see Fig. �.�), which is known to be a geometric FE without anomalous
Z⇤ . Interestingly, we found the same parabolic polarization vs. strain behavior as
for NaMnF3 where the B-term is approximately ��% of the linear one.
We also report in the figure the polarization vs. strain (P vs. ") of cubic BaTiO3 ,
cubic CaTiO3 and cubic NaMnF3 (i.e. in the absence of the AFD motions in the
last two of them). We can see that in BaTiO3 the P vs. " coupling is linear while
�.� non-collinear magnetism and me coupling in pna� 1 NaMnF 3 ��
Figure �.�: Polarization (out-of-plane) vs. epitaxial strain in the ab-plane for the geo-
metric FE BaCoF4 , for cubic BaTiO3 , cubic NaMnF3 and cubic CaTiO3
case in which we simply relax the polar distortion without the AFD motions,
and for improper-FE PbTiO3 /SrTiO3 and CaTiO3 /SrTiO3 superlattices.
We note that the strong non-linear polarization/strain coupling is present
mostly in steric geometric FE crystals such as NaMnF3 and BaCoF4 .
it is strongly non-linear for NaMnF3 (which means that the octahedral motion
do not play a determinant role for this mechanism). We note that for CaTiO3 ,
the P vs. " is linear but it also exhibits some non-linearity. We would thus con-
clude that BaTiO3 is driven by a charge transfer mechanism while in CaTiO3 a
mixed charge transfer and steric geometric effect are present. To go beyond, we
also performed the same analysis in improper-FE superlattices PbTiO3 /SrTiO3
(PT/ST) and CaTiO3 /SrTiO3 (CT/ST) and we found that in PT/ST the P vs.
" is linear and slightly non-linear in CT/ST. Improper-FE is thus not necessar-
ily related to a geometric effect in the sense of our steric definition and thus
NaMnF3 has a unique response in the perovskite family. This reinforces the idea
that the linear vs. non-linear polarization-strain coupling can be a key ingredient
to distinguish between a charge transfer FE and a steric geometric FE but further
studies would be welcome to confirm this hypothesis. Nonetheless, it seams that
this large non-linear piezoelectric response is very unique to the steric geometric
origin of the ferroelectricity.
�� induced multiferroic and magnetoelectric behavior in n a m n f 3
�.� Conclusions
Using density functional theory we studied the effect of the epitaxial strain on the
structural, electronic, and magnetic properties in the fluoride perovskite NaMnF3 .
A multiferroic-induced state was predicted for NaMnF3 under epitaxial strain at
compressive or tensile strain. A non-linear behavior of ferroelectric instability,
and then, a non-linear piezoelectric response as a function of epitaxial strain was
found. The later completely different as it has been reported in perovskite oxides.
Additionally, an out-of-plane polarization, with respect to the epitaxial strain,
was observed again in opposite trend when comparing to oxides. A cooperative
FE ordering between Na and Mn sites was found for compressive strain. The
later against a pure A-site geometrically driven ferroelectricity at tensile values of
orthorhombic ac-strain. Magnetic ordering reveals a non-collinear ground state
under the Gz Ax Fy representation. And even more interesting, and non-linear
magnetoelectric coupling was found in this material under the strained Pna�1
ground state. We hope that our results can motivate to experimentalist in the
search of this material which is, to our knowledge, the first perovskite-fluoride in
which non-linear piezoelectric and magnetoelectric response co-exist.
Part IV
In addition to the studies in bulk fluorides, much less is known about the mag-
netoelectric properties of interfaces and surfaces of complex fluorides. There are
a few reported examples of complex fluorides thin films and the corresponding
interface characterization. These include Fe/KFeF3 bilayers [���] and the trilay-
ered system Fe/KNiF3 /FeF2 [���] that has a dynamic and rotatable exchange
anisotropy. These systems exhibit exchange bias phenomena, which results in a
center shift of the hysteresis loop of the ferromagnetic Fe layer away from H =
0, as a result of the magnetic interface interaction with the antiferromagnetic
material.
Recently, several efforts have been dedicated to study and to understand mixed
systems by considering the effect of fluorine replacement by oxygen ions in ox-
��
�� �deg, �dhg, fe, and magnetism in oxyfluorides abo 3 /abf 3 interfaces
Introduction: Electronic phenomena at interfaces
Charge transfer
KNiF3
Strain
KTaO3 Charge
Ch
harge
–
Orb
Orbital
Frustration +
Symmetry
Figure 2
The rich behavior of transition metal oxides arises from the complex interactions between their charge,
orbital, spin, and lattice degrees of freedom. These interactions can be modified at interfaces between
different oxides through the effects of local symmetry breaking, charge transfer, electrostatic coupling,
strain, and frustration, leading to fascinating new phenomena.
Figure �.�: Some of the physical phenomena that have lead to several multifunctional LaAlO3/SrTiO3 system. The close proximity in energy of different structural phases and
ide/fluoride (KTaO3 and KNiF3 in the figure) interfaces. The diagram prop- popular method for controlling the lattice distortions, which in turn influence the elec-
tronic degrees of freedom, providing a powerful tool for adjusting materials’ properties
and stabilizing alternative ground states. The advances in the field of strain engineering of
erties has been taken from Zubko et al. [���]. oxide heterostructures are outside the scope of this review, and the interested reader is
referred to References 34–37.
When polar materials are involved, electrostatic interactions also come into play, and the
P. Zubko et. al, Annu.
boundary Rev. Condens.
conditions Matter
at the interfaces lead toPhys. 2, 141–165
an electrostatic coupling(2011).
between the layers,
which is discussed in some detail in Section 4.5. In addition, competition between different
interactions at interfaces and frustration, for example of spins, can also result in qualitatively
Previous reports in oxyfluorides: Moon, E. J. et al.. J. Am. Chem. Soc. 136, 2224–2227 new(2014).
structural and electronic phases. These concepts of symmetry breaking, charge transfer,
tion” has been successfully applied to epitaxial growth of SrFeO3-x Fx (x from An important issue to consider is that of length scales. Although a particular phenomenon
may have its origin at the interface, it is not necessarily confined to it. Unscreened Coulomb
0 to 1) thin films [���] where substantial changes in the electric conductivity www.annualreviews.org ! Complex Oxide Interfaces 147
Table �.�: Electronic configuration of pseudopotentials used in the analysis of this chap-
ter. Valences and semicore electrons taken into account for calculations.
A-site B-site X-site
K (�p6 �s1 ) Ta (�p6 �d5 �s2 ) O (�s2 �p4 )
Na (�p6 �s1 ) Zn (�d10 �s2 ) F (�s2 �p5 )
La (�s2 �p6 �d1 �s2 ) Ni (�p6 �d8 �s2 )
Mn (�p6 �d5 �s2 )
Sc (�p6 �d1 �s2 )
CRYSTAL��— Due to the fact that some of the interesting interfaces count
with a large number of atoms and in order to describe correctly the strongly cor-
related nature of the electrons, hybrid calculations have been performed using the
linear combination of atomic orbitals (LCAO) as implemented in the CRYSTAL
code [��]. This code uses a Gaussian-type basis set to expand the Kohn-Sham or-
bitals, including polarization effects, with the advantage that less number of the
basis is required, which impact into the requirements of the computational time
used to calculate strongly correlated materials when compared to plane wave ba-
sis set as it is the case of VASP. Additionally, all the electrons have been included
in K, Ni, Zn, O, and F while a pseudopotential method has been used for the
case of Ta. For the exchange-correlations term, the hybrid functionals taken into
account are the B�WC [��], B�LYP [���, ���], and PBEsol0 (derived from the
PBE0 functional but based in the PBEsol approximation) [��] all of them with a
�� �deg, �dhg, fe, and magnetism in oxyfluorides abo 3 /abf 3 interfaces
��%, ��%, and ��% of exact exchange-correlation term respectively (see Chapter
�).
KTaO2.5F0.5 KNiF2.5O0.5
Figure �.�: Localized-DOS for the �d-Ni and �d-Ta at the KNiF2.5 O0.5 and
KTaO2.5 F0.5 oxyfluorides respectively. This DOS was computed by means
of the hybrid functional HSE06.
in SrFeO3-x Fx thin films [���]. Then, the tetragonal distortion and the change
in the electronic occupation of the metal cations can induce a tetragonal (T )
splitting of the d-orbitals.
In order to understand in more detail, we computed the localized DOS by means
of the hybrid functional HSE06 (Heyd-Scuseria-Ernzerhof) [��]. This functional
is known to be more precise in the electronic structure calculations correcting the
problems due to the a lack in the description of the exchange-correlation energy
(see Fig. �.�). In KTaO2.5 F0.5 the lowest orbitals are the dxz/yz followed by dxy .
These orbitals exhibit a partial occupation that is in agreement with the expected
Ta oxidation state as discussed previously. The electronic configuration of the d-
orbitals for the Ta are depicted in Fig �.�. In the KNiF2.5 O0.5 oxyfluoride we
can see that the lowest occupation belongs to the dxy orbital, followed by the
dxz/yz orbitals. It is important to note that the dxz/yz is slightly occupied,
which is induced by the d(8- ) Ni state. Next, we can observe a half occupied
dx2 -y2 orbital and right after, around -� eV, the dz2 orbital. This electronic
configuration corresponds to the T -compression configuration as depicted in Fig.
�.�. A magnetic moment is observed in the Ni case and the ordering remains
G-type AFM as expected from the former KNiF3 . Nonetheless, in the Ta-based
case, we do not observe a magnetic moment for the partial �d occupation.
Electronic Density of States
c = 4.127 Å
1.5
1.0 Electronic Density of States
dyz , dxz
−3.0
4.5
KTaO3 KTaO2.5F0.5
−6.0
c = 4.127 Å
1.5
−6.5
1.0
0.5
−10.0 −8.0 −6.0 −4.0 −2.0 0.0
0.0 ElectronicEnergy
Density
(eV)
of States
−0.5
d(2+d)+
Density (states/eV)
5d:t2g
−2.0
75.0
KTaO3 KTaO2.5dF
−6.0
30.0
20.52
x -y (b1g)
−6.5
a = b = 4.005 Å
25.0
dx2-y2 dz2
20.0 −10.0 −8.0 −6.0 −4.0 −2.0 0.0
15.0 ElectronicEnergy
Density
(eV)
of States
10.0
ddz2xy(a(b1g2g) ) −35.0
50.0
−40.0
45.0
(Ni-O)z: 1.813 Å
−45.0
dx2-y2 (b1g)
−60.0
a = b = 4.005 Å
25.0
−65.0
dx2-y2 dz2
20.0
−70.0
15.0
KNiF3 KNiF2.5O0.5
−75.0
10.0
c = 3.837 Å 5.0 −7.0 −6.5 −6.0 −5.5 −5.0 −4.5 −4.0 −3.5 −3.0 −2.5 −2.0 −1.5 −1.0 −0.5 0.0 0.
0.0 Energy (eV)
−5.0
z
−60.0
−65.0
−70.0
KNiF3 KNiF2.5O0.5
−75.0
−7.0 −6.5 −6.0 −5.5 −5.0 −4.5 −4.0 −3.5 −3.0 −2.5 −2.0 −1.5 −1.0 −0.5 0.0 0
Energy (eV)
x
y
��� �deg, �dhg, fe, and magnetism in oxyfluorides abo 3 /abf 3 interfaces
(KZnF3)m/(KTaO3)n
Insulator
m/n: 1/1
Metal m/n
Ta+5O2-2/K+1F-1 Zn+2F2-1/K+1O-2
m/n: 7/7 n-type p-type
B
Figure �.�: �/� and �/� superlattices at which the insulator-to-metal transitions takes
place. The TaO2 /KF and ZnF2 /KO interfaces can be observed.
KZnF3
ZnF2 /KO
TaO2/KF
ZnF2/KO
TaO2 /KF
KTaO3
KTaO3
KZnF3
Figure �.�: Atomic displacements of the B-sites along the z-axis (left) with respect to
the oxygen apical plane for the �/� superlattice. K–K c-lattice parameter in
which, a large expansion and compression can be observed at the TaO2 /KF
and ZnF2 /KO interfaces respectively (right).
(KZnF3)1/(KTaO3)1 (KNiF3)1/(KTaO3)1
Figure �.�: Localized density of states for M = Zn and Ni in the �/� superlattice ratio
obtained considering the HSE06 hybrid functional.
obtained based in the hybrid functional calculation. We found that the correct
electronic structure and magnetic moment is expected for U values above �.� eV.
Besides, without the use of an U value within the GGA calculations, the unoccu-
pied Ni states are located between the valence oxygen states and the conduction
Ta bands, which represents a wrong description of the electronic structure when
is compared to the one obtained from the use of the hybrid functional.
Within the hybrid functional calculations, we also found a weak magnetic mo-
ment of m = �.��� µB /atom localized on the Ta cations. These moments are
following a FM interaction along z with respect to the Ni sites ordered antifer-
romagnetically with m = �.��� µB /atom in agreement with a Ni+2 . This is ex-
plained by the strongest FM interaction between half-occupied �d–Ni eg and the
empty �d–Ta [���] rather that an AFM interaction driven by the super-exchange
interaction between the Ni–O–Ta or Ta–F–Ni.
The dynamical and vibrational stability of the �/� superlattice was investi-
gated by computing the phonon-dispersion curves. We found no unstable mode
in the full Brillouin zone, which suggests a stable structure of the superlattice
configuration as shown in Fig �.�.
The basic superlattice structure is shown in Fig. �.� for B = Zn where alternating
n- and p-type interfaces can be appreciated. The insulator-to-metal transition
(IMT) is a function of the number of layers taken into account in the calcula-
tions. This is in full agreement with the “polar catastrophe” scenario expected
from unbalanced charge at the Zn+2 F-1 2 /K
+1 O-2 and Ta+5 O-2 /K+1 F-1 inter-
2
faces [���]. Then, after several number of layers we found the presence of two-
dimensional electron gases (�DEG) and a two-dimensional holes gases (�DHG)
in KTaO3 /KBF3 superlattices induced by the polar catastrophe and the local
electronic density reconstruction which, is discussed in detail below.
First, we found a decreasing band gap with an increasing of the number of layers
in each compound with a stoichiometry (KTaO3 )n /(KBF3 )l , where n and l are
the number of atomic layers in each compound. This occurred mainly due to the
electrons and holes donated by each layer at the interface and the compensation
of the electrostatic potential generated by the additional electrons in the system
due to the O–B–F bondings. The insulator-to-metal transition was observed for
the n = l = � in both cases. It is important to note that the transition is achieved
with less polar/non-polar layers than in oxides systems such as LaAlO3 /SrTiO3
in which, the IMT is observed for values above ��/�� ratio [���]. Another point
that is important to remark is the localization of the �DEG and �DHG (Fig. �.�)
that, according to our results comes from the dxy –Ta and O p-orbital respectively.
Thus, the fluoride film is not actively having an effect close to the Fermi energy in
the electronic structure and the electrons and holes are only present at the oxide
film. Then, the fluoroperovkite film is acting as a “condensating” layer retaining
the �DEG and �DHG at the oxide layers. Additionally, from the electrostatic point
�.� KT aO 3 / KBF 3 ( B = Zn and Ni ) superlattices ���
(KZnF3)7/(KTaO3)7 (KNiF3)7/(KTaO3)7
Zn+2F2-1/K+1O-2 Ni+2F2-1/K+1O-2
p-type p-type
KZnF3 KNiF3
G-type AFM
TaO2-2/K+1F-1 Ta+5O2-2/K+1F-1
n-type n-type
KTaO3 KTaO3
~ 2eV ~ 2eV
Figure �.�: Layered revolved DOS for the KTaO3 /KZnF3 and KTaO3 /KNiF3 (�/�)
superlattices where a �DEG and �DHG are observed at the TaO2 /KF and
MF2 /KO.
of view a larger electric field across the entire system was found in both cases
in comparison to the LAO/STO oxide/oxide counterpart. For example, when B
= Zn: E = ��� mV/Šagainst E = �� mV/Š[���]. At the metallic interface we
observed that the effective masses in both cases (B = Zn and Ni) are close to
�.�me for the electrons at the �DEG that lies in the lowest �dxy :Ta band. In
contrast, heavier hole carriers with m⇤ of �.�me at the �DHG localized at the
p-orbital of the oxygens at the KO layer. Then, in the �DEG we have a high
mobile carriers at the interfaces, lighter that in the SrTiO3 /LaAlO3 [���]. In
order to explain the electrostatic behavior and the formation of the �DEG and
�DHG we presented the Band-edge diagram built for this system in Fig. �.�. We
can see that the conduction fluoride band-edge lies well above the Fermi level,
besides, the valence band is almost ��� meV below the Fermi energy. Then, it is
again observed that the �DEG and the �DHG are localized at the oxide layer.
Now, we study the magnetic ordering as another degree of freedom in the oxyflu-
orides system for the case of B = Ni. Our results suggest a non-collinear mag-
netic ordering at NiF2 /KO interface, which resulted in a large net magnetization
around �.� µB /atom near the interface for the �/� superlattices at low U values.
Nonetheless, for U values above the fitted one (U > � eV) the magnetic ordering
converged to a pure G-type AFM.
The O–B–F bonding is one of the main points and differences with respect
to other oxides-based superlattices that completely introduces a new degree of
freedom to the oxyfluoride interfaces never explored before. This bonding induces
strong differences in the electronic structure and orbital distribution at the in-
��� �deg, �dhg, fe, and magnetism in oxyfluorides abo 3 /abf 3 interfaces
CB
E VB
2DHG 2DEG
EF
KTaO3 KZnF3
KTaO3/KZnF3 U=0eV
Figure �.�: Band-edge
IMPORTANT:diagram for the electrostatic behavior of the system where the
electric field profile and the �DEG (�DHG) formation are presented at the
Even with an U = 6eV in the KTaO inducing an increase of 1eV in the gap, the system is still metal.
3
Then, the polar catastrophe is still taking place and the electric field goes from 105 to 140 meV/Å. (in
oxyfluoride
comparison to TaO
57 meV/Å /KF
2 in LaAlO(ZnF
3 /SrTiO ). /KO) interface. This was obtained based on the
3
2
layered resolved
Can we say DOS
that the fluorides (see
is still Fig
insulator but �.�).
the oxideCB, VBmaterial?
is the metal an EF denotes the conduction-
band, valence-band and the Fermi energy respectively.
terfaces as it can be observed in Fig. �.�. In this figure we can see the breaking
of the degeneracy of the d-orbitals due to the tetragonal (T ) like distortions. Be-
sides, the d-orbital configurations for both types of interfaces can be the source of
the particular magnetic behavior. Due to the strong electronic correlation in the
partially occupied �d Ta orbitals and the T Ni states at the interfaces as shown
previously, we can see that the magnetic ordering at the interface is changing
with respect to the U value. We can observe, for the �/� superlattice, that the
magnetic state goes from a pure G-type AFM at U =� eV to a C-type AFM plus a
FM layer at the interface for U = �.� eV. To get a closer look in the disagreement,
calculations based on hybrid functionals were carried out.
In order to corroborate the results related to the magnetic ordering and because
of the limitations in the computational efficiency of the calculation with hybrid
functionals within the plane-wave expansion, we used the CRYSTAL code. Using
this code, we selected three different approaches within the hybrid approximation:
B�WC, B�LYP, and PBEsol0. Several magnetic orderings were tested, G-, C-, A-
AFM respectively and mixed G-AFM+FM layer (at the NiF2 immediately after
the TaO2 /KF interface and at the NiF2 /KO layer). Nevertheless, according to
the calculations the magnetic ordering always converges to the G-type AFM as
ground state in the KNiF3 film. These results are in full agreement with those
obtained with the GGA+U approximation for U values above � eV. However,
it is important to mention that the mixed G-AFM+FM-layer (at the NiF2 /KO)
state is very close in energy, E = �.�� meV/atom, closer than other magnetic
ground state (i.e. A- and C-AFM).
Interestingly, there is a magnetic-moment modulation across de entire thickness
of the magnetic counterpart in all the superlattices taken into account. This
modulation can be observed in Fig �.�� where variations between a range of �.��
µB /atom are expected.
Even if a charge-transfer, charge-ordering, and orbital ordering is not achieved
in these interfaces due to the small valence deviations at the interfaces (as dis-
�.� KT aO 3 / KBF 3 ( B = Zn and Ni ) superlattices ���
Figure �.��: Magnetic moment per layer at the KTaO3 /KNiF3 superlattices. A moment
modulation can be appreciated when n/l ratio is varied from �/� to �/�.
cussed in the oxyfluoride section) we have to keep in mind that fluorides can be
better candidates for a charge-ordering state where a larger charge-disproportion
cation is observed [���] in comparison to their oxide counterparts.
With the main idea to understand the isolated p- and n-type interfaces, we have
studied both terminations. A KO/ZnF2 and TaO2 /KF terminated computational
cells were used respectively where a KO plane is set at the center as a mirror
symmetry plane of the full heterostructure. Then a �� u.c. was used, �–KTaO3
and �–KZnF3 respectively. After full relaxation we have found magnetism and
interestingly, this magnetism is coming from two different sources. In the p-type
interface, p-orbital magnetism has been found in agreement with the electron
vacancies and the �DHG formed at the interface (see Fig �.��). The magnetic
structure shows a AFM-type ordering with mT a = �.���µB /atom and mO =
�.���µB /atom. On the other hand, in the n-type interface a �DEG is clearly
formed as shown in Fig �.��. This in addition to the d-orbital magnetism coming
from the Ta:dxy (see inset in Fig �.��) orbital in a A-type AFM with a magnetic
moment mT a of �.���µB /atom at the surface and a �.���µB /atom at the layer
right below. From these interfaces it can be noted as well that the effective electron
masses for hole and electron carriers are �.�me and �.�me respectively. Thus, the
n-carriers are lighter and with a high mobility perfect for electronic applications
[���]. It is important to note that the fluoride layer (i.e. in this case KZnF3 )
is not actively participating in the phenomena, however, is acting as insulating-
confining layer.
We presented in Fig. �.�� the bands for the case of the magnetic fluoride layer.
We can observe that the carriers in the n-type interface are slightly heavier in
comparison to the non-magnetic case (i.e. m⇤ = �.�me against m⇤ = �.�me ).
For the p-type case, the hole carrier present almost the same behavior with m⇤ =
�.�me . Additionally and surprisingly, in the n-type interface we observed a charge-
ordering at the third lower layer with respect to the surface. This is concluded
��� �deg, �dhg, fe, and magnetism in oxyfluorides abo 3 /abf 3 interfaces
Figure �.��: Bands structure for the p- (left) and n-type (right) interfaces along the
X–M– –X path in the Brillouin zone for the KTaO3 /KZnF3 system. The
formation of the �DHG and �DEG at the p- and n- interfaces can be
appreciated.
from a change in the magnetic moment from -�.��� µB /atom to �.��� µB /atom
suggesting an d8 – d6 occupation. This could explain the deeper penetration of
the �DEG below the Fermi level, indicating a larger carrier density in the B =
Ni n-type interface than in the B = Zn one.
In the p-type interface Fig. �.��, we do not observe a charge-ordering. We
found that the magnetism at the KNiF3 remains in its G-type AFM ordering in
addition to the p-orbital magnetism. The later is also observed in the oxyfluoride
interface explained previously with a non-magnetically active B-site.
In comparison, even when the carriers in the n-interface are heavier for B = Ni,
it is important to mention that the �DEG is more pronounced, which in this case
lies at -�.� eV against -�.�� eV for B = Zn. Knowing that the in-plane lattices
parameters and polarization across the interface are quite the same, it can be
suggested that the charge missing at the KNiF3 layer is transferred directly to
the �DEG increasing the amount of density carriers at the �DEG in the interface.
However, further studies and analysis are needed to explain and exploit this
particular behavior in the KTaO3 /KNiF3 n-interface oxyfluoride.
Figure �.��: Bands structure for the p- (left) and n-type (right) interfaces along the
X–M– –X path in the Brillouin zone for the KTaO3 /KNiF3 system. The
�DHG and �DEG at the p- and n- interfaces can be appreciated. Besides,
and empty Ni:eg band is observed between �.� to �.� eV respectively.
showing a similar behavior in SrTiO3 - and KTaO3 -based transistors [���, ���].
Then, for these single perovskite-based devices, a k-cubic dependence of the split-
ting was found (see Eq �.�) in contrast to the linear plus cubic dependence of the
Rashba splitting at the LaAlO3 /SrTiO3 interface due to the dxy – dxz/yz multior-
bital nature of the lowest bands [���]. Additionally, in the case of Au(���) metal
surface, it has been shown that even when the surface states exhibit a mainly
p–orbital character, the d-orbitals drives the Rashba splitting and dictates the
direction of the spin by means of the orbital-angular momentum (OAM) [���, ���].
Therefore, the transition metal d-orbitals at perovskite structures can be the key
parameter for a large and controllable electron band splitting in future spintronic
and spin-orbitronic devices [���, ���]. This later property gives an additional de-
gree of freedom to the already rich field of multifunctional perovskites. Then, in
order to look for the spin-texture and spin-orbit coupling effects in the oxyfluoride
interfaces, the analysis of the electronic structure is done by taking into account
the SOC. Due to the large computational cost, the study of the bands was per-
formed in the �/� superlattice. Notice that these results can be traced up to the
�/� system where, the �DEG and �DHG is formed and the IMT takes place as
shown previously. Thus, the results presented here can be extrapolated to larger
n/l superlattices. in Fig �.�� it can be observed the band structure along the
kx path in which, the sx , sy , and sz spin components are presented separately
and a colors notation of the up (red) and down (blue) spin was adopted. First, it
can be noted that the electronic structure and the spin-texture at the conduction
bands above to Fermi level (and thus, close to the TaO2 /KF interface at the �/�)
are entirely �D spin-polarized at the kx –ky plane. Therefore, no band-dispersion
is observed for the sz spin component. Surprisingly, a large spin-splitting of the
��� �deg, �dhg, fe, and magnetism in oxyfluorides abo 3 /abf 3 interfaces
bands is appreciated for an energy around E = �.� eV. Additionally, the SOC
effect is observed at the heavy dxz/yz crossing with the dxy lighter bands (see
the upper inset in Fig �.��).
Δk = 0.08 π/a
ΔE = 61 meV
sx sy sz
2π/a 2π/a 2π/a
Figure �.��: Band structure for the k x path in which, s x , s y , and s z spin components
are taken into account. The up and down orientation of the spin are repre-
sented by red and blue colors respectively.
where ↵3 is the cubic-Rashba coupling term which measures the strength of the
splitting (in units of eV·Å3 ), ± = �/�( x ± i y ) and x,y are the Pauli matrices.
The Ez is the electric field perpendicular to the electron plane movement denoted
by k± = (kx ± iky ) [���]. Here, the E vs. k profile for the bands is represented
by the relationship:
h2 k2
E± (k) = ± ↵3 k3 . (�.�)
2m⇤
This type of spin-splitting can be understood like an opposite rotation of the
lowest spin–up band with respect of the spin-down around an axis fixed at the
vertices of the parabolic d-bands. This is in contrast with the linear horizontal dis-
placement (or offset) of the spin–up and spin–down bands respectively observed
in the most common linear-Rashba effect� present in materials such as BiTeI
[���].
� The linear-Rashba effect has the E vs. k form: E± (k) = (h2 k2 /2m⇤ ) ± ↵|k| coming from the
Hamiltonian of the form HR = ↵Ez i(k- + - k+ - )
�.� cubic-rashba spin-splitting at oxyfluoride interface ���
Table �.�: Layer-by-layer spontaneous polarization along the z-axis (Pz in µC·cm-2 ) in
the �/� superlattice is shown. Besides, the oxide and fluoride polarizations
are PKT aO3 = -��.�� µC·cm-2 and PKZnF3 = �.�� µC·cm-2 respectively.
KTaO3 KZnF3
Layer � � � � � � � �
Pz -�.�� -�.�� -�.�� -�.�� �.�� �.�� �.�� �.��
From the magnetic superlattice, KTaO3 /KNiF3 we noticed that the same spin-
texture remains at the TaO2 /KO. Nonetheless, the spin-texture is becoming more
complex due to the overlapping of the �d:Ta band with the unoccupied �d(eg ):Ni
states above the Fermi level.
It is important to take advantage of the fact that the magnetic properties
can be conserved in the KNiF3 . Then, a mixed multifunctional properties of
��� �deg, �dhg, fe, and magnetism in oxyfluorides abo 3 /abf 3 interfaces
LaScO3 NaMnF3
-2.6 μB/Mn
-4.6 μB/Mn
4.6 μB/Mn
-4.6 μB/Mn
4.6 μB/Mn
4.6 μB/Mn
4.6 μB/Mn
-4.6 μB/Mn
4.6 μB/Mn
-4.6 μB/Mn
2.6 μB/Mn
-4.6 μB/Mn
4.6 μB/Mn
-4.6 μB/Mn
A. H. Romero, [email protected] PACIFICHEM, December 15-20, Hawaii
Figure �.��: LaScO3 /NaMnF3 (�/�) superlattice. The AFD distortions and rotations
can be observed in the upper figure in which the KO/MnF2 and ScO2 /KF
interfaces can identified. In the lower figure it is shown the spin-polarized
charge density, red and blue colors denote the up and down magnetic mo-
ment directions. The G-type AFM ordering is observed as well as the CO
state.
The results shows a charge-ordering (CO) state found at the middle Mn layer
in the fluoride NaMnF3 part (see Fig. �.��). The CO state is demonstrated by a
change in the magnetic moment from �.� µB /atom to �.� µB /atom as shown in
Fig. �.��. This moment difference suggests an d5 – d3 Mn occupation. Nonethe-
less, for a U = � eV (value normally used in Mn at oxides and fluorides [�, �, �]),
the CO-state remains but now, the magnetic moment at the middle of the layer
goes from �.� µB /atom to �.� µB /atom. The later claims for a deeper study and
understanding of the exchange-correlation in this particular oxyfluoride. We ob-
served that this CO state remains after several U values (from � to � eV). Besides,
the main magnetic ordering exhibit a G-type AFM after testing several possible
states (i.e. G-, C-, and A-type AFM and FM).
The analysis of the spontaneous polarization reveal values of Pz = ��.�� µC·cm-2
and Pxy = �.�� µC·cm-2 . Here, the Pz is induced by the symmetry breaking at
the KO/MnF2 and ScO2 /KF interfaces and the F–Mn–O and O–Sc–F bondings
respectively as can be seen in Fig. �.��. This polarization cannot be reversed,
nonetheless, can be used to enhance the polar catastrophe phenomena and tune
the �DEG as demonstrated in SrTiO3 /PbTiO3 [���], in this interface the polar
catastrophe is neglected due to the neutral character of the atomic layers. How-
��� �deg, �dhg, fe, and magnetism in oxyfluorides abo 3 /abf 3 interfaces
ever, the ferroelectric polarization along the PTO layer induce a �DEG. It is
important to note that possibly this polarization can be coupled with the strain-
induced Pz in the NaMnF3 as shown in the Chapter �. The Pxy component is
explained in terms of the hybrid-improper FE phenomena where different A-sites
at the orthorhombic interfaces brakes the anti-polar displacements and a net po-
larization can be obtained in the xy-plane [���]. This type of displacements can
be observed in the Table �.� in which are presented the atomic displacements per
layer at the xy-plane
Table �.�: A -site displacements per layer in the xy -plane taking the high-symmetry
position in the cubic symmetry as reference. Displacements in Å units.
LaScO3
-�.�� +�.�� -�.�� +�.�� -�.�� +�.�� -�.��
NaMnF3
+�.�� -�.�� +�.�� -�.�� +�.�� -�.�� +�.��
This preliminary results are very promising but further studies are needed in
order to unveil the geometric+hybrid-improper ferroelectricity in these superlat-
tices. Besides, the possible coupling of the FE with the CO-state and the magnetic
degrees of freedom needs to be investigated.
The n/l ratio effect for several n and l values and the study of the �DEG above
the �/� superlattice ratio is on their way.
�.� Conclusions
We have analyzed by means of first-principles calculations the electronic and struc-
tural properties of oxyfluorides (KTaO3 )n /(KBF3 )l B = Zn and Ni interfaces. We
found that the orbital levels symmetry at the interfaces is strongly modified as a
function of the O–B–F coordination. The polar catastrophe phenomena also takes
place in the oxyfluoride interfaces, however, we found that less number of layers
are needed in order to achieve the insulator-to-metal transition. We observed that
the magnetism in the KTaO3 /KNiF3 exhibits a moment magnitude modulation.
Nevertheless, the magnetic structure keeps the G-type AFM such as in the bulk
former compound. Surprisingly, we observed a large k3 -Rashba type splitting in
at the oxyfluoride interfaces, at least four times larger than the one reported in
SrTiO3 /LaAlO3 interface and twice of the KTaO3 -based transistor.
Looking for novel oxyfluoride multiferroics, the LaScO3 /NaMnF3 was studied.
This interface show promising preliminary results where charge-ordering, G-type
AFM, hybrid-improper ferroelectricity and possible tuning of the electronic band
�.� conclusions ���
���
��� conclusions and perspectives
We have also analyzed from first-principles the electronic and structural proper-
ties of oxyfluorides (KTaO3 )n /(KBF3 )l superlattices (B = Zn and Ni). We found
that the orbital splitting at the interfaces is strongly dependent of the O–B–F co-
ordination. Moreover, a polar catastrophe phenomena, similar to that reported at
the LaAlO3 /SrTiO3 interface, also takes place in the oxyfluoride polar interfaces.
We observed that the thickness required to achieve the insulator-to-metal transi-
tion is however smaller than in the oxides. We observed that the magnetism in the
KTaO3 /KNiF3 system exhibits a modulation of the magnetic moment amplitude
across the entire superlattices. Surprisingly, we reported a large k3 -Rashba-type
splitting at the oxyfluoride interfaces, at least four times larger than that reported
in SrTiO3 /LaAlO3 interface and twice that of the KTaO3 -based transistor.
As demonstrated for the multiferroics properties in the NaMnF3 and the elec-
tronic properties in superlattices and interfaces of KTaO3 /KBF3 , in fluoride-
based systems, the tuning of properties can be achieved independently. For ex-
ample, the FE in the NaMnF3 is tuned by the A-site and the magnetism from
the B-site avoiding the d0 -ness rule. Another example of this behavior is the mag-
netic properties that remains in the KNiF3 and the �DEG and �DHG localized at
the KTaO3 in the KTaO3 /KNiF3 system. Then, an additional effort was carried
out in order to tune the magnetic and electric properties. This is the case of the
LaScO3 /NaMnF3 systems in which, MF/ME and �DEG can be expected from
the fluoride and oxide layer respectively.
On the side of what we presented in the main manuscript we also performed side
projects that are presented in the Appendices. In Appendix C, we present the
study of the multiferroic behavior of BaCoF4 thin films under epitaxial strain
where we found a non-linear dependence of the spontaneous polarization as a
function of the external constrain. Besides, several magnetic transitions were ex-
plained in terms of the non-collinear magnetic ordering. In Appendix D we present
a study of the electronic and particularly, spin-texture properties of the SrBO3
(���) with B = Ti and Hf surfaces. We found a tangible Rashba-type splitting
for B = Hf, besides, despite the reported giant-splitting for B = Ti, we did not
found such spin-splitting. Finally, in Appendix E we report the insulator-to-metal
transition control and tuning at the LaAlO3 /SrTiO3 interface by means of OH
and H molecules implantation at the LaAlO3 surface.
in order to study the surface effects and the isolated interfaces that are formed
depending on the chosen substrate. For example future studies can be focused on
SrTiO3 (���)/NaMnF3 (���) and LaMnO3 (���)/NaMnF3 (���) thin films.
Finally, we believe that the discovery of new multiferroic materials could lead
to new possible solutions for the challenges in the field, even though the present
result does not fully achieve the “ideal” target of room temperature magneto-
electrics. Moreover, we hope that the results we have presented in this thesis will
motivate further theoretical and experimental studies in the search of new and
exotic multifunctional properties in fluoride and oxyfluoride-based perovskites.
CONFERENCES AND MEETINGS
In the next section are listed the conferences and meeting in which, the thesis
results were presented:
���
Part V
APPENDIX
The next appendix list condenses the work related to scientific col-
laborations that lead to an important learning in my PhD formation
but, that are not directly related with the main topic of the thesis.
Nonetheless, the study keep its focus on the electronic structure of
fluorides and some perovskites.
O V E RV I E W O F T H E M A G N E T I C
A
INTERACTIONS IN
FLUOROPEROVSKITES
Table A.�: Magnetic canted moment Fy (µB /atom) at low temperature in Pnma ABF3
systems. The large magnetic canting for fluorides compounds can be noted
in the Gz Ax Fy configuration. The ionic radius of the magnetic B-site (rB )
was taken into account the valence state and the bonding coordination [��].
Compound Fy Structural parameters
ABX3 U =� eV Exp. [���] d-Occ. rB [pm] a [deg]
���
��� overview of the magnetic interactions in fluoroperovskites
Table A.�: Experimentally reported exchange interaction (J) constants in ABF3 com-
pounds [���]. The components of the DM vector are in meV.
ABX3 J1 J2 Structure
KMnF3 �.��� �.��� Pnma
KFeF3 �.��� — Pm3̄m
KNiF3 �.��� — Pm3̄m
KCuF3 ��.��� �.��� I�/mcm
NaFeF3 �.��� �.��� Pnma
NaMnF3 �.��� �.��� Pnma
LaFeO3 �.�� �.�� Pnma
BiFeO3 �.�� �.�� R�c
occupation, B-site radii, and the tilt angle, . For fluorides it can be noted that
the rotation angle are almost twice of the ones found in the oxides (see Fig. A.�).
Additionally, the weak-FM follows a trend with the octahedral angle: higher is the
angle, larger is the Fy magnetic moment. At the same time, the octahedral rota-
tion increases when the B radii does it is shown in Table A.�. The later magnetic
behavior is in agreement with the Dzyaloshinskii-Moriya interaction [���, ���],
then: Dij / x̂i ⇥ x̂j where x̂i,j are the unitary vector along the M–X–M bondings.
We can observe the reported weak-FM component for the NaFeF3 is rather large,
around �.��� µB /Fe, following the rotation and radii trend. Therefore, the ex-
pected Fy component in the NaMnF3 should be around �.� µB /Fe, which is one
order of magnitude large that the ones found in perovskite oxides. Unfortunately,
the computed values are not in well agreement with the reported experimental
moments. Then, a careful research within the DFT+U +J needs to be carried
out in future studies. The RbFeF3 compound shows the largest Fy component,
nonetheless, its ground state belongs to an hexagonal structure and the structural
differences needs to be taken into account in the analysis.
As expected from the weak covalent bonding in fluorides, comparing with ox-
ides, the exchange-constants, Ji are smaller by at least one order of magnitude
as it is presented in Table A.�. The KCuF3 compound exhibit a quite large J1
component due to the Jahn-Teller distortion presented in the partially occupied
d9 -orbital. Surprisingly, the KNiF3 shows large exchange interaction. These low
exchange interactions can explain the low magnetic transition temperatures ob-
served in fluoride perovskites. Therefore, even when the wake ferromagnetism
tend to be larger than in oxides, the temperature regime for this particular mag-
netic behavior seems to be very low. The later statement needs to be taken into
account in the application of these type on materials in further technological appli-
cations. Nonetheless, compound with several phenomena, such as the Jahn-Teller
overview of the magnetic interactions in fluoroperovskites ���
Figure A.�: Octahedral rotation for several oxides and fluorides compounds. The
fluoroperovskites are ordered according to the B-site radii.
The first step into the analysis is to identify the space group and its point group.
In the case of the post-perovskites we have a Cmcm space group (No. ��) with a
D2h point group, this point group have the symmetry operations: I, -I, 2z , 2y ,2x ,
mz , my , and mx extracted from the Bilbao Crystallographic Server [���, ���].
Then, we have:
h i
I = 1 1 1
h i
-I = -1 -1 -1
h i
2z = -1 -1 1 + ⌧(0, 0, 1/2)
h i
2y = -1 1 -1 + ⌧(0, 0, 1/2)
h i
2x = 1 -1 -1
h i
mz = 1 1 -1 + ⌧(0, 0, 1/2)
h i
my = 1 -1 1 + ⌧(0, 0, 1/2)
h i
mx = -1 1 1
���
��� group theory analysis to the non-collinear magnetism in ppv
Now, we analyze how the magnetic sites transform under the allowed symmetry
operations it can be observed that, for example, the M1 cation transforms to the
M2 under the 2z operation as follows:
2 3 2 3 2 3 2 3
-1 0 0 0 0 0
6 7 6 7 6 7 6 7
2z · M1 = 4 0 -1 05 · 405 + 4 0 5 = 4 0 5 ⌘ M2 (B.�)
0 0 1 0 1/2 1/2
Therefore we have that all the magnetic sites transforms under the symmetry
operations as follows obtaining the Table B.�:
Table B.�: Transformation of the magnetic sites under the symmetry operations. The
numbers in the table are the i index of the Mi magnetic cation.
I 2z 2y 2x -I mz mx mx
M1 � � � � � � � �
M2 � � � � � � � �
M3 � � � � � � � �
M4 � � � � � � � �
+ ⌧(�/�, �/�, �)
M1 � � � � � � � �
M2 � � � � � � � �
M3 � � � � � � � �
M4 � � � � � � � �
Here, at this point, we know how the magnetic sites transform under the sym-
metry operations. We now proceed with the analysis of the magnetic moments
transformation. For this purpose, we have to define the possible magnetic config-
urations into the collinear level. The selected notation is F for ferromagnetic, and
A, C, and G for A-, B- and G-type antiferromagnetic orderings. Then, we have:
In Table B.� the + and - signs are for up and down spin orientation respectively.
We proceed by calculating how the spin transforms under the same symmetry
operations such as the magnetic sites does. Therefore, we have to include the
non-collinear orderings by defining the spin pseudo-vector such as s = (sx , sy ,
sz ). Due to the pseudo-vector character of the spin and in order to understand
how it transforms, we need to be defined such pseudo-vector as the product of A
(where A = (x, y, z)) and B (with B = (x 0 ,y 0 ,z 0 )) as follows: s = A⇥B. Then, s
= (yz 0 - zy 0 , zx 0 - xz 0 , xy 0 - yx 0 ).
group theory analysis to the non-collinear magnetism in ppv ���
Table B.�: Definition of the different magnetic orderings. The G, A-, and C-type anti-
ferromagnetic orderings and the ferromagnetic F state are defined for each
magnetic cation where + and - denoted up and down spin orientations.
M1 M2 M3 M4
F + + + +
A + + – –
C + – + –
G + – – +
To exemplify how proceed, we show how the spin will transform under the 2z
operation as follows:
2 3 2 3 2 3 2 3
-1 0 0 yz 0 - zy 0 -yz 0 + zy 0 -1
6 7 6 0 7 6 7 6 7
2z · s = 4 0 -1 05 · 4 zx - xz 0 5 = 4 -zx 0 + xz 0 5 ⌘ 4-15 · s (B.�)
0 0 1 xy 0 - yx 0 xy 0 - yx 0 +1
Then, after applying to all the operations we have that the spin transforms as
condensed in Table B.�. It is important to note that the s transforms under the
-I as -s = B⇥A. Nonetheless, the spin must be exchanged and undistinguible,
therefore, it transforms in the same way for the -I as for the I symmetry operation.
The later statement is also valid for the mz , my , and mx .
Table B.�: Transformation sx , sy , and sz spin orientations under the symmetry opera-
tions. The transformation of each magnetic ordering under each symmetry
operation is labeled as +� and -�, which indicates when the ordering is invari-
ant under the related transformation (i.e. +�) or when it is reversed under
the application of the symmetry operation (i.e. -�).
I 2z 2y 2x I mz my mx
sx +� –� –� +� +� –� –� +�
sy +� –� +� –� +� –� +� –�
sz +� +� –� –� +� +� –� –�
Table B.�: Transformation of the magnetic ordering under the allowed symmetry oper-
ations.
I 2z 2x 2x -I mz my mx
Fx +� –� –� +� +� –� –� +� B3g
Fy +� –� +� –� +� –� +� –� B2g
Fz +� +� –� –� +� +� –� –� B1g
Ax +� –� –� +� +� –� –� +� B3g
Ay +� –� +� –� +� –� +� –� B2g
Az +� +� –� –� +� +� –� –� B1g
Cx +� +� +� +� +� +� +� +� Ag
Cy +� +� –� –� +� +� –� –� B1g
Cz +� –� +� –� +� –� +� –� B2g
Gx +� +� +� +� +� +� +� +� Ag
Gy +� +� –� –� +� +� –� –� B1g
Gz +� –� +� –� +� –� +� –� B2g
By re-grouping the orderings into the same symmetry character we obtain the
Table B.�.
Table B.�: Magnetic ordering allowed in the Cmcm structure according to the D2h point
group [���, ���].
Magnetic Ordering I 2y 2y 2x -I mz my mx
Fx , Ax +� –� –� +� +� –� –� +� B3g
Fy , Ay , Cz , Gz +� –� +� –� +� –� +� –� B2g
Fz , Az , Cy , Gy +� +� –� –� +� +� –� –� B1g
Cx , Gx , +� +� +� +� +� +� +� +� Ag
To obtain the possible non-collinear states we take into account the superex-
change interaction in order to avoid unprovable orderings. Besides, we neglect the
ferrimagnetic orderings. Then, finally we have that the allowed states are Cz Fy ,
Cz Ay , Cy Fz and Gy Az .
For more details in the group theory analysis applied to the non-collinear mag-
netism the reader can check the Ref. from Bertaut [���].
BaCoF4 THIN FILMS
C
The orthorhombic fluoride family BaMF4 , with M = Mn, Fe, Co, Ni, and Zn
[���, ���, ���] was intensively studied four decades ago [��] due to their interest-
ing multifunctional properties and physics discovered. These fluorides exhibit si-
multaneously antiferromagnetic and ferroelectric long-range orders. While the an-
tiferromagnetic order has a pronounced two-dimensional character, the ferroelec-
tricity mechanisms are fundamentally different from the “classical” charge-transfer
phenomena widely known in ferroelectric perovskites oxides such as BaTiO3 and
PbTiO3 [��, ��]. Recently, new theoretical analysis [��, ���] that includes a pre-
dicted magnetoelectric effects present in some of those materials (for example in
BaNiF4 ) motivated further experimental studies on those multiferroic systems.
As commented before, the BaMF4 compounds crystallize in orthorhombic lattice
with the space group Cmc�1 No. �� [��]. In this structure, a shared-corners dis-
torted fluorine octahedra M–F6 in the ac-plane form puckered sheets, which are
separated by Ba layers along the crystallographic b-axis [��, ���]. A few years ago
was theoretically described the ferroelectric long range order for the family mem-
bers M = Mn, Co, Fe and Ni which has its origin in a combination of octahedral
rotations in M–F6 and polar displacements of Ba cations, in a so-called geomet-
ric proper ferroelectricity [���] making it very appealing to tune multifunctional
properties. Previous experimental reports have revealed ferroelectric switching in
two members of the family, BaCoF4 and BaNiF4 with ferroelectric saturation
polarization along the a-axis and Ps = �.� µC·cm-2 and �.� µC·cm-2 , respec-
tively [���]. However, no switching has been observed for BaFeF4 and BaMnF4 ,
nonetheless, the analysis by Ederer and Spaldin [��] proposed that the ferroelec-
tric switching in those materials is likely possible by means of higher electric
fields. From the magnetic point of view, magnetic moments of M2+ have their
strongest antiferromagnetic interactions in the ac-plane, with nearly linear M–F–
M–F chains along the c-axis, and a zig-zag M–F–M–F configuration along the
a-axis [���, ���]. Since the interplanar spacing along the b-axis is much larger
(i.e. about �� Å) than the intraplanar ones (about �-� Å), the antiferromagnetism
in BaMF4 has a pronounced two-dimensional character. Spins are oriented along
the b-axis [���] except in BaCoF4 , where they lie along the a-axis. It is important
to remark that two antiferromagnetic phases have been found in case of BaCoF4
[���]: phase A and phase B being distinguished by the sign of the interplanar cou-
pling along the b-axis, negative and positive, respectively (see Fig. C.�). Then,
in this appendix we present the calculations and analysis that support the exper-
imental synthesis of multiferroic BaCoF4 fluoride thin films which are for first
���
��� b a c o f 4 thin films
time achieved. For more details in this work please check the Ref. of Borisov et al.
[�].
In the next sections we will discuss the results regarding the effect on the
xz-plane due to an applied epitaxial strain in BaCoF4 thin films. Thus the multi-
ferroic state is detailed analyzed by means of the employed theory. In particular,
the spontaneous polarization as the non-collinear magnetism as a function of the
strain is studied.
ide perovskites and originally paraelectrics, which under epitaxial strain becomes
ferroelectric. The fluoride case is rather different due to its larger sensitivity to
strain in which a paraelectric to ferroelectric transition is achieves for values of
strain less than �.�%. Then, keeping fixed the ac-lattice plane to the experimen-
tal values, internal coordinates and z-axis are allowed to fully relax. After this
structural relaxation a new phase for the strained BaCoF4 was found, going from
the Cmc�1 (symmetry group No. ��) to Pca�1 (symmetry group No. ��). In this
new phase, the ferroelectric behavior remains unchanged due to the polar nature
of the structure as in the Cmc�1 phase. This transition in driven by atomic dis-
tortions (represented as red arrows) according to the 1 mode as shown in Fig.
C.�.
Figure C.�: BaCoF4 [���] supercell in which the red arrows represent the distortion to
go from the Cmc�1 to Pca�1 space groups. Atoms in green, blue, and grey
colors for Ba, Co, and F species respectively. This figure was elaborated
with the support of the VESTA code [���].
As a next step and with the goal to test the thermal stability of the new BaCoF4
phase, phonon dispersion calculations were performed in the relaxed strained
structure with lattice parameters such as found in thin films (approx. +�.�%).
No imaginary modes traduced in unstable modes were found, which suggest a
stable Pca�1 phase structure for this value of strain. After the later analysis, we
proceeded to compute the ferroelectric and magnetic properties in the relaxed
BaCoF4 structures for strains between -�% to +�% to complete the landscape
of the multiferroic properties as a function of the epitaxial strain such as in
experimentally growth thin films.
Fig. C.�, we present the behavior of spontaneous polarization as a function of
the ac-plane strain. It can be observed the non-linear coupling of the polarization
as a function of the strain, which suggest a non-linear piezoelectric effect in this
fluoride compound, which is quite rare in these materials [�]. As expected, the
polarization develops along the crystallographic b-axis such as has been found
in bulk compounds [���]. Nonetheless, small components of polarization along
x and z are observed but are not shown here (with an order of �.� µC/cm2 ).
��� b a c o f 4 thin films
Figure C.�: (Top) Spontaneous polarization along the b-axis. a, b, and c lattice param-
eters (Middle) and behavior of the octahedral angle bonding (Co–F–Co)
along the x-axis (bottom) as a function of the ac epitaxial strain respec-
tively.
This can be due to the additional distortion of the CoF6 octahedra induced
by the epitaxial strain. For values of an epitaxial strain close to those observed
in the BaCoF4 /Al2 O3 (����) thin films, which is close to �.�%, we compute a
polarization value of Ps = �.�� µC·cm-2 . Additionally, Fig. C.� also shows that
the b lattice parameter decrease as a and c increase as can be expected from the
Poisson ratio.
Regarding the magnetic ordering in the collinear approach, three different con-
figurations have been proposed. These configurations are labeled as A, B and C
respectively as depicted in the Fig. C.�,. We found that the B-phase is the most
stable, however, the A and C are at �.�� and �.�� µeV/f.u. larger in energy. This
fact suggests an easy switching of the ordering achieving by small fields and even
temperature with that takes an AFM to AFM magnetic transition.
C.� noncollinear magnetism and multiferroic behavior ���
Figure C.�: From left to right, the A, B and C collinear magnetic orderings. In green,
blue navy and grey are Ba, Co and F ions respectively. Magnetic moments
in each atom are located along the y-axis. In the figures can be observed
the unit cell in black lines as well as the computation magnetic cell that
takes into account �� atoms.
ment with the experimental findings (Gz Cy Fx ). It is important to remark that the
weak-AFM ordering exhibit two different magnetic moment components which
are in agreement with the different octahedral tilting of the CoF6 in the xz-layer
such as can be observed in the Fig. C.� and in the values presented in the Table
C.�.
Table C.�: Magnetic moment per atom (or Co atom) along each axis (x, y and z) for
each magnetic ordering where the Gz is the main collinear AFM order and
Cy and Fx are the weak-AFM and weak-FM spin orders, respectively. All
strains are isotropic in the ac-plane, except for the �.�% strain result which
was calculated with strain applied only along the z-axis, per experimental
results.
Strain mi [ µB /atom]
ac-plain Gz Cy Fx
-�% �.��� �.���,�.��� �.���
�% �.��� �.���,�.��� �.���
+�.�% (exp.) �.��� �.���,�.��� �.���
+�% �.��� �.���,�.��� �.���
c.� Conclusions
We have studied by means of first-principles calculations, the behavior of the
magnetic and ferroelectric properties of the BaCoF4 as a function of the ac-plane
epitaxial strain. We found a phase transition from the Cmc�1 to Pca�1 phase
as a function of the epitaxial strain was found. A non-linear polarization-strain
coupling was observed suggesting a second order dielectric behavior allowed in
the BaCoF4 fluoride. An enhancement of the polarization is observed for values
of strain > �%. For the value of epitaxial strain correlated to the experimental
thin films, a non-collinear ground state was found. Besides, a possible path of
magnetic transitions such as Gz Cy Fx ) Gz Cy ) Gz is observed according to
the total energy differences. Non trend in the weak-FM as a function of the
strain was observed even when the Co–F–Co bonding exhibit a change from �
to �� going from -�% to +�% strain respectively. The full understanding of the
non-collinear magnetism and its origin is still unclear and needs further studies.
SPIN-TEXTURE IN SrBO3 (B) = Ti
D
A N D H f ( � � � ) S U R FA C E S
Oxide perovskites exhibit a very large variety of physical properties, which allow
their use in a wide range of applications such as spintronics or thermoelectricity
[���, ���, ���]. In particular, during the last decade it has been observed new and
exotic functional properties at the surface of oxide crystals with no parallel to
their bulk counterpart [���, ���]. In this field of multifunctional oxide materials,
SrTiO3 (STO) has attracted a lot of attention due to its wide range of physical
and chemical properties that goes from superconductivity [���] to engineered fer-
roelectricity [��]. In various studies a two-dimensional electron gas (�DEG) has
been reported at the (���) surface of STO crystal, both theoretically[���] and ex-
perimentally by using angle-resolved photoemission spectroscopy (ARPES)[���,
���], which appeared to be induced by the presence of oxygen vacancies and ex-
plained by an electronic charge redistribution at the surface. The dynamics of the
Ti–d electron quasiparticle on this surface have been experimentally character-
ized and the appearance of cooperative effects of the electrons has been reported,
giving rise to an orbital ordering itself leading to a Rashba-like splitting at the
dxy - dxz/yz bands crossing close to the Fermi level[���]. Recently, a mixed char-
acter of the confined electrons of the (���) surface of STO displaying a complex
spin texture at the Fermi surface has been reported. This unexpected result was
ascribed to a Zeeman interaction plus Rashba-like splitting of the lowest dxy owed
as well to oxygen vacancies[���].
As commented previously, similar Rashba-effect[���] in oxides has been ob-
served in the LaAlO3 /SrTiO3 interface [���, ���, ���, ���] or in SrTiO3 - and
KTaO3 -based transistors [���, ���]. For these single perovskite-based devices, a
k-cubic dependence of the splitting was found in contrast to the linear+cubic
dependence of the Rashba splitting at the LaAlO3 /SrTiO3 interface due to the
dxy – dxz/yz multiorbital nature of the lowest bands [���]. As in the case of
Au(���) surface, where it has been shown that even when the surface states ex-
hibit a mainly p–orbital character, the d–orbitals drives the Rashba splitting
and dictates the direction of the spin by means of the orbital-angular momen-
tum (OAM) [���, ���]. Therefore, the d-orbitals of the transition metal atom of
perovskite oxides can be the key parameter for a large and controllable electron
band splitting in future spintronic devices [���, ���]. This later property gives an
additional degree of freedom to the rich field of multifunctional perovskites. How-
ever, the recent findings of a giant spin-splitting at the STO surface, opens some
questions about the nature of the spin-splitting in this type of materials, such as:
What is the origin of this giant-splitting? Could it be manipulated or observed
���
��� spin-texture in s r bo 3 (b) = t i and h f (���) surfaces
12Å Vacuum
O
Ti
OV Sr
Octahedral
AFD-distortions
y
x’
y’
x
z z
y x
Figure D.�: Schematic view of the SrTiO3 : TiO2- (���) surface used in the calcula-
tions with an oxygen vacancy (OV) at the yz- (left) and xz-plane (right)
respectively. The OV is located at the ( 14 , 14 , z) position with respect to the
computational x 0 - and, y 0 - coordinates). The amplification of the AFD dis-
tortions due to the oxygen vacancy can be clearly appreciated. Additionally,
a surface induced off-centering of the Ti along the z direction is observed
with respect to the oxygens (rumpling), inducing an electric dipole. In the
right bottom part we present a top view of the unit cell in the xy-plane.
This figure was elaborated with the support of the VESTA code [���].
D.� computational details ���
SrTiO3 SrHfO3
sx sx
Energy [eV]
Energy [eV]
sy sy
sz sz
kx ky kx ky
Figure D.�: Band structures around -point for SrTiO3 and SrHfO3 with an oxygen
vacancy on the BO2 terminated surface. Red and blue denote positive and
negative projection of spin vector ~s on Cartesian axes, units of the k-vector
in units of ⇡/a. The non spin-polarized bands are depicted in dark grey.
The effect of the SOC and a Rashba-type splitting is observed in the dxy
and dxz/yz bands crossing (noted with the dashed circles) as well as the
vertical energy shift ⌘ EZeeman + El . For SHO a spin-splitting of the
lower bands is observed for the x and y components around E = –�.� eV.
Figures were created by using the PyProcar script [���].
Ti–�d orbitals present in the TiO6 octahedra are expected to split into t2g (dxy ,
dxz , dyz ) and eg (dz2 , dx2 -y2 ) due to the crystal-field splitting. Here, we remark
that small distortions in the octahedral geometry, as the one observed at the sur-
face, generate a small shifting in the energy-values of the electron levels ( El ) and
a possible mixing of the t2g and eg states, mainly due to the symmetry breaking
and the uncompensated TiO6 octahedra at the surface [���]. In the oxygen va-
cancy STO surface we found that the lower d–levels bands corresponds to the dxy
orbitals in the paths that runs from -�.�kx – – +�.�kx and -�.�ky – – +�.�ky .
In addition, we also found a magnetic moment on the Ti atoms close to the sur-
face. The amplitude of these magnetic moments is around �.��� µB /Ti-atom in
the surface and �.��� µB /Ti-atom in the sub-surface layer. These results are in
agreement with the experimental observation of magnetism in STO with oxygen
vacancies [���] and other theoretical works of the same system [���, ���, ���].
The spin polarization is analyzed by plotting the sx , sy and sz components of the
magnetization for each path, where up and down spin orientations are depicted
by red and blue colors respectively. We can appreciate that lower light bands close
to the Fermi level are separated vertically due to an induced Zeeman splitting
( EZeeman ) as expected from the magnetism found at the surface and the en-
ergy splitting El . This magnetic moment is induced by the effect of the electron
surface density rearrangement as a response of the uncompensated orbitals, thus
oxygen vacancy acts as a magnetic impurity [���, ���].
From Fig. D.� we can also see the parabolic-like shape of the d-bands, which is
a characteristic of a �DEG, as reported in previous studies [���, ���, ���]. Here,
it is important to note that this �DEG is induced by the electronic reconstruction
related to the confinement of the electrons at the surface [���, ���], rather than
a polar catastrophe origin as mentioned for the SrTiO3 /LaAlO3 interface [���].
Thus, we can expect that the behavior of the conducting electrons in STO differs
from the one found in the oxide/oxide interfaces.
The effect of the SOC can be clearly appreciated in the Rashba-like spin-
splitting at the dxy and dxz/yz crossing (pointed out by the dashed circles in
Fig. D.�), similarly to what has been previously reported in this surface [���].
For the dxy state, we obtain an effective mass of �.� me for the orbital located at
the second layer of TiO2 . This calculated effective mass is two times larger than
the one reported in Ref.[���] (�.�� me ), which is too small for d like states [���].
However, a spin-splitting such as in a classical Rashba picture is not observed (i.e.
crossing of the spin up and down bands) between the light branches or a spin
inversion within the dxy bands for any of the spin-orientations along the kx and
ky paths. This result is in contrast with the interpretation of the experimental
measurements of Santander-Syro et al. [���] who reported a giant spin-splitting
in the lowest d-bands of the (���) STO surface. One possible explanation of this
splitting is that it can be a consequence of the surface ionic vibrations. In that
respect, we have performed ab–initio molecular dynamics calculations at constant
temperature of ��� K for the slab geometry. Selected uncorrelated frames were
taken from the simulation and electronic band structures were calculated for each
��� spin-texture in s r bo 3 (b) = t i and h f (���) surfaces
case. As for the � K ground state, we did not observe any Rashba splitting like
the one reported in Ref [���], which indicates that at least for low temperatures
there is no coupling between the ionic vibrations and the electronic response.
Another possibility could be that it comes from anharmonic vibrations at much
higher temperatures, a possibility we do not explore in this work. The absence
of a large spin-splitting can be understood from the small SOC present in the
�d–Ti orbitals. To test this hypothesis, we performed a computer experiment
by replacing all the Ti atoms by a heavier cation with much larger SOC such
as Hafnium. Taking the advantage that the Hf atom has the same electronic
configuration as the Ti but with �d orbitals � . In these calculations cell lattice
parameters and internal coordinates of the optimized geometry of the STO surface
were fixed in order to keep the same structure. The �DEG generated at the
HfO2 surface with oxygen vacancies presents the same main features as in the
STO case. However, a large change in the spin-splitting state is observed in the
electron dispersion as shown in Fig. D.�, and somehow closer to the experimental
suggestion of A. F. Santander-Syro et al. [���], who claimed for a spin-inversion
direction within the lowest d-band. In both, kx and ky paths in Fig. D.� and at
the Fermi level the bands are mainly composed by sz spin direction. The spin
switching is observed for the light dxy Hf band at E = –�.� eV and reciprocal
points close to –�.� and �.� for sx and sy spin components respectively. This
suggests that a strong surface polarization is not enough to induce a Rashba-like
spin-splitting and the presence of high-Z cations in the structure is also necessary
to achieve this spin state as previously reported for Ag, Cu, and Au surfaces [���].
The Rashba-like splitting in the SrHfO3 (SHO) is only observed for the sx and
sy spin components of the xy-plane perpendicular to z-axis as expected from the
Rashba Hamiltonian HR = ↵R (~ez ⇥ ~k) · ~s [���].
In addition to the band analysis, we have also extracted the in-plane spin-
texture as shown in the Fig. D.�. In Fig. D.� we can see two concentric rings for
up and down spin polarization respectively for the STO (at E = � eV) and one
single ring for SHO (E = –�.� eV) surface. The rings in the STO appears as a
result of the vertical splitting due to EZeeman + El , where El is induced by
the lowering of the dxy energy bands in comparison to dyz and dzx orbitals as
discussed before. In the case of SHO (Fig. D.�), spin precession in the lower ring
is observed in addition to a spin switching of the dxy band in agreement with the
bands in Fig. D.�) as expected for the Rashba-type splitting.
The absence of the spin-splitting in the STO surface could be also linked to
magnetism, as observed at the the Gd(����) surface [���] where the presence of
magnetism and a Rashba-type splitting has been reported. For this system the
splitting is lower than �.� Å-1 even with the large SOC present in Gd (around
� times larger than in Ti)[���]. Moreover, due to the presence of EZeeman ,
there is no crossing of the bands and the direction of the magnetization vector
(M
~ ) determine the presence of the spin-splitting. This is due to the fact that the
� The spin-orbit coupling in Hf is around �� times larger than in Ti based in the (Ze2 /hc)2
expression
D.� ST O : T iO 2- surface with oxygen vacancies ���
sx sy
SrTiO3 sz
ky ky ky
kx kx kx
SrHfO3
sx sy sz
ky ky ky
kx kx kx
Figure D.�: Band structure at the kx -ky plane (in ⇡/a) where the spin-texture is cal-
culated for each spin component along the cartesian axis (i.e. sx , sy and
sz ). For SrTiO3 compound (above, E = � eV) simple concentric rings for
up (red) and down (blue) spin polarizations are observed. In contrast a
Rashba-like splitting is observed by replacing the Ti cation by the heavier
Hf in the SrHfO3 (bottom, E = – �.� eV). Figures were created by using
the script PyProcar [���].
Rashba-splitting vanishes for ~k k M~ and only EZeeman remains [���]. The later
statement needs to be carefully taken into account in the STO case for which we
find a magnetic moment components at the surface mx =�.���, my =�.��� and
mz =�.��� µB /Ti-atom respectively and thus after considering several possible
magnetic orientations. Different starting magnetic configurations, along x- y- and
z-axis respectively were taken into account and full structural and energy mini-
mizations were applied. Nonetheless, the system always relaxed toward the same
magnetic ground state in which a diagonal direction of the magnetic moment is
observed with respect to the surface. No traces of magnetization were observed
at the oxygen sites confirming the Ti–localized magnetic moments. Thus, this
rather large in- and out-of-plane magnetic moments have the tendency to pre-
vent the possibility to have a giant Rashba splitting, in addition to the small
SOC discussed previously.
In the case of SHO, we obtain that the Hf magnetic moments are strongly re-
duced to mx =�.���, my =�.��� and mz =�.��� µB /Hf–atom respectively close
to the surface, which strongly reduced the EZeeman . Fig. D.� shows that the
lowest d-band in SHO only exhibit spin components along the xy-plane. Nonethe-
less for higher energy values, the bands are mainly spin-polarized along the z-axis,
demonstrating a very complex spin structure with a unique spin-texture. There-
��� spin-texture in s r bo 3 (b) = t i and h f (���) surfaces
fore, only a deep understanding of this subtle competition between the polar
distortions, the symmetry conditions, SOC, and magnetism would allow us to
tune the Rashba-splitting at the (���) BO2 surface.
In order to test the validity of our results against the exchange-correlation
functional, we also carried out the same calculations with both the Local Den-
sity Approximation [���] (LDA) and the GGA PBEsol+U functionals. Using the
LDA approximation, we did not find substantial differences with the GGA PBEsol
calculations. Due to the strong correlation of the d-orbitals that drives the mag-
netism, we might expect the DFT+U formalism [���] to enhance the electron cor-
relation and change the physical picture discussed previously. Indeed, above U =
� eV, we observe an insulating state of the STO surface with the oxygen vacancy,
which do not agree with the experimental measurements [���, ���, ���]. Then,
for whatever U value below � eV, we do not observe a large Rashba+Zeeman
(RZ) type spin-splitting. In order to test the oxygen vacancy position effects, we
performed a second set of calculations considering the oxygen vacancy in the sec-
ond inner layer from the surface instead of the first layer and, as before, we did
not observe a Rashba spin-splitting in the STO surface within this configuration.
sx OV
Energy [eV]
sy
sz
ky
kx ky
Figure D.�: Band structures around the -point along kx and ky for SrHfO3 surface
with an oxygen vacancy on the top of the HfO2 layer for each direction of
the spin. Red and blue colors represent the up and down spin directions
respectively. A larger Rashba-type splitting of the bands is observed for
the z component as appreciated for the ky – sz extracted image. The SHO
slab structure is also shown in which the large octahedral tilting can be
observed. Figures were created by using the script PyProcar [���].
�.�� µC·cm-2 found along the x-axis (Px ) and y-axis (Py ) respectively. This in-
plane polarization in orthorhombic perovskites is quite well known and is related
to the presence of anti-polar motion of the A-site in the xy-plane, which, in the
case of the surface symmetry breaking, do not cancel out from layer to layer
and thus resulting in an effective spontaneous polarization along the xy-plane
[���, ���]. Interestingly, and in contrast to the STO case, we do not find any
presence of magnetism in the Pbnm SHO surface with oxygen vacancies and at
the same time the Rashba-splitting is considerably enhanced, confiming the rule
discussed above. Interestingly, when replacing the Hf by Ti in the relaxed Pbnm
SrHfO3 surface, the magnetism appears again and the Rashba splitting dissapear.
��� spin-texture in s r bo 3 (b) = t i and h f (���) surfaces
All of that confirms the influence of the magnetic ordering and the SOC strength
of the B-site at the spin-texture of SrBO3 surfaces with oxygen vacancies.
���
��� role of surface water adsorptions in s r t i o 3 /l a a l o 3 interface
Figure E.�:
Figure 5. Simulated lattice distortions
Relaxed geometries in 3uc LaAlO
for the pristine, 3/SrTiO
and the heterostructureand
OH, 3 �.�H+�.�OH, with
H
different surface H and OH densities. Tiltation of oxygen octahedra becomes more
groups. It can be observed the octahedral tilting and rotation increase when
significant asthe
themolecules
surface proton density
are placed atincreases.
the surface.
We have noticed that the octahedral rotations and tilts across the LAO layer
—dissipating up to the STO substrate— increases when H and OH are placed
at the LAO surface (see Fig E.�). These AFD distortions have a direct effect on
the interfacial properties and specifically, at the carriers density at the interface
as we have observed in a computational experiment where, the carried density
was computed with and without full atomic relaxations. Here, it is important to
remark that those distortions are present even in the pristine bulk structure due
to the rhombohedral R�c ground state of the LaAlO3 [���].
In the �uc pristine case, the electric field of �� mV/Š(computed from the layer
resolved DOS) is clearly not enough to overcome the existing electron gap in the
STO/LAO system, leading to the obtained insulated behavior (Eg = �.� eV).
The water dissociation and the later removal of the (OH)- groups with the AFM
tip leaves only absorbed hydrogen atoms on the surface which, according to our
results, act as electron-donors. These additional electrons, under the effect of the
electric field existing in the LAO film, are transferred to the interface, which then
generate a �DEG with dxy character from the Ti orbitals. This electron transfer
explains the induced IMT in the �uc LAO, as depicted in Fig. E.�. Consequently,
the absorbed hydrogen act as electron-donor and provide additional charge, which
is enough to compensate the surface charge responsible for the polar discontinuity.
The later clearly alters the out-of-plane electrostatic potential. For example, when
��� role of surface water adsorptions in s r t i o 3 /l a a l o 3 interface
the surface is fully hydrogen saturated, the field takes a value of -�� mV/Šthat
is in opposite direction to the potential in the (AlO2 )- pristine surface. This field
crosses through � mV/Å for �.�H (see Fig. E.�), which demonstrates the direct
effect of the electron-donors in the electrostatic response, as was expected. On the
other hand, if we consider �uc of LaAlO3 on SrTiO3 (���), the electric field is now
close to ��� mV/Šacross the entire LAO film as shown in Fig. E.�. This value is
sufficient to overcome the present electronic band gap of the STO material and
the system becomes metallic, as it is widely known. Additionally, a �DEG and a
�DHG are created at the interface and surface respectively. Now, with the same
process of water absorption and later (OH)- removal, as in the previous �uc case,
the hydrogen electron-donor alters the amount of electrons in the system. These
additional electrons are transferred from the surface to the interface, similarly to
the case of �uc of LaAlO3 , and the amount of carriers density increase drastically
as observed in Fig. E.�, as it is also supported from experimental results.
Figure E.�: Layer-by-layer density of states for the �uc in which the electric field profile
can be observed. A insulator-to-metal transition can be observed when going
from the pristine to the H-doped structures.
Figure E.�: Layer-by-layer density of states for �uc for the pristine structure and the
water dissociated at the surface.
E.� water dissociation in the surface ���
Our theoretical analyses also reveal the appearance of magnetic moments at the
interface for the �uc case, as expected from the unpaired �d-Ti orbitals occupation
as it has been already reported [���]. In the �uc case, magnetic moments are only
observed at the �d-Ti when the interface becomes metallic due to the bonded
hydrogen. When only (OH)-groups are in the LAO surface, the magnetism is
only present at the top layer. This is due to the electron holes at the localized
oxygen p-orbitals as it has been shown in some d0 semiconductor [���].
The fact we are considering oxides in which, strong electron correlations takes
place, makes important to investigate the effect in the results of the exchange
correlation functional by enhancing the electron-electron repulsion in the so called
DFT+U [��]. This method allows to increase the electron correlation lacking in
the d-orbitals and then improving the electron band gap, in closer agreement with
experimental predictions. This correction can also affect negatively the searched
properties because it changes at the same time some other orbitals. Thus, we have
considered initially the �uc case as a prototypical example. All results here were
re-computed with a U = �.� and �.� eV for Ti and La respectively and J = � eV
[���]. We found that the general description holds but the valence carrier density
increases by more than one order of magnitude, which is not in agreement with
the reported experiments. Therefore, our results reported here only correspond
to the case of U = J = � eV.
Figure E.�: Computed carrier density for �uc and �uc respectively. Here we considered
several molecules as a dopant agent at the surface.
Now, if we think in the idea to reverse the voltage in the tip, and then, the
H+ is removed from the surface, the effect can be seen in the reduction of the
electric field in the LAO film. This field goes from ��� mV/Što �� mV/Šfor
the pristine structure to the OH saturated surface. Nonetheless, according to our
results, this field reduction is not having a major effect in the carrier density as
can be appreciated in Fig. E.�. This can be explained by the fact that in the
�uc, the charge is already confined at the interface and there is no charge transfer
going back to the surface, which is traduced in an increase of the holes at the
��� role of surface water adsorptions in s r t i o 3 /l a a l o 3 interface
surface. In the �uc case, the electric field is also reduced to �� mV/Šand the
system remains insulator.
We considered the possibility of reversing the voltage on the AFM tip, H+ is
removed and OH- are kept on the surface, such that there is a reduction of the
electric field in the LAO film. This field goes from ��� mV/Što �� mV/Šfor the
pristine structure to the OH saturated surface respectively. However, according
to our results, this field reduction is not having a major effect in the carrier
density as can be appreciated in Fig. E.�. This can be explained by the fact that
in the �uc, the charge is already confined at the interface and there is no charge
transfer going back to the surface, which is already saturated with OH. The later
is traduced in an increase of the holes at the surface. In the �uc case, where the
heterostructure is already an insulator, the electric field is also reduced to ��
mV/Å and the system kept a finite electronic band gap.
LIST OF FIGURES
and R+ 4 . modes. . . . . . . . . . . . . . . . . . . . . . . . . ��
Figure �.� Polar 4- mode, into the Bethe’s notation [��], as a func-
tion of the A- and B-site in ABF3 perovskites, the imagi-
nary modes are presented as negative values. . . . . . . . ��
Figure �.� Ferroelectric instability in LiBF3 and NaBF3 (B = Ca,
Mn, V, Zn, and Ni) as a function of the tolerance factor, t. ��
Figure �.� Selected IFC interactions in the �⇥�⇥� supercell. A- B-
and F-sites are depicted in green, violet and grey respec-
tively. Additionally, the atoms are labeled to correlate the
IFCs shown in the Table �.�. . . . . . . . . . . . . . . . . ��
Figure �.� Double-well energy obtained at the FE-mode freezing in
NaBF3 compounds. The gain in energy is in agreement
with previous results that shows that the bigger is the
B-site, bigger is the FE instability and with this, the pos-
sibility to be ferroelectric. ± denotes the full amplitude
of the freezing. . . . . . . . . . . . . . . . . . . . . . . . . ��
Figure �.� Double-Well energy curves and schematic representation
of atomic displacements in different models of ferroelec-
tricity for oxides a) BaTiO3 and b) NaMnF3 . . . . . . . . ��
Figure �.� Unstable FE mode at the cubic Pm3̄m structure in ABX3
compounds. . . . . . . . . . . . . . . . . . . . . . . . . . . ��
Figure �.� Eigendisplacements for NaCaF3 compound showing the
⌘ transition also found in B-site: V, Mn, and Zn (left).
Schematic representation of eigendisplacement in the flu-
orides compounds for strain values above and below of the
relaxed structure (right). . . . . . . . . . . . . . . . . . . . ��
Figure �.� Eigendisplacements for BaTiO3 showing the change of the
mode as cause of suppression of ferroelectric instability as
a function of isotropic pressure. . . . . . . . . . . . . . . . ��
���
��� List of Figures
Figure �.�� On-site IFC for A-, B- and F-sites (x=y, and z ) as a func-
tion of strain for each of sites in NaBF3 with B = Ca, Mn,
V, Zn, and Ni. . . . . . . . . . . . . . . . . . . . . . . . . ��
Figure �.�� Contribution of modes in the phase transition from Pm3̄m
to Pnma in the NaBF3 fluoroperovskites. . . . . . . . . . ��
Figure �.�� Behavior of FE and R-point instabilities as a function of
hydrostatic pressure for NaNiF3 . . . . . . . . . . . . . . . ��
Figure �.�� Eigendisplacements for all polar modes in NaNiF3 . . . . . ��
Figure �.�� Squared frequency at the vibrational instabilities. We can
observe the FE-mode at and AFD-modes at R-point as
a function of hydrostatic strain for LiNiF3 (left). Eigendis-
placements of FE-mode instability in LiNiF3 are also shown
(right). . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ��
Figure �.� (a) Total energy of NaMnF3 as a function of the volume
per atom for all different crystal structures. Cmcm and
Pnma structures show the lowest energy values for lower
volumes. (b) H/atom as a function of external pressure in
NaMnF3 structures. Pnma to Cmcm structural transition
can be seen at the critical pressure of � GPa. . . . . . . . ��
Figure �.� (a) Schematic view of the phase transition from Pnma to
Cmcm. Na, Mn and F atoms are represented in green, vio-
let and grey colors respectively. The breaking of bondings
for the transition from corner-shared to edge-shared oc-
tahedra can be appreciated. (b) Simulated XRD patterns
for Pnma and Cmcm structures. The room temperature
Pnma phase is in good agreement with the experimental
reports (open triangles) [���, ���, ���]. (c) Octahedra tilt
angle as a function of pressure in NaBF3 . For B = Zn,
Ni, Co and Mg experimental results are taken from Yusa
et al. [���]. All of these compounds present the transi-
tion from the Pnma to the Cmcm structures at a tilting
angle close to �� . Experimental results for NaMnF3 at
low pressure were taken from Katrusiak et al. [��] Full
agreement between experimental (filled squares) and our
theoretical calculations (empty squares) can be appreciated. ��
Figure �.� Vibrational behavior as a function of isotropic pressure for
Pnma and Cmcm pPv. The open and filled symbols are
the Raman and the Infra-red modes respectively. . . . . ��
List of Figures ���
Figure �.� (a) Energy differences between the different magnetic or-
derings of the pPv phase of NaMnF3 versus the U pa-
rameter. (b) Sketch of the Cy Fz and the Cz Ay magnetic
orderings in the pPv phase of NaMnF3 . For U < �.� eV
the system exhibits a AFM+FM magnetic state with a
high ferromagnetic canting. For U > �.� eV, the system
transit toward a canted AFM state with no weak FM. . . ��
Figure �.� Magnetic canting dependence as a function of U and J pa-
rameters for NaMnF3 . Canting angle (') varies from �.��
to ��.�� when U goes from �.� to �.� eV. The canting
angle goes from ��.�� to �.�� when J varies from �.� to
�.� eV. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ��
Figure �.� Projected local density of states (LDOS) of the Cmcm
phase of NaMnF3 . A regular F-p bonding Mn-d anti-
bonding is observed with small Mn:3d – F:2p overlaping
DOS is responsible for the small superexchange interac-
tion present in NaMnF3 . . . . . . . . . . . . . . . . . . . . ��
Figure �.� Simulated XRD patterns with K↵:Cu for CaRhO3 (above)
and NaNiF3 (below) in the Cmcm phase at atmospheric
pressure are observed. . . . . . . . . . . . . . . . . . . . . ��
Figure �.� The non-collinear magnetic orderings in the pPv Cmcm
symmetry labeled as Cz Fy , Cz Ay , Cy Fz and Gy Az are
presented in the figure. Two of them with pure AFM
(Cz Ay and Gy Az ) behavior and two AFM+weak-FM (Cz Fy
and Cy Fz ) can be observed (Fig. elaborated with the VESTA
software [���]). . . . . . . . . . . . . . . . . . . . . . . . . ��
Figure �.� Schematic representation of the atomic sites and magnetic
interactions in ABX3 (A, B and X -sites in blue, yellow
and red respectively). Exchange constants labeled as Ja ,
Jc and Jab and DM interactions represented as Da and Dc
equivalent to D12 and D13 respectively are shown. (Fig.
elaborated with the VESTA software [���]). . . . . . . . . ��
Figure �.� Ferroelectric instability T1u at the cubic Pm3̄m as well as
the frequency values for the softest B2u mode in the Pnma
NaBF3 (B = Mn, V, Zn, and Ni) compounds. . . . . . . . ��
Figure �.� a) Epitaxially strain induced B2u instability in Pnma NaMnF3
in which a nonlinear coupling is observed, imaginary fre-
quencies are represented as negative values by notation. b)
Eigendisplacements of the polar B2u along b-axis (⌘b ) for
each site (Na in green squares, Mn in violet circles and F? ,
Fk in upper and down triangles respectively) in NaMnF3
in Pna�1 . c) Ferroelectric mode freezing in which the tran-
sition from the paraelectric Pnma to a polar Pna�1 phase
is observed at equivalent �% in the ac-plane. . . . . . . . ��
��� List of Figures
Figure �.� Localized-DOS for the �d-Ni and �d-Ta at the KNiF2.5 O0.5
and KTaO2.5 F0.5 oxyfluorides respectively. This DOS was
computed by means of the hybrid functional HSE06. . . . ��
Figure �.� Electronic configuration corresponding and crystal-field
splitting of the KTaO3 oxide and the KTaO2.5 F0.5 oxyflu-
oride. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ��
Figure �.� Electronic and magnetic configuration corresponding to
the crystal-field splitting for the KNiF3 fluoride and the
KNiF2.5 O0.5 oxyfluoride. . . . . . . . . . . . . . . . . . . ��
Figure �.� �/� and �/� superlattices at which the insulator-to-metal
transitions takes place. The TaO2 /KF and ZnF2 /KO in-
terfaces can be observed. . . . . . . . . . . . . . . . . . . . ���
Figure �.� Atomic displacements of the B-sites along the z-axis (left)
with respect to the oxygen apical plane for the �/� super-
lattice. K–K c-lattice parameter in which, a large expan-
sion and compression can be observed at the TaO2 /KF
and ZnF2 /KO interfaces respectively (right). . . . . . . . ���
Figure �.� Localized density of states for M = Zn and Ni in the �/�
superlattice ratio obtained considering the HSE06 hybrid
functional. . . . . . . . . . . . . . . . . . . . . . . . . . . . ���
Figure �.� Layered revolved DOS for the KTaO3 /KZnF3 and KTaO3 /KNiF3
(�/�) superlattices where a �DEG and �DHG are observed
at the TaO2 /KF and MF2 /KO. . . . . . . . . . . . . . . . ���
Figure �.� Band-edge diagram for the electrostatic behavior of the
system where the electric field profile and the �DEG (�DHG)
formation are presented at the oxyfluoride TaO2 /KF (ZnF2 /KO)
interface. This was obtained based on the layered resolved
DOS (see Fig �.�). CB, VB an EF denotes the conduction-
band, valence-band and the Fermi energy respectively. . . ���
Figure �.�� Magnetic moment per layer at the KTaO3 /KNiF3 super-
lattices. A moment modulation can be appreciated when
n/l ratio is varied from �/� to �/�. . . . . . . . . . . . . . ���
Figure �.�� Bands structure for the p- (left) and n-type (right) inter-
faces along the X–M– –X path in the Brillouin zone for
the KTaO3 /KZnF3 system. The formation of the �DHG
and �DEG at the p- and n- interfaces can be appreciated. ���
Figure �.�� Bands structure for the p- (left) and n-type (right) in-
terfaces along the X–M– –X path in the Brillouin zone
for the KTaO3 /KNiF3 system. The �DHG and �DEG at
the p- and n- interfaces can be appreciated. Besides, and
empty Ni:eg band is observed between �.� to �.� eV re-
spectively. . . . . . . . . . . . . . . . . . . . . . . . . . . . ���
��� List of Figures
Figure D.� Band structures around -point for SrTiO3 and SrHfO3
with an oxygen vacancy on the BO2 terminated surface.
Red and blue denote positive and negative projection of
spin vector ~s on Cartesian axes, units of the k-vector in
units of ⇡/a. The non spin-polarized bands are depicted
in dark grey. The effect of the SOC and a Rashba-type
splitting is observed in the dxy and dxz/yz bands crossing
(noted with the dashed circles) as well as the vertical en-
ergy shift ⌘ EZeeman + El . For SHO a spin-splitting
of the lower bands is observed for the x and y components
around E = –�.� eV. Figures were created by using the
PyProcar script [���]. . . . . . . . . . . . . . . . . . . . . ���
Figure D.� Band structure at the kx -ky plane (in ⇡/a) where the spin-
texture is calculated for each spin component along the
cartesian axis (i.e. sx , sy and sz ). For SrTiO3 compound
(above, E = � eV) simple concentric rings for up (red) and
down (blue) spin polarizations are observed. In contrast a
Rashba-like splitting is observed by replacing the Ti cation
by the heavier Hf in the SrHfO3 (bottom, E = – �.� eV).
Figures were created by using the script PyProcar [���]. . ���
Figure D.� Band structures around the -point along kx and ky for
SrHfO3 surface with an oxygen vacancy on the top of the
HfO2 layer for each direction of the spin. Red and blue col-
ors represent the up and down spin directions respectively.
A larger Rashba-type splitting of the bands is observed for
the z component as appreciated for the ky – sz extracted
image. The SHO slab structure is also shown in which
the large octahedral tilting can be observed. Figures were
created by using the script PyProcar [���]. . . . . . . . . ���
Figure E.� Relaxed geometries for the pristine, and the OH, �.�H+�.�OH,
and H groups. It can be observed the octahedral tilting
and rotation increase when the molecules are placed at
the surface. . . . . . . . . . . . . . . . . . . . . . . . . . . ���
Figure E.� Layer-by-layer density of states for the �uc in which the
electric field profile can be observed. A insulator-to-metal
transition can be observed when going from the pristine
to the H-doped structures. . . . . . . . . . . . . . . . . . . ���
Figure E.� Layer-by-layer density of states for �uc for the pristine
structure and the water dissociated at the surface. . . . . ���
��� List of Figures
Figure E.� Computed carrier density for �uc and �uc respectively.
Here we considered several molecules as a dopant agent at
the surface. . . . . . . . . . . . . . . . . . . . . . . . . . . ���
L I S T O F TA B L E S
Table �.� Available physical quantities through the first- and second-
order energy derivatives with = ⌧, ⌘, and E, (see the Eq.
�.��). Table adapted from Ref. [��]. . . . . . . . . . . . . ��
Table �.� Electronic configuration of pseudopotentials. Valences and
semicore electrons taken into account for calculations. . . ��
Table �.� Summary of most unstable high symmetry point phonon
eigenvalues, relaxed lattice constants for cubic Pm3̄m sym-
metry, ionic A and B sites radii and reported ground state
symmetry for all compounds [��]. HS: High Spin, LS: Low
Spin. All modes in ( 4- , R+4 , R5 , X5 and M3 ) in cm
+ + + -1 . . ��
Table �.� Eigendisplacements (⌘) for all compounds at relaxed cubic
structure. X? and Xk are the sites located perpendicular
and parallel to the M-F bonding respectively . . . . . . . ��
Table �.� Born Effective Charges (Z⇤ , e- units) in for all compounds
computed at relaxed cubic structure computed according
to the Eq. �.�. The X? and Xk are the components per-
pendicular and parallel to the B–X bonding. . . . . . . . . ��
Table �.� On-site interatomic constants (on-site IFC, in eV/Å2 units)
for the ABX3 flouroperovskites. The on-site IFC for the
BaTiO3 and PbTiO3 are also presented for comparison.
The Fx=y sites lies in A-F plane meanwhile the Fz lies
along B–F bonding. . . . . . . . . . . . . . . . . . . . . . . ��
Table �.� Selected IFC according to interactions illustrated in the
Fig.�.� for the A-, B-, and F-sites respectively. . . . . . . ��
Table �.� On-site IFC for all NaMF3 compounds at same lattice
parameter equivalent to expansion. . . . . . . . . . . . . . ��
Table �.� Lattice parameters of the pPv phase (a, b, and c) and
related temperature (T) and pressure (P) transition pa-
rameters from the Pnma phase to the pPv phase of fluo-
ride and oxide perovskites. In the pPv phase, the A and
B cations occupy the 4a and 4c Wyckoff positions and the
O/F anions occupy the 4c and 8d positions. Most of the
oxides have higher pressure and temperature transitions
than the fluorides. . . . . . . . . . . . . . . . . . . . . . . ��
Table �.� Calculated Raman, Infra-red (IR) and silent modes fre-
quencies (cm-1 ) of the Pnma phase of NaMnF3 (at zero
pressure). In brackets are presented the experimental val-
ues of the Raman [���] and the IR [���] modes. . . . . . ��
���
��� List of Tables
���
��� references
[��] D. S. Sholl and J. A. Steckel, “Dft calculations for simple solids,” in Density
Functional Theory (John Wiley and Sons, Inc., ����) pp. ��–��.
[��] M. Fiebig, Journal of Physics D: Applied Physics ��, R��� (����).
[��] A. Poole, B. Roessli, O. Zaharko, and K. W. Kramer, Journal of Physics:
Condensed Matter ��, ������ (����).
[��] G. Nenert, Journal of Physics: Condensed Matter ��, ������ (����).
[��] G. Nenert and T. T. M. Palstra, Journal of Physics: Condensed Matter ��,
������ (����).
[��] W. Cochran, Phys. Rev. Lett. �, ��� (����).
[��] P. Ghosez, X. Gonze, and J.-P. Michenaud, EPL (Europhysics Letters) ��,
��� (����).
[��] R. Cohen, Nature ���, ��� (����).
[��] P. Ghosez, J.-P. Michenaud, and X. Gonze, Phys. Rev. B ��, ���� (����).
[��] P. Ghosez, E. Cockayne, U. V. Waghmare, and K. M. Rabe, Phys. Rev. B
��, ��� (����).
[��] M. Guennou, P. Bouvier, G. Garbarino, J. Kreisel, and E. K. H. Salje,
Journal of Physics: Condensed Matter ��, ������ (����).
[��] P. Edwardson, L. Boyer, R. Newman, D. Fox, J. Hardy, J. Flocken, R. Guen-
ther, and W. Mei, Physical Review B ��, ���� (����).
[��] Z. Friedman, M. Melamud, J. Makovsky, and H. Shaked, Phys. Rev. B �,
��� (����).
[��] D. Babel, in Structure and Bonding (����).
[��] J. W. Flocken, R. A. Guenther, J. R. Hardy, and L. L. Boyer, Phys. Rev.
B ��, ���� (����).
[��] W. Zhong, D. Vanderbilt, and K. M. Rabe, Phys. Rev. B ��, ���� (����).
[��] D. Vanderbilt and W. Zhong, Ferroelectrics ���, ��� (����).
[��] N. A. Benedek and C. J. Fennie, The Journal of Physical Chemistry C ���,
����� (����).
[��] S. Amisi, E. Bousquet, K. Katcho, and P. Ghosez, Phys. Rev. B ��, ������
(����).
[��] A. Togo, F. Oba, and I. Tanaka, Phys. Rev. B ��, ������ (����).
��� references
[���] S. Yakovlev, M. Avdeev, and M. Mezouar, J. Solid State Chem. ���, ����
(����).
[���] F. L. Bernal, K. V. Yusenko, J. Sottmann, C. Drathen, J. Guignard, O. M.
Lovvik, W. A. Crichton, and S. Margadonna, Inorganic Chemistry ��,
����� (����), pMID: ��������, http://dx.doi.org/��.����/ic������y .
[���] F. Birch, Phys. Rev. ��, ��� (����).
[���] X. Gonze, Phys. Rev. A ��, ���� (����).
[���] X. Gonze, Phys. Rev. A ��, ���� (����).
[���] P. Daniel, M. Rousseau, A. Desert, A. Ratuszna, and F. Ganot, Phys. Rev.
B ��, ����� (����).
[���] A. Lane, D. Sharp, and J. Barraclough, J. Chem. Soc. ��, �� (����).
[���] V. Gava, A. L. Martinotto, and C. A. Perottoni, Phys. Rev. Lett. ���,
������ (����).
[���] J. R. Shane, J. Appl. Phys. ��, ���� (����).
[���] M. Akaogi, Y. Shirako, H. Kojitani, T. Nagakari, H. Yusa, and K. Yamaura,
Physics of the Earth and Planetary Interiors ���, ��� (����).
[���] S. Tateno, K. Hirose, N. Sata, and Y. Ohishi, Phys. Earth and Planetary
Interiors ���, �� (����).
[���] M. O’keeffe, B. Hyde, and J. Bovin, Phys. Chem. Minerals ���, ��� (����).
[���] C. D. Martin and J. B. Parise, Earth and Planetary Science Letters ���,
��� (����).
[���] B. Beznosikov and K. Aleksandrov, Crystallography Reports ��, ���
(����).
[���] D. Babel and A. Tressaud, in Inorganic Solid Fluorides, edited by P. Ha-
genmuller (Academic Press, ����) pp. �� – ���.
[���] M. Aroyo, L, J. M. Perez-Mato, C. Capillas, E. Kroumova, S. Ivantchev,
G. Madariaga, A. Kirov, and H. Wondratschek, Zeitschrift fur Kristallo-
graphie - Crystalline Materials ���, �� (����).
[���] M. I. Aroyo, A. Kirov, C. Capillas, J. M. Perez-Mato, and H. Wondratschek,
Acta Cryst ��, ��� (����).
[���] P. Anderson, Physical Review �� (����).
[���] J. Kanamori, Journal of Physics and Chemistry of Solids ��, �� (����).
references ���
[���] R. Zorin, B. Al’shin, and D. Astrov, Soviet Physics JETP ��, ��� (����).
[���] M. Bibes and A. Barthélémy, Electron Devices, IEEE Transactions on ��,
���� (����).
[���] S. Fusil, V. Garcia, A. Barthélémy, and M. Bibes, Annual Review of Ma-
terials Research ��, �� (����).
[���] L. Martin, Y.-H. Chu, and R. Ramesh, Materials Science and Engineering:
R: Reports ��, �� (����).
[���] N. Reyren, S. Thiel, A. D. Caviglia, L. F. Kourkoutis, G. Hammerl,
C. Richter, C. W. Schneider, T. Kopp, A.-S. Rüetschi, D. Jaccard,
M. Gabay, D. A. Muller, J.-M. Triscone, and J. Mannhart, Science ���,
���� (����).
[���] J. Shen, H. Lee, R. Valentí, and H. O. Jeschke, Phys. Rev. B ��, ������
(����).
[���] A. F. Santander-Syro, O. Copie, T. Kondo, F. Fortuna, S. Pailhès, R. Weht,
X. G. Qiu, F. Bertran, A. Nicolaou, a. Taleb-Ibrahimi, P. Le Fèvre, G. Her-
ranz, M. Bibes, N. Reyren, Y. Apertet, P. Lecoeur, a. Barthélémy, and
M. J. Rozenberg, Nature ���, ��� (����).
[���] P. D. C. King, S. McKeown Walker, A. Tamai, A. de la Torre, T. Eknapakul,
P. Buaphet, S.-K. Mo, W. Meevasana, M. S. Bahramy, and F. Baumberger,
Nature Communications �, � (����).
[���] A. F. Santander-Syro, F. Fortuna, C. Bareille, T. C. Rödel, G. Landolt,
N. C. Plumb, J. H. Dil, and M. Radović, Nat. Mater ��, ���� (����).
[���] M. Dawber and E. Bousquet, MRS Bulletin ��, ���� (����).
[���] N. C. Bristowe, P. Ghosez, P. B. Littlewood, and E. Artacho, Journal of
Physics: Condensed Matter ��, ������ (����).
[���] Z. Q. Liu, C. J. Li, W. M. Lü, X. H. Huang, Z. Huang, S. W. Zeng, X. P.
Qiu, L. S. Huang, A. Annadi, J. S. Chen, J. M. D. Coey, T. Venkatesan,
and Ariando, Phys. Rev. X �, ������ (����).
[���] A. H. Romero and F. Munoz, “Pyprocar,” http://sourceforge.net/p/
pyprocar/PyPROCAR/HEAD/tree/ (����).
[���] G. Khalsa, B. Lee, and A. H. MacDonald, Phys. Rev. B ��, ������ (����).
[���] W. D. Rice, P. Ambwani, M. Bombeck, J. D. Thompson, G. Haugstad,
C. Leighton, and S. a. Crooker, Nature materials ��, ��� (����).
[���] H. Choi, J. D. Song, K.-R. Lee, and S. Kim, Inorganic Chemistry ��, ����
(����).
[���] L. P. Gor’kov, Journal of Physics: Condensed Matter ��, ������ (����).
[���] C. Lin and A. A. Demkov, Phys. Rev. Lett. ���, ������ (����).
[���] W. Meevasana, P. D. C. King, R. H. He, S.-K. Mo, M. Hashimoto, a. Tamai,
P. Songsiriritthigul, F. Baumberger, and Z.-X. Shen, Nature materials ��,
��� (����).
[���] G. Hautier, A. Miglio, D. Waroquiers, G.-M. Rignanese, and X. Gonze,
Chemistry of Materials ��, ���� (����).
[���] F. Reinert, G. Nicolay, S. Schmidt, D. Ehm, and S. Hüfner, Phys. Rev. B
��, ������ (����).
[���] O. Krupin, G. Bihlmayer, K. Starke, S. Gorovikov, J. E. Prieto, K. Döbrich,
S. Blügel, and G. Kaindl, Phys. Rev. B ��, ������ (����).
[���] J. P. Perdew and A. Zunger, Phys. Rev. B ��, ���� (����).
[���] A. I. Liechtenstein, V. I. Anisimov, and J. Zaanen, Phys. Rev. B ��, R����
(����).
[���] B. J. Kennedy, C. J. Howard, and B. C. Chakoumakos, Phys. Rev. B ��,
���� (����).
[���] D. Vanderbilt, Journal of Physics and Chemistry of Solids ��, ��� (����).
[���] N. A. Benedek, J. M. Rondinelli, H. Djani, P. Ghosez, and P. Lightfoot,
Dalton Trans. ��, ����� (����).
[���] J. Sánchez-Barriga, A. Varykhalov, J. Braun, S.-Y. Xu, N. Alidoust, O. Ko-
rnilov, J. Minár, K. Hummer, G. Springholz, G. Bauer, R. Schumann, L. V.
Yashina, H. Ebert, M. Z. Hasan, and O. Rader, Phys. Rev. X �, ������
(����).
[���] C.-H. Park and S. G. Louie, Phys. Rev. Lett. ���, ������ (����).
[���] C. Jozwiak, C.-H. Park, K. Gotlieb, C. Hwang, D.-H. Lee, S. G. Louie, J. D.
Denlinger, C. R. Rotundu, R. J. Birgeneau, Z. Hussain, and A. Lanzara,
Nature Physics �, ��� (����).
[���] J. M. Rondinelli and N. A. Spaldin, Phys. Rev. B ��, ������ (����).
��� references