Structurally triggered metal-insulator transition in rare-earth nickelates.

Rare-earth nickelates form an intriguing series of correlated perovskite oxides. Apart from LaNiO3, they exhibit on cooling a sharp metal-insulator electronic phase transition, a concurrent structural phase transition, and a magnetic phase transition toward an unusual antiferromagnetic spin order. Appealing for various applications, full exploitation of these compounds is still hampered by the lack of global understanding of the interplay between their electronic, structural, and magnetic properties. Here we show from first-principles calculations that the metal-insulator transition of nickelates arises from the softening of an oxygen-breathing distortion, structurally triggered by oxygen-octahedra rotation motions. The origin of such a rare triggered mechanism is traced back in their electronic and magnetic properties, providing a united picture. We further develop a Landau model accounting for the metal-insulator transition evolution in terms of the rare-earth cations and rationalizing how to tune this transition by acting on oxygen rotation motions.

Introduction

First synthesized in 1971[1], rare-earth nickelates (RNiO3, with R = rare earth) are appealing for various applications[2, 3], and the possibility to tune their properties in epitaxial films and heterostructures[4] has recently fueled an even larger interest[5-7]. RNiO3 compounds belong to the family of perovskite oxides with a reference cubic structure (Fig. 1a), nevertheless not directly observed. Apart for LaNiO3, which always keeps a metallic phase and will not be further discussed here, all RNiO3 adopt at reasonably high temperature a metallic Pbnm orthorhombic phase[8]. This phase, rather ubiquitous[9] amongst perovskites with small Goldschmidt tolerance factor[8], t, corresponds to a distortion of the cubic structure arising from the appearance of combined anti-phase rotations of the oxygen octahedra along the x and y directions, R (Fig. 1b), with the same amplitude Q R and in-phase rotations of the oxygen octahedra along z, M (Fig. 1c), with amplitude Q M.

Fig. 1

Reference cubic perovskite structure and most relevant atomic distortions. a Sketch of the reference cubic perovskite cell of RNiO3 compound with R at the corner, Ni at the center, and O atoms at the middle of the faces, forming corner-shared octahedra surrounding the B atoms. Starting from this reference, three main atomic distortions drive the system successively to the Pbnm and then P21/n phases: b anti-phase rotations of oxygen octahedra of same amplitudes about x- and y-axis (R ), c in-phase rotations of oxygen octahedra about z-axis (M ), d breathing of the oxygen octahedra (B OC). The energy wells associated to the freezing of individual distortion of increasing amplitude in the cubic cell are shown for different R cations, associated to distinct tolerance factor t, and either a ferromagnetic (full line) or E′-type antiferromagnetic (dashed line, only for R = Y) spin arrangement: YNiO3 (t = 0.920, blue); GdNiO3 (t = 0.938, orange) and SmNiO3 (t = 0.947, green). The atomic distortions are normalized to their amplitude in the P21/n AFM-E′ ground state of YNiO3. Calculations are done for each compound in a cubic cell that has the same volume as the ground state

In this phase, all Ni atoms are equivalent and formally in a Jahn-Teller active d 7 (likely low spin) configuration. Surprisingly, although compatible with the Pbnm symmetry, cooperative Jahn-Teller distortions do not appear. Instead, at a temperature T MI, which strongly evolves with the R cation (i.e., with t)[10], RNiO3 compounds exhibit an electronic metal-insulator transition (MIT), which was shown to be concurrent with a structural transition from Pbnm to monoclinic P21/n symmetry[11]. This symmetry lowering is accompanied with the appearance of a breathing distortion of the oxygen octahedra, B OC (Fig. 1d), which alternatively expand and contract with amplitude Q B, according to a rock-salt pattern. This gives rise to two types of Ni sites with long and short Ni–O bonds, respectively.

At low temperature (100–200 K), RNiO3 compounds finally show a magnetic phase transition toward an unusual E′-type antiferromagnetic (AFM) spin order identified by a Bragg vector q = (1/4, 1/4, 1/4) in pseudocubic notation[11-13]. For large cations (R = Nd and Pr), the Néel temperature T N = T MI and the system goes directly from paramagnetic metal (PM-M) to AFM insulator (AFM-I). For smaller R cations, T N is much lower than T MI; the two phase transitions are decoupled and the system goes through an intermediate paramagnetic insulating phase (PM-I).

The origin of the MIT has been widely debated in the literature[14-19]. It was sometimes interpreted as a charge disproportionation at Ni sites[20]: 2(d 7) → d 8 + d 6. However, the importance of Ni 3d-O 2p hybridization—i.e., transfer of electrons from O to Ni and formation of oxygen holes keeping Ni in a d 8 configuration (i.e., )—was evoked early on[15]. Nowadays, the MIT is usually viewed as a charge ordering of the type [17, 21, 22]. In this scenario, B OC appears important to stabilize the charge ordering and open the gap. As suggested in ref. [23] and confirmed from statistical correlation techniques[24], R and M are also expected to play an active role. However, a complete picture linking electronic, structural and magnetic properties is yet to emerge.

Unlike recent theoretical studies, which were focusing specifically on the electronic properties[17–19, 21], we investigate here self-consistently the electronic and structural degrees of freedom of RNiO3 compounds from density functional theory calculations (DFT, see Methods). Specific care was given to the validation of our approach, which turns out to provide an unprecedented agreement with experimental data. Focusing on YNiO3, we show (Supplementary Notes 1–4) that not only the atomic structure but also the AFM-E′ ground state, the estimated T N and the electronic bandgap of the insulating phase are very accurately reproduced, making therefore our approach a method of choice to shed light on the interlink between these different features.

Results

Energetics of individual lattice distortions

We start from the reference cubic structure. Inspection of the phonon dispersion curves (Supplementary Note 5) reveals dominant structural instabilities at R (ω R = 310i cm−1) and M (ω M = 278i cm−1) points of the Brillouin zone (BZ), which are associated, respectively, to the R and M distortions responsible for the Pbnm phase. These imaginary frequencies ω i are linked to a negative energy curvature α i at the origin and to a typical double-well (DW) shape of the energy when freezing R and M distortions of increasing amplitude within the cubic structure, as illustrated in Fig. 1. These wells are nearly independent of the spin order but strongly evolve with the R cation: they become shallower when R size increases, consistently with a reduction of the related distortion amplitudes in the Pbnm phase.

In contrast, the B OC motion, corresponding to another phonon at R, is stable and extremely stiff (in fact the stiffest mode with ω B = 700 cm−1), in line with the single-well (SW) shape illustrated in Fig. 1. Decreasing R cation size tends to stabilize slightly further B OC, in apparent contradiction with the observed increase of T MI. As illustrated for YNiO3, switching from ferromagnetic (FM) to AFM-E′ spin order reduces slightly the curvature but does not reverse it; instead it shifts the SW to lower energy[13], yielding a finite Q B at the minimum. Although B OC tends to make the system insulating, this amplitude (25% of ground state’s value) is however not large enough to open a gap (more than 75% would be required). This shows that B OC and the magnetic order alone cannot explain the MIT by themselves.

Mode coupling and triggered mechanism

Our central result is presented in Fig. 2 where we report in panel a the evolution of the B OC energy well of YNiO3 at various fixed amplitudes of oxygen rotation motions. It highlights that although initially stable (SW), B OC is progressively destabilized (DW) by the appearance of R and M . As illustrated in panel b, α B is renormalized into , which evolves linearly with and . The slope associated to Q R is twice as large as that related to Q M, attesting that each of the three individual rotations similarly affects B OC. This behavior arises from the presence in the energy expansion of cooperative (λ < 0) bi-quadratic coupling terms between B OC and oxygen rotations (, i = M,R), which, being the lowest-order couplings allowed by symmetry, should give rise to the appearance of B OC through a triggered phase transition according to Holakovsky[25]. The same behavior is observed independently of the magnetic order (Supplementary Note 8). From now we focus on representative FM results while coming back to the role of magnetism later.

Fig. 2

Triggered mechanism. a Evolution of the energy E in terms of the amplitude of the breathing distortion (Q B) for fixed amplitudes of oxygen rotations (Q R = Q M from 0 to 110%) in the FM cubic cell of YNiO3 (same volume as the ground state). It highlights the softening of the energy well associated to B OC, triggered by the oxygen rotations R and M . b Linear evolution of the energy curvature at the origin, along Q B, in terms of the square of the amplitude of the individual distortions Q R (orange) and Q M (blue)

Landau model and phase diagram

To further assess the relevance of such a triggered mechanism in nickelates, we built a Landau model, including R , M , and B OC degrees of freedom[24], restricting ourselves to lowest coupling terms and assuming temperature dependence of the oxygen rotations only:Within this model, R and M appear at T 0R and T 0M, respectively. On cooling, they progressively develop within the Pbnm phase and renormalize the energy curvature α B of B OC as made clear when grouping the terms in Eq. (1):When reaching a critical amplitude at which , they trigger the appearance of B OC and produce concurrent structural and metal-insulator phase transitions. The phase transition appears to be second order within this model, which is however too simple to be conclusive on this point (Supplementary Note 10).

All parameters and their evolution with R were directly fitted from first principles; only Curie temperatures were uniformly scaled to reproduce the experimental T MI of YNiO3 (Supplementary Note 8). The phase diagram of nickelates as predicted within this model is reported in Fig. 3. This figure demonstrates that the cooperative coupling of B OC with R and M is a key mechanism that, by itself, accounts for the experimentally observed evolution of T MI with the tolerance factor.

Fig. 3

Nickelate phase diagram. Phase diagram of RNiO3 compounds in terms of their tolerance factor t and the temperature T. It includes three phases: a metallic Pbnm paramagnetic phase (PM-M, yellow area); an insulating P21/n paramagnetic phase (PM-I, green area); and an insulating P21/n E′-type AFM phase (AFM-I, magenta area). The yellow line and dots show the experimental evolution of T MI with the tolerance factor t. The blue line and squares is the prediction of the simple Landau model fitted on the first-principles data (FM order). The dashed blue line is the fit of the Landau expression of T MI (t) on experimental data. The blue star is the magnetic phase transition predicted for YNiO3 from first principles. The blue arrows indicate the correction to be applied on T MI for large cations when properly incorporating the change of magnetic order

Hence, the MIT in nickelates turns out to be a concrete example of triggered phase transition[26, 27], a kind of transition never identified before in simple perovskites. Indeed, although bi-quadratic interactions are generic in this class of compounds, different distortions usually compete and exclude each other[9]. The cooperative coupling of B OC with oxygen rotations pointed out here is therefore an unusual and intriguing feature, whose origin is now traced back in the electronic band structure.

Electronic origin of the triggered mechanism

In the cubic phase, as expected from the formal Ni 3d 7 occupancy, the Fermi energy E f crosses levels of dominant Ni 3d-e g character (i.e., anti-bonding Ni 3d-O 2p states); such levels form an isolated and rather dispersive set of two e g bands, shifted above the t 2g levels by the crystal field (Fig. 4). Forcing into this cubic structure a B OC distortion, associated to a phonon at q R = (1/2, 1/2, 1/2), can open a gap at q c = (1/4, 1/4, 1/4) within the e g bands but well above E f and without any direct impact on the metallic character and the occupied states. Nevertheless, the oxygen rotation motions substantially affect the e g bands (Fig. 4), reducing their bandwidth and yielding a progressive down shift of the e g levels at q c. As the rotations gradually increase and the bandwidth decreases, B OC more substantially pushes down the electronic states around E f, providing a progressive gain of electronic energy, which, in turn, results in the softening of ω B. The critical rotation amplitude at which B OC becomes unstable is precisely that at which the e g levels at q c cross E f. At these amplitudes, the electronic system itself becomes unstable; the appearance of B OC is favored and opens a gap at E f, making the system insulating. As such, the MIT can therefore be interpreted as a Peierls instability but one that is not initially present and has been triggered by oxygen rotations.

Fig. 4

Electronic properties. a Electronic dispersion curves of YNiO3 along different high-symmetry line of the Brillouin zone of the phase (FM case, majority spins): Γ = (0, 0, 0), X = (1/2, 0, 0), M = (1/2, 1/2, 0), and R = (1/2, 1/2, 1/2). The Ni 3d e g bands are highlighted in green. The Fermi energy corresponds to the horizontal blue line. The point q c = (1/4, 1/4, 1/4) is located by vertical dashed blue lines. b–g Electronic dispersion curves around the Fermi energy E f (FM order, majority spins) along high-symmetry points (coordinates in pseudocubic notations) in the Brillouin zone of the Pbnm or P21/n 20-atom cell, in which bands have been folded respect to a. Consecutive panels show the evolution of the dispersion curves when freezing into the phase increasing amplitudes of oxygen rotations (Q R = Q M = b 0%, d 50%, e 75%, or f 90%, lines) and eventually adding the breathing distortion (Q B = 100% and Q R = Q M = c 0% or g 90%, dashed lines)

Discussion

For compounds with small R cations, Q R and Q M are large and able to produce the MIT at relatively high temperatures, well above T N. For large cations (R = Nd, Pr), oxygen rotations are reduced and, from our Landau model (built on FM results), no more sufficient to trigger the MIT (Fig. 3). However, as previously mentioned, the AFM-E′ spin order is compatible by symmetry with B OC and induces its appearance as an improper order (linear shift of B OC SW, Fig. 1d). Hence, although not opening a gap in the cubic phase, the onset of the AFM-E′ order in the Pbnm phase of NdNiO3 and PrNiO3 promotes the occurrence of the MIT almost triggered by the rotations. In these compounds, we have therefore T MI = T N; the transition is more abrupt and first-order[8]. Such active role of magnetism for large cations is supported by the experimental results and discussion in ref. [28]. It is also confirmed by our first-principles calculations on NdNiO3, showing that the predicted T MI is rescaled when including the change of magnetic order: while the system prefers to stay in the Pbnm metallic phase when imposing a FM order, it switches to the P21/n phase when adopting a AFM-E′ order. The cooperative effect of the magnetic order remains true for small cations but without any impact on T MI (>T N).

In conclusion, the concurrent electronic and structural transitions at T MI in RNiO3 compounds take the form of a Peierls instability, which, primarily, is structurally triggered by the oxygen rotation motions R and M and, eventually, is further assisted by the appearance of the E′-type AFM magnetic ordering. Our Landau model, and its possible extension to incorporate explicitly strain degrees of freedom neglected here for simplicity, provides a simple and useful quantitative tool to estimate and interpret how T MI can be tuned toward the monitoring of oxygen rotation motions R and M when making solid solutions[10], applying pressure[8], or playing with the epitaxial strain and the orientation of the substrate in thin films[29]. Our findings are relevant to other families of perovskites like A2+ Fe4+ O3 compounds[30]. For instance, they can explain why CaFeO3, which exhibits oxygen rotations, undergoes a MIT while SrFeO3 and BaFeO3, which stay cubic, remain metallic. In addition, the same physics is also inherent to manganites like LaMnO3, suggesting a close competition between charge and orbital orderings in this family compounds. However, the situation is slightly different in a bismutate like BaBiO3, in which B OC is intrinsically unstable in the cubic phase[31].

Methods

First-principles calculations

First-principles calculations were performed in the framework of DFT[32, 33] using a Projected Augmented Wave (PAW) approach[34] as implemented within the ABINIT package[35-38]. The calculations relied on the Generalized Gradient Approximation using the PBEsol[39] exchange-correlation functional. We worked within a collinear spin approximation. We included a Hubbard correction U = 1.5 eV[40] on the 3d orbitals of Ni atoms. A special care was devoted to the determination of the appropriate U parameter (Supplementary Notes 1– 4).

We made use of JTH atomic potentials[41]. For the wavefunctions, we used an energy cutoff of 24 Ha (38 Ha for the second grid in the PAW spheres), which guarantees a convergence better than 1 meV on the total energy. The BZ was sampled with k-point meshes equivalent to a 12 × 12 × 12 grid in the five-atom unit cell. During structural relaxations, thresholds of 10−5 Ha/bohr on the maximum force and of 10−7 Ha/bohr3 on the maximum stress have been considered.

Structural analysis

The Goldschmidt tolerance factor[42], t = d R−O/ d Ni−O, of RNiO3 compounds were determined using Nicole Benedek’s tool[43] relying on a bond valence model[44] to calculate d R−O and d Ni−O, respectively, the ideal R–O and Ni–O bond lengths in the cubic perovskite structure.

Symmetry-adapted mode analysis have been performed with AMPLIMODE[45, 46]. The modes are normalized to their amplitude in the P21/n AFM-E′ ground state. This normalization is such that in cubic phase (volume of the ground state) Q R = 100% corresponds to rotation angles ϕ  = ϕ  = 11.33° (Ni–O–Ni angle of 157.33°), Q R = 100% corresponds to a rotation angle ϕ  = 12.12° (Ni–O–Ni angle of 155.75°) and Q B = 100% corresponds to oxygen displacements d O = 0.0358.

Landau model

The Landau model parameters have been fitted for YNiO3, GdNiO3, and SmNiO3 on first-principles data using in each case a FM cubic phase (volume of the ground state) and interpolated for the other compounds. T MI was determined analytically (Supplementary Notes 6–9).

Data availability

All relevant data are available from the authors.