Simulating Charge, Spin, and Orbital Ordering: Application to Jahn

Dec 24, 2017 - The degrees of freedom associated with orbital, spin, and charge ordering can strongly affect the properties of many crystalline solids...
2 downloads 10 Views 1MB Size
Subscriber access provided by READING UNIV

Article

Simulating charge, spin, and orbital ordering: application to Jahn-Teller distortions in layered transition-metal oxides Maxwell D Radin, and Anton Van der Ven Chem. Mater., Just Accepted Manuscript • DOI: 10.1021/acs.chemmater.7b03080 • Publication Date (Web): 24 Dec 2017 Downloaded from http://pubs.acs.org on December 24, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Chemistry of Materials is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 14 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

Simulating charge, spin, and orbital ordering: application to Jahn-Teller distortions in layered transition-metal oxides Maxwell D. Radin and Anton Van der Ven* Materials Department, University of California, Santa Barbara, Santa Barbara, California 93106, United States ABSTRACT: The degrees of freedom associated with orbital, spin, and charge ordering can strongly affect the properties of many crystalline solids, including battery materials, high-temperature superconductors, and naturally occurring minerals. This work reports on the development of a computational framework to systematically explore the ordering of electronic degrees of freedom, and presents results on orbital ordering associated with Jahn-Teller distortions in four layered oxides relevant for Li- and Na-ion batteries: LiNiO2, NaNiO2 , LiMnO2, and NaMnO2. Our calculations reveal a criterion for the stability of orbital orderings in these layered materials: each oxygen atom must participate in two short and one long transition-metal/oxygen bond. The only orderings that satisfy this stability criterion are row orderings, such as the “zigzag” ordering. The near degeneracy of such row-orderings in LiNiO2 suggests that boundaries between domains with distinct but symmetrically equivalent Jahn-Teller distortions will be relatively low in energy. Based on this result, we speculate that a microstructure consisting of nanoscale Jahn-Teller domains could be responsible for the apparent absence of a collective distortion in experiments on LiNiO2.

INTRODUCTION The properties of many crystalline materials are strongly affected by electronic degrees of freedom associated with each site, resulting in charge ordering, spin ordering, and orbital ordering. (As discussed below, orbital ordering is often intimately related to the Jahn-Teller effect.) Magnetite (Fe3O4) has played an important role historically: it was the first material found to be magnetic1 and , some two millennia later, the first material found to be charge ordered.2–4 The ground-state charge and orbital ordering, however, has only recently been determined.5 The “stripe” charge-ordering in perovskite manganates6,7 is another well-known example of simultaneous charge, spin, and orbital ordering; similar stripe ordering phenomena are thought to play an important role in cuprate high-temperature superconductors.8 Layered alkali transition-metal oxides, widely studied because of their ubiquitous application as battery electrode materials (e.g., LiCoO2),9,10 also exhibit complex charge, spin, and orbital ordering phenomena. In particular, the layered Ni11–14 and Mn15–18 oxides exhibit strong Jahn-Teller activity associated with the eg electrons of Ni3+ and Mn3+ which results in orbital ordering.19 However, a number of questions remain unanswered. For example, the nature of the ground state of LiNiO2 remains controversial, despite a large number of experimental12,20–24 and theoretical studies.19,25–27 Furthermore, the energy spectrum of excited states in the layered oxides (which deter-

mines the finite-temperature behavior) is not well understood. Unanswered questions about charge, spin, and orbital ordering motivate the use of first-principles modeling to explore electronic degrees of freedom. While the ordering of chemical species is in some sense straightforward to study computationally, the ordering of electronic degrees of freedom is less amenable to simulation. This is primarily because imposing a particular charge, spin, or orbital ordering within first-principles calculations is non-trivial: one cannot coerce electrons into a particular state. Additionally, orbital ordering is more complex than chemical ordering because different possible electronic states of an atom may be related by symmetry. (In contrast, two distinct chemical species are not related by symmetry.) A number of methods have been developed to navigate the energy landscape associated with electronic degrees of freedom in first-principles simulations. Techniques for biasing density-functional theory calculations towards particular magnetic states include the initialization of atomic magnetic moments, constraints on the total spin in the simulation cell, and energy penalties for atomic moments.28 A common strategy for imposing charge/orbital ordering is to perturb the initial atomic positions26,29,30 based on an educated guess. Charge ordering can also be induced by the application of different U parameters to different sites.31–33 Gradually increasing the U parameter (U ramping) has also been suggested as a route towards finding the global minimum of the energy landscape associated with electronic degrees of freedom.34

1

ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Arguably the most general strategy, however, is the imposition of occupation matrices (occupation matrix control)35,36, as in principle it can be used to impose spin, orbital, and charge ordering. However, a systematic methodology for exploring electronic ordering is absent. This report details a new framework that combines occupation matrix control (OMC) with classical lattice alloy theory to enable the systematic study of orderings of electronic degrees of freedom with first-principles calculations. As a case study, we apply this occupation-matrix lattice-alloy framework to orbital, spin, and charge ordering in the four layered transition-metal oxides: LiNiO2, LiMnO2, NaNiO2, and NaMnO2. In these systems, spin ordering can arise from the unpaired electrons on Ni and Mn; charge ordering can result from the disproportion of Ni3+ into Ni2+ and Ni4+; and orbital ordering occurs as a consequence of the Jahn-Teller activity of Ni3+ and Mn3+ in an octahedral environment. (The partial deintercalation of these compounds can also result in charge ordering;37,38 however the calculations presented here are restricted to the fully intercalated phases.) We have chosen to focus on LiNiO2, LiMnO2, NaNiO2, and NaMnO2 because of the extensive interest in them from both the solid-state physics community17,18,23,24,39 and the battery community.9,10,40–42 Much of the interest from the physics community stems from the potential for these materials to realize exotic ground states associated with magnetic frustration and orbital degeneracy. For the battery community, layered LiNiO2 and LiMnO2 represent end members of the LiNixMnyCo1−x−yO2 (NMC) and LiNixCoyAl1−x−y O2 (NCA) families, which are extensively used for intercalation electrodes in commercial Li-ion batteries.9,10 The first main result is a criterion for the stability of orbital orderings associated with Jahn-Teller distortions in the Ni and Mn layered oxides: each oxygen must participate in two short M-O bonds, and one long M-O bond. Configurations which did not satisfy this criterion were found to be locally unstable: they would spontaneously relax to different orderings in the absence of any constraints. The second main result is the large number of nearly degenerate row orderings of orbital occupancy in LiNiO2. This suggests that the twin boundary energy is quite low, and that the apparent absence of a collective Jahn-Teller distortion in LiNiO2 might be due to a highly twinned nanostructure. The high stability of such boundaries appears to be related to interactions between the Li+ ions and occupied eg* orbitals.

Page 2 of 14

A further subtlety is that the ordering of cations breaks the cubic symmetry of the rocksalt host, and as a result the octahedra are slightly compressed or stretched along the (111) rocksalt axis. This reduces the point group symmetry of the cation sites from Oh to D3d, resulting in a further splitting of the t2g states into states with a1g and eg symmetry. (The two eg states of the Oh point group, on the other hand, remain degenerate when the octahedra are compressed.) Because the orbital degeneracy of Ni3+ and Mn3+ does not directly involve t2g states, their splitting does not appear to play a major role in the compounds under consideration presently. (However, the splitting of t2g states does play an important role in layered Ti43,44 and V45 oxides, where the t2g manifold is only partly occupied.) The presence of a single electron in the doubly degenerate eg level leads to a local Jahn-Teller distortion in the Ni and Mn layered oxides.19 In essence, the distortion causes one of the eg levels to go down in energy and the other to go up; this leads to an overall lowering of the energy because only one of the eg states is occupied. The distortion observed in layered Ni and Mn oxides is a “positive” Jahn-Teller distortion: two metal-oxygen bonds lengthen while the other four contract. There are three symmetrically equivalent positive distortions, as shown in Figure 2. Figure 2 also shows the shapes of the eigenstates in the distorted octahedron that are derived from the eg orbitals. This picture is complicated by the hybridization of transition-metal eg orbitals with oxygen 2p orbitals, which results in the formation of bonding and antibonding states with eg symmetry. In the Ni and Mn oxides, the bonding states are entirely filled and the Jahn-Teller effect arises from the antibonding states (often referred to as eg*).46 Antibonding states with lobes pointed along lengthened bonds will be lower in energy than those with lobes pointed along contracted bonds, as indicated in Figure 2.

CRYSTALLOGRAPHY The structures of LiNiO2, LiMnO2, NaNiO2, and NaMnO2 are all derived from the α-NaFeO2 structure (shown in Figure 1), also referred to as O3. This structure is a cation-ordered rocksalt, with the alkali and transition metals residing in alternating (111) planes. All of the cations are coordinated octahedrally by oxygen. In a perfectly octahedral environment, the five transition-metal d levels are split into two degenerate eg states whose lobes point towards the anions, and three degenerate t2g states whose lobes point between anions. 2

Figure 1. Crystal structure of a layered oxide AMO2 with the O3 (α-NaFeO2) structure.

A contraction of two bonds and lengthening of four (not shown) is referred to as a negative distortion. Such

ACS Paragon Plus Environment

Page 3 of 14 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

negative distortions have been reported, for example, in the case of degeneracies in the t2g states in layered NaVO2.45 However, since degeneracies in the eg states generally lead to a positive distortion in octahedral environments,19,47–49in this study we do not consider negative distortions.

Figure 2. The three symmetrically equivalent positive JahnTeller distortions of an octahedral complex.

The six orbitals in Figure 2 are not linearly independent: they span a vector space of dimension two. To illustrate this graphically, Figure 3 shows how any of these six ψ n (r) orbitals can be expressed as

ψ n ( r ) = aψ 1 ( r ) + bψ 2 ( r ) , where ψ1 ( r ) and ψ 2 ( r ) are the two orbitals associated with distortion 1. The arrows indicate the basis formed by ψ1 ( r ) and ψ 2 ( r ) . The points in Figure 3 that correspond to the six ψ n ( r ) of Figure 2 are marked with black dots; the six points corresponding to −ψ n ( r ) are also marked.

3

Figure 3. Illustration of how any eg orbital can be constructed as a linear combination of two reference eg orbitals. The circle represents orbitals that have a norm of unity.

Figure 2 and Figure 3 illustrate the relationship between the Jahn-Teller distortions and orbital ordering. The term distortion refers to a displacement of atoms, whereas orbital ordering refers to the occupation of different eg orbitals at each site. One describes the state of the ions while the other describes the state of the electrons. When the electrons reside in their ground state, there is a oneto-one mapping between distortions and orbital occupancy, as illustrated in Figure 2. In this work, we assume that electrons reside in the ground state with regards to orbital occupancy; this is justified when the energy splitting induced by the distortions is significantly larger than the thermal energy. Figure 4 shows several important orbital orderings (or equivalently, orderings of Jahn-Teller distortions) in the layered crystal structure. The collinear ordering, shown for a single MO2 layer in Figure 4a, has the same orbital occupancy at each transition-metal site; this corresponds to the local environment of each site being distorted in the same manner. The simplest non-collinear ordering within a layer is the zigzag (i.e., herringbone) ordering, which consists of alternating rows of two different orientations as shown in Figure 4b. Note that there are actually two zigzag orderings which exhibit the in-plane ordering of Figure 4b; these two differ in the relative alignment of Jahn-Teller orderings from one layer to the next. Another important ordering is the trimer ordering, Figure 4c.

ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 14

Figure 4. Examples of Jahn-Teller orderings within the layers of a layered oxide. Only the lengthened M-O bonds are shown. The dash black lines indicates the unit cell of each ordering, while the short solid black lines indicate the orientation of the JahnTeller distortion at each transition-metal site.

tent any electrons with spin σ possess both φi and φ j

COMPUTATIONAL METHODS The orbital, charge, and spin orderings are considered here in the context of lattice alloy theory.50–53 Within this model, each possible orbital, charge, or spin state is considered as a discrete species. For example, six species that represent distorted high-spin Mn3+ in an octahedral environment arise from the three orientations for a positive Jahn-Teller distortion and the two possible spin states, s = 2 and s = −2. (These spin states arise from the four unpaired electrons each with spin 1 2 .) The use of such a methodology requires a mapping between first-principles simulations and the lattice model: we must be able to impose (or at least bias) first-principles calculations towards an arbitrary ordering of electronic states, and also must be able to determine if, after relaxation, the configuration has retained this state or relaxed into a different one. In the present work, we employ occupation matrices as the foundation for such a mapping. Occupation matrices characterize the electronic state of each site in terms of the single-particle orbitals obtained from electronic structure calculations. For a given orbital angular momentum channel l , the spin-up and spindown occupation matrices are each described by a ( 2l + 1) × ( 2l + 1) matrix defined by54 (1) where

nijσ = ∑ f bσk ψ bσk φi φ j ψ bσk b, k

φi

and

φj

are reference atomic orbitals, σ

σ = ± 1 is the spin channel, ψ bk is the Kohn-Sham wavefunction of band b with spin σ at a particular k point,

fbσk is the occupation (between zero and one) of that state. The value of nijσ in essence measures to what ex4

character. We adopt the following notation: while n is a ( 2l + 1) × ( 2l + 1) × 2 array (i.e., a pair of matrices), nσ will refer to the ( 2l + 1) × ( 2l + 1) matrix corresponding to spin channel σ . Within the lattice alloy picture, each possible “species” (i.e., electronic state of an atom) is defined by a pair of occupation matrices describing the spin-up and spindown electron configuration. We will denote the pair of occupation matrices for a species p as n p . The application of a symmetry operation Sˆ that is respected by the host structure can transform one species p into a different ˆ p = n p ′ . In other words, unlike a convenspecies p′: Sn tional lattice alloy model, the species considered here can be anisotropic. For example, consider the effect of a threefold rotation operation on the six different Jahn-Teller orientations and spin states of Mn3+ described above: this rotation will cyclically permute the three spin-up states and also cyclically permute the three spin-down states (cf. Figure 2). The anisotropy of the species fundamentally changes the enumeration of inequivalent orderings: configurations that are equivalent in an alloy of isotropic components may be inequivalent in an alloy of anisotropic components, and orderings that are equivalent in the alloy of anisotropic components may be inequivalent in an alloy of isotropic components. For example, consider the volume 1 configurations obtained by decorating the αNaFeO2 structures with the six states of Mn3+. There are six such configurations, but all are symmetrically equivalent; therefore only one needs to be calculated. On the other hand, if six unique isotropic species (e.g., different chemical elements) are permitted to occupy the transi-

ACS Paragon Plus Environment

Page 5 of 14 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

tion-metal site in the α-NaFeO2 structure, then there are six unique configurations with the periodicity of the volume 1 primitive cell. The enumeration of symmetrically unique orderings can be accomplished once it is determined how the symmetry operations permute the species. The action of symmetry operations on the species is described by the ˆ p = n p ′ and permutation representation: P Sˆ = 1 if Sn pp ′

( )

( )

( )

ˆ p ≠ n p ′ . P Sˆ is a M × M matrix, where M Ppp ′ Sˆ = 0 if Sn is the number of species. For example, a threefold rotation acting on the six states of Mn3+ would be represented by:

0  1 0 P Cˆ 3 =  0 0  0

0 0 0  0 0 0 0 0 0  0 0 1 1 0 0  0 0 0 1 0  0 1

0 1 0 0

( )

(2)

0 0 0 0

where the first three species are the spin-up states and the second three the spin-down states. The permutation representation can be determined from the action of the symmetry operations on the elements of the occupation matrices. As discussed in the Supporting Information, the occupation matrices transform according to

( ) ˆ Sn

(3)

σ

( )

( )

= D Sˆ −1 nt Sˆ σ D Sˆ ( )

( )

where D Sˆ is the ( 2l + 1) × ( 2l + 1) representation that describes how the reference atomic orbitals are transformed by symmetry operations, and t Sˆ = ±1 indicates

()

whether or not a symmetry operation includes time reversal (which has the effect of swapping spin-up and spindown states). Algorithms for the determination of the permutation representation from the action of symmetry on n and use of the permutation representation in enumerating unique orderings will be described in a forthcoming publication.55 In order to check whether electronic degrees of freedom changed during relaxation, we define a distance metric between occupation matrices:

d ( n, n′ ) = (4)

∑σ ( n σ − n′σ ) ij

i, j ,

ij

2

 2 = tr ∑ ( nσ − nσ′ )  σ 

Although other distance metrics can be defined, the metric employed here is physically motivated because it is invariant with respect to rotations and reflections (see Supporting Information). We assign the occupation of each site after relaxation based on which reference occupation matrix is closest to the relaxed occupation matrix. In the calculations presented in this study, we considered three types of lattice alloy models. The first case represents different orbital orderings of M3+ ions with ferromagnetic spin ordering; this is similar to a ternary alloy 5

with the three species representing the three symmetrically equivalent orientations of the positive Jahn-Teller distortion. The second case allows for different magnetic orderings, and has six species arising from the three orbital orientations and two spin states. (We presently assume the spin to be collinear.) The third case considered was a disproportionated model, corresponding to a binary alloy between M2+ and M4+. The symmetrically inequivalent orderings were enumerated with the Clusters Approach to Statistical Mechanics (CASM) code.56–59 This enumeration accounted for the permutation of species upon application of symmetry operations to eliminate symmetrically equivalent orderings. Time-reversal symmetry was not explicitly included in this process, but enumeration was restricted to orderings where the number of spin-down sites was no greater than the number of spin-up sites. Once enumerated, the electronic states were imposed in DFT calculations using occupation-matrix control (OMC), as implemented by Allen and Watson.35,36 OMC adds a non-local bias potential (of the same form as the U on-site Coulomb correction) that invites the electrons to relax into a target electronic configuration, specified in terms of an occupation matrix. The target occupation matrices used for the different Jahn-Teller orientations in each compound were taken to be the occupation matrices obtained by fully relaxing the three variants of the collinear Jahn-Teller ordering of that compound. Because energies obtained with the OMC bias potential on are not physically meaningful, a second DFT calculation was performed after relaxing with the bias potential. We tried both performing full relaxations starting from the geometries and wavefunctions obtained with OMC, and performing static calculations keeping the OMC geometry fixed. For structures that retained the same electronic ordering upon full relaxation, we found that the single-point and fully relaxed energies were nearly identical, as shown in Figures S1 and S2 of the Supporting Information. DFT calculations were performed with the Vienna ab initio Simulation Package (VASP)60–63 with projector augmented-wave (PAW) pseudopotentials.64 All calculations were spin polarized, employed plane-wave basis sets with 600 eV energy cutoffs, and sampled the Brillouin zone with a k-point mesh density of 26 Å or more. Relaxations used Gaussian smearing of width 0.1 eV and were converged to within a force convergence criteria of 0.02 eV/Å; these were followed by single-point calculations using the tetrahedron method with Blöchl corrections65 in order to obtain accurate energies. Except where otherwise stated, calculations employed the PBE GGA exchange-correlation functional66 with the rotationally-invariant on-site Coulomb corrections of Dudarev et al.67 A U − J value of 5 eV was used, representing typical values applied to Mn and Ni in oxides.67–69 For comparison, the energies of Jahn-Teller orderings in LiNiO2 with a smaller U − J of 2 eV are provided in the Supporting Information. The smaller value of U − J has little effect on the relative energies of different orderings,

ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 14

although it does make the distorted structures somewhat less stable compared to the undistorted geometry; this result is in agreement with prior calculations examining the effect of U on Jahn-Teller distortions in LiNiO2.25 A few calculations were performed with the HSE hybrid functional70,71 as spot checks, as discussed below. These employed a slightly coarser k-space grid density of 18 Å and lower plane-wave cutoff of 530 eV.

RESULTS AND DISCUSSION We now apply this framework to study spin and JahnTeller orderings in two nickel oxides (LiNiO2 and NaNiO2) and two manganese oxides (LiMnO2 and NaMnO2). For the nickel oxides, we consider different orderings of Jahn-Teller orientations, but restrict the magnetic ordering to ferromagnetic because the intra-plane interactions are ferromagnetic in LiNiO272,73 and NaNiO2.13,74–76 On the other hand, LiMnO219,77–79 and NaMnO217,18,39,80,81 possess antiferromagnetic in-plane interactions and so it is necessary to include spin ordering in order to capture the ground state. LiNiO2 Figure 5 shows the single-point energies for structures that were relaxed with OMC. The colors indicate whether these structures are locally stable when the OMC bias potential is removed: blue indicates that the structure retained its Jahn-Teller ordering, while orange indicates that it relaxed to a different ordering. (For the orderings that are stable, the single-point and relaxed energies are essentially identical; see Figure S1 in the Supporting Information.) Figure 5 shows that there are many competitive low energy states that involve zigzag-type orderings, with rows of Ni being distorted in the same direction. Many of the enumerated orderings have the same in-plane ordering but differ in how the in-plane ordering is stacked along the c-axis of the layered crystal. DFT+U predicts the lowest energy ordering to be a zigzag ordering that repeats every two rows; after full relaxation, this is 12 meV/LiNiO2 lower in energy than the collinear ordering, consistent with prior DFT+U studies.25,26 Other zig-zag-type orderings involve repeat units larger than two. For example, the second lowest energy structure found for LiNiO2 repeats every three rows.

Figure 5. Energies of Jahn-Teller orderings in LiNiO2 relative to the undistorted state. The energies represent the singlepoint energy at the geometry obtained with occupationmatrix control (OMC). Blue marks represent structures that retained the ordering when further relaxed without OMC, whereas orange represents structures that spontaneously relaxed into different orderings

The high stability of the zigzag orderings relative to the collinear ordering appears to be connected to interactions between eg* states and lithium ions. This Ni-O-Li interaction has previously been observed to play an important role in the relative energies of lithium/vacancy orderings in LixNiO2,37 and the site-projected density of states (not shown) confirms that the eg* states hybridize to some degree with Li s and p states. The first point to consider is how Li 2s and 2p orbitals can couple occupied eg* states in different layers. In the collinear structure, the two eg* states that bond with a Li site are 180° apart (cis configuration); as shown in Figure 6, this can lead to a favorable bonding interaction between the oxygen and lithium orbitals. A similar 180° interaction is present in the ground state zigzag structure. There is a higher energy structure which has the same in-plane zigzag ordering, but differs in the stacking of the layers. Because this structure has a 90° Ni-O-Li-O-Ni interaction (trans configuration), we expect bonding between layers to be less favorable. (We shall refer to these two structures as cis- and trans-zigzag orderings.) The importance of Ni-O-Li bonding is supported by the fact that the trans-zigzag structure is higher in energy than the collinear and cis-zigzag orderings, as shown in Table 1. Further geometric considerations suggest that the ciszigzag ordering may have more favorable Ni-O-Li interactions than the collinear ordering. Table 1 shows the lengths of Li-O bonds in the collinear, cis-zigzag, and trans-zigzag structures. In the collinear case, the Li-O bonds aligned with the occupied eg* orbitals are slightly longer than those aligned with the unoccupied eg* orbitals (2.139 Å vs. 2.102 Å). That the bonds aligned with the occupied orbitals are longer can be attributed to the strain in the plane of the layers arising from the collective

6

ACS Paragon Plus Environment

Page 7 of 14 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

Jahn-Teller distortion. In contrast, in the zig-zag orderings the Li-O bonds aligned with the occupied eg* states are significantly shorter than those aligned with unoccupied eg* states. We speculate that the shorter Li-O distances along occupied eg* orbitals may contribute to the stability of the cis-zigzag structure, as these shorter bonds will enhance Ni-O-Li bonding.

quite close in energy to the lowest energy Jahn-Teller distorted phase. After full relaxation, the lowest energy disproportionated structure is 7 meV/LiNiO2 higher in energy than the lowest energy Jahn-Teller distorted structure.

Table 1. Bond lengths in selected LiNiO2 Jahn-Teller orderings. Collinear

Zigzag (cis)

Zigzag (trans)

Energy (meV/LiNiO2)

−108

−120

−103

Ni-O distance (Å)

2.137, 1.892

2.116, 1.900

2.105, 1.907

Li-O distance, (Å)

Along occupied eg*

2.139

2.035

2.033

Along unoccupied eg*

2.102

2.096, 2.238

2.101, 2.243

2+

4+

Figure 7. Energies of disproportionated Ni /Ni orderings in LiNiO2 relative to the undistorted state having a uniform 3+ distribution of Ni . The energies represent the single-point energy at the geometry obtained with occupation-matrix control (OMC). Blue marks represent structures that retained the ordering when further relaxed without OMC, whereas orange represents structures that spontaneously relaxed into different orderings.

Figure 6. Illustration of interaction between Li and eg* orbitals in the collinear Jahn-Teller ordering. The dotted lines indicate the Li-O bonds which desire to contract because of interactions between the Li site and the electron in the eg* state.

For the collinear and cis-zigzag orbital orderings, we additionally considered several antiferromagnetic spin orderings (seven in the collinear structure and nine in the cis-zigzag structure). For a given orbital ordering, the antiferromagnetic spin orderings were within 3 meV/LiNiO2 of the ferromagnetic spin ordering – an energy difference that is comparable to the numerical errors in the calculation. None of the antiferromagnetic configurations was found to be lower in energy than the ferromagnetic cis-zigzag configuration.

To further resolve the relative energies of competing groundstates in LiNiO2, we calculated the energies of selected structures with the HSE functional. 70,71 Like PBE+U, HSE predicts the collinear and zigzag Jahn-Teller orderings to be nearly degenerate. However, while PBE+U predicts the disproportionated cell to be nearly degenerate with these orderings, HSE predicts the disproportionated structure to be nearly 70 meV higher in energy than the Jahn-Teller structures. That the disproportionated structure is less stable than the Jahn-Teller distorted structures is supported by neutron and EXAFS experiments, which show that there are more short Ni-O bonds than long ones.11,12 (A positive Jahn-Teller distortion results in two long bonds and four short, while disproportion results in an equal number of short and long.) However, thermal excitations may nevertheless induce the formation of a substantial concentration of Ni2+ and Ni4+. Furthermore, crystallographic defects such as Li vacancies or Ni in the Li layer may induce the formation of Ni2+ and Ni4+.82,83

We also explored the previously proposed hypothesis that LiNiO2 could disproportionate, resulting in a mixture of Ni2+ and Ni4+.25 Figure 7 shows the energies calculated for all Ni2+/Ni4+ orderings with unit cells containing up to four formula units. Our findings confirm the results of Chen et al.:25 a zigzag arrangement of Ni2+ and Ni4+ is 7

ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Table 2. Relative energies of selected LiNiO2 structures. Relative energy (meV/LiNiO2) PBE+U, U = 5 eV

HSE06

0

0

Zigzag Jahn-Teller

−11.9

−0.2

Disproportionated

−5.2

68.9

Collinear Teller

Jahn-

NaNiO2 DFT calculations indicate that the ordering of JahnTeller distortions in NaNiO2 is rather different from LiNiO2. As shown in Figure 8, the collinear distortion is predicted to be the ground state of NaNiO2, in agreement with experimental observation.13 The zigzag orderings are substantially higher in energy than the collinear ordering, and unlike LiNiO2, the trans-zigzag is preferred over the cis-zigzag by 11 meV/NaNiO2. This suggests that the bonding between sodium and occupied eg* states is weaker in NaNiO2 than LiNiO2, although other factors such as intercalant size could play a role as well.84

Page 8 of 14

the experimentally observed crystal structure.15,16 The lowest energy ordering is predicted to be antiferromagnetic, consistent with the experimental observation of antiferromagnetism in LiMnO2.79 In the ground state, the two short Mn-Mn bonds are antiferromagnetic while half of the four long Mn-Mn bonds are antiferromagnetic and half are ferromagnetic; this is in agreement with previous first-principles calculations.19,77,78 Compared to the LiNiO2 and NaNiO2, the strength of the Jahn-Teller distortion is much larger in LiMnO2 in the sense that the difference in energy between the ground state and undistorted state is much larger. The lowest energy zigzag configuration in LiMnO2 (excluding zig-zig-zag orderings) is 59 meV/LiMnO2 higher in energy than the ground state. It has an antiferromagnetic spin ordering, and Mn-O-Li bonds in the cis configuration.

Figure 9. Energies of Jahn-Teller orderings in LiMnO2 relative to the ferromagnetic undistorted state. The energies represent the single-point energy at the geometry obtained with occupation-matrix control (OMC). Blue marks represent structures that retained the ordering when further relaxed without OMC, whereas orange represents structures that spontaneously relaxed into different orderings. Figure 8. Energies of Jahn-Teller orderings in NaNiO2 relative to the undistorted state. The energies represent the singlepoint energy at the geometry obtained with occupationmatrix control (OMC). Blue marks represent structures that retained the ordering when further relaxed without OMC, whereas orange represents structures that spontaneously relaxed into different orderings.

LiMnO2 Figure 9 shows the relative energies of orbital and magnetic orderings of LiMnO2 versus the average magnetic moment. (An average magnetic moment of zero corresponds to an antiferromagnetic spin ordering, while a magnetization of 4 μB/LiMnO2 corresponds to a ferromagnetic ordering.) The calculations predict the collinear Jahn-Teller distortion to be the most stable orbital ordering for any given net magnetization, in agreement with 8

NaMnO2 Figure 10 shows the energies of spin and orbital orderings in NaMnO2. Like LiMnO2, NaMnO2 favors a collinear orbital ordering with antiferromagnetic in-plane magnetic ordering, in agreement with prior experiments17,18,39,85,86 and calculations.80,81 This magneto-elastic coupling breaks the equivalence between two nearest-neighbor Mn-Mn distances. In the calculated ground state structure, these two distances differ by 0.17%, a value similar to that obtained by low-temperature neutron diffraction (0.23%).17 The higher energy configurations of NaMnO2, however, differ from those of LiMnO2. The lowest energy zigzag configuration in NaMnO2 (excluding zig-zig-zag orderings) is 103 meV/NaMnO2 higher in energy than the ground state. Unlike the lowest energy zigzag structure in LiMnO2, this structure has a ferromagnetic spin ordering and Mn-O-Na bonds in the trans configuration.

ACS Paragon Plus Environment

Page 9 of 14 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials Table 3. Energies of Jahn-Teller and magnetic orderings, relative to the ferromagnetic undistorted (rhombohedral) state. Material

Figure 10. Energies of Jahn-Teller orderings in NaMnO2 relative to the ferromagnetic undistorted state. The energies represent the single-point energy at the geometry obtained with occupation-matrix control (OMC). Blue marks represent structures that retained the ordering when further relaxed without OMC, whereas orange represents structures that spontaneously relaxed into different orderings.

GENERAL DISCUSSION We have introduced an approach that enables a systematic enumeration of orbital, charge and magnetic orderings in crystalline solids. This has enabled us to perform an in-depth study of the interactions between different localized electronic degrees of freedom in layered transition metal oxides. Applied to LiNiO2, NaNiO2, LiMnO2 and NaMnO2, we have been able to identify several general trends. First, the strength of the Jahn-Teller distortion associated with orbital ordering varies substantially with chemistry. As shown in Table 3, the strength of a Jahn-Teller distortion is much stronger in terms of energy in Mn oxides vs. Ni oxides. Second, orderings where an oxygen site does not participate in exactly two short M-O bonds and one long M-O bond are mechanically unstable. This trend is consistent with a bond valence sum picture: too many long bonds makes an oxygen anion undercoordinated; too many short bonds makes it overcoordinated. As shown in the Supporting Information, the only orderings satisfying this oxygen bond criterion are row orderings consisting of any combination of “zigs” and “zags”, as illustrated in Figure 11.

Figure 11. Examples of structures that satisfy the twoshort/one-long criterion for M-O bonds. The lines show the Jahn-Teller axis at each transition-metal site within a layer.

9

Ground state JahnTeller Ordering

Ground state inplane magnetic ordering

Distortion energy (meV/ AMO2)

LiNiO2

Zigzag

FM

−120

−12

NaNiO2

Collinear

FM

−131

22

LiMnO2

Collinear

AFM

−354

59

NaMnO2

Collinear

AFM

−350

103

E[zigzag] – E[collinea r] (meV/ AMO2)

Although DFT predicts collective Jahn-Teller distortions to be thermodynamically favorable in LiNiO2, no collective distortion (zigzag or collinear) has been experimentally observed in LiNiO2, even at temperatures as low as 1.4 K.87,88 While thermal fluctuations resulting in a dynamic Jahn-Teller distortion can explain the absence of a collective distortion at room temperature, the apparent absence of a collective distortion at low temperatures is a long-standing mystery. (In contrast, NaNiO2, LiMnO2, and NaMnO2 are experimentally observed to exhibit collective Jahn-Teller distortions, even at room temperature.13,15–17,85) Even though no collective distortion occurs in LiNiO2, the observation of two distinct Ni-O bond lengths with EXAFS20 and neutron diffraction12 experiments indicates that local Jahn-Teller distortions occur. Several possible explanations for the absence of a collective distortion at low temperatures have been hypothesized: 1.

The disruption of long-range orbital ordering by point defects, especially Ni atoms substituting on the Li layer. This defect chemistry results in a composition of Li1–zNi1+zO2, with z typically being on the order of several percent.82,83

2.

The disproportionation of Ni3+ into Ni2+ and Ni4+, rather than JT activity.25

3.

Trimer Jahn-Teller orderings that, at least locally, preserve three-fold symmetry.12,22

The first possibility (point defects) lies outside the scope of this work and is not considered further. The second possibility (disproportionation) seems unlikely to explain the absence of a collective distortion: as discussed above, DFT calculations suggest that Jahn-Teller distortions are more favorable than disproportionation. Furthermore, experimental data shows that there are more short Ni-O bonds than long ones;12,20 this is not consistent with full disproportionation, which would result in equal numbers of short and long bonds. The third possibility, trimer orderings, can be ruled out based on the results presented here. Figure 5 shows that the trimer model is locally unstable, and furthermore has the highest single-point energy of all orderings considered. The instability of the trimer ordering is consistent with prior DFT calculations.25,26 The trimer orderings represent the most egregious violation of the criterion for MO bonds: instead of having one long and two short bonds,

ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

each oxygen in the trimer ordering has either three long bonds or three short bonds. It was previously suggested such trimers may be more stable when arranged into anti-phase domains.12 To test this trimer-domain model, we constructed a large supercell model containing 12 formula units as shown in Figure S6 of the Supporting Information. While relaxation with OMC did result in the formation of anti-phase trimer domains, this ordering was found to be unstable when the OMC bias potential was removed. This result indicates that LiNiO2 will not form a trimer-domain structure. Based on the calculated energies for Jahn-Teller orderings, we hypothesize another possible explanation for the apparent absence of a collective Jahn-Teller distortion in LiNiO2: a nanostructure with Jahn-Teller domains too small to be readily observed by diffraction. Figure 12a illustrates a hypothetical boundary between two zigzag domains. The near degeneracy of different zig-zag-type orderings suggests that the energy to create such boundaries in LiNiO2 is quite low. Similar domain formation could also occur if a collinear distortion, rather than zigzag ordering, were the ground state; Figure 12b illustrates such a boundary between two collinear domains. (Note that the interface in Figure 12a is an antiphase boundary because the two domains are related by a translation, while the interface in Figure 12b is a twin boundary because the two domains are related only through rotation/reflection.) In addition to explaining the absence of an observable collective distortion, the formation of such Jahn-Teller domains could explain the anomalous behavior observed in neutron PDF wherein long-range ordering becomes weaker as the temperature is lowered.12 To further explore the structure of such domain boundaries, explicit models between zigzag domains in LiNiO2 and collinear domains in NaNiO2 were calculated, as discussed in the Supporting Information. In these models, each domain consisted of four rows. We considered several different boundary orientations, and found the lowest energy domain superstructure in LiNiO2 to be 10 meV/supercell higher in energy than the LiNiO2 ground state (a single cis-zigzag domain). In contrast, the lowest energy domain superstructure in NaNiO2 was 65 meV/supercell higher in energy than the NaNiO2 ground state (a single collinear domain). That the domain boundary energy is significantly lower in LiNiO2 than NaNiO2 supports the hypothesis that the absence of a collective Jahn-Teller distortion in LiNiO2 is due to a Jahn-Teller domain nanostructure Recently, an analogous phenomenon has been suggested to occur in NaMnO2: the formation of antiferromagnetic domains within collinearly distorted crystallites.86 Once the direction of the collinear Jahn-Teller distortion is fixed, there are four symmetrically equivalent ground state magnetic orderings; these arise from the two equivalent directions in which the ferromagnetic chains can run (resulting in different magnetoelastic distortions), and then two translational variants for each of these (i.e., variants where the spins are reversed). The microstructure of NaMnO2 was previously suggested to consist of domains of translational variants,86 whose interfaces represent an10

Page 10 of 14

tiphase boundaries as shown in Figure 12c. We speculate that domains with differing magnetoelastic distortions may also coexist. The interfaces between such domains represent twin boundaries, as shown in Figure 12d. To test the feasibility of antiphase- and twin-boundary formation in NaMnO2, we calculated the energy of a zigzag spin ordering superimposed on a collinear Jahn-Teller distortion, as shown in Figure S7 of the Supporting Information. This structure is within 1 meV/NaMnO2 of the ground state structure (the row ordering of spins), indicating that the energy to create the boundaries shown in Figure 12c and Figure 12d is very small. The near degeneracy of the zigzag and row orderings of spins can be rationalized by noting that in both structures, each Mn exhibits 2 antiferromagnetic short Mn-Mn bonds, 2 antiferromagnetic long Mn-Mn bonds, and 2 ferromagnetic long Mn-Mn bonds.

Figure 12. Hypothetical boundaries between Jahn-Teller distorted domains. (a) An antiphase boundary between zig-zag domains. (b) A twin boundary between collinear domains. (c) An antiphase boundary between antiferromagnetic domains that have the same collinear Jahn-Teller distortion. (d) A twin boundary between antiferromagnetic domains that have the same collinear Jahn-Teller distortion but different variants of antiferromagnetic order. In (c) and (d), the dark and light lines represent spin-up and spin-down transitionmetal sites.

A final point of consideration is the nature of topological defects in the ordering of Jahn-Teller distortions, which can manifest as junctions between three or more domains. Figure 13 illustrates two hypothetical examples of topological defects. The left panel shows a junction between three collinear domains. In this structure, every

ACS Paragon Plus Environment

Page 11 of 14 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Chemistry of Materials

oxygen atom participates in two long bonds and one short bond, except for the oxygen site at the center of the junction; this site has three short bonds. The right panel shows a junction between six collinear domains. In this structure, every oxygen retains two short and one long MO bond, but the transition-metal ion at the center of the junction (represented by a black dot) possesses no long M-O bonds.

transformation from a high symmetry phase to the low temperature, low symmetry phase. We conclude by discussing the applicability of the occupation-matrix lattice-alloy framework presented here to other problems. The results above show that the energy landscape of the layered oxides has many local minima that were not previously identified (e.g., the trans-zigzag ordering); that the simulations were able to find these minima suggests that the occupation-matrix lattice-alloy framework would be effective in exploring electronic degrees of freedom in other materials. Another natural extension of this work would be to construct a lattice-alloy cluster expansion to predict the energies of charge/spin/orbital orderings.

ASSOCIATED CONTENT Figure 13. Hypothetical triple and sextuple junctions between collinear Jahn-Teller domains in a layered oxide. The black dot in the right panel represents a transition-metal site with no Jahn-Teller distortion.

We speculate that topological defects such as those shown in Figure 13 would tend to form around chemical point defects. For example, the junction shown in the right panel of Figure 13 might be centered on a transitionmetal ion that is not Jahn-Teller active. In the case of LiNiO2, this could be a Ni2+ which serves to charge compensate Ni defects in the Li layer,82,83 or a Ni4+ which serves to charge compensate Li vacancies. In the case of an alloy, other chemical species may reside at the center of topological defects. For example, previous DFT calculations on Li(NixCoyAlz)O2 found the ordering of Co and Al to be coupled to the orientation of Jahn-Teller distortions on nearby Ni sites.89

CONCLUSION A systematic investigation of charge, spin, and orbital orderings in the layered oxides reveals two criteria for the stability of Jahn-Teller orderings. First, each transition metal must exhibit two long bonds and four short, due to the well-known preference for positive Jahn-Teller distortions in octahedral environments.19,47–49 The second criterion is that each oxygen will exhibit one long bond and two short. These two constraints do not uniquely determine the orbital ordering: there is an infinite number of zig-zag-like row orderings that satisfy both criteria. Whether the ground state will be a zigzag ordering or the collinear ordering will depend on the details of the system. Importantly, in LiNiO2 the zigzag-type orderings are very close in energy to the collinear ordering, possibly due to strong interactions between eg* orbitals and Li+ ions. The near degeneracy of LiNiO2 orbital orderings suggests that LiNiO2 crystallites may exhibit a complex microstructure with nanoscale domains. Such a nanostructure could explain why a collective Jahn-Teller distortion is not observed by diffraction in LiNiO2. Similar phenomena have been widely studied in the context of martensitic transformations in alloys,90 where twinned microstructures emerge to accommodate the strain accompanying the

Supporting Information. Derivation of symmetry properties of occupation matrices, comparison of single-point and relaxed energies, results obtained with different values of U, topological analysis of Jahn-Teller orderings. This material is available free of charge via the Internet at http://pubs.acs.org.

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected]

Notes The authors declare no competing financial interest.

ACKNOWLEDGMENTS M. Radin would like to thank John C. Thomas for his assistance and Ram Seshadri for valuable discussion. This work was supported as part of the NorthEast Center for Chemical Energy Storage (NECCES), an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Basic Energy Sciences under Award # DE-SC0012583. Computational resource support was provided by the Center for Scientific Computing at the CNSI and MRL: an NSF MRSEC (DMR-1121053) and NSF CNS-0960316. Computational resources provided by the National Energy Research Scientific Computing Center (NERSC), supported by the Office of Science and U.S. Department of Energy, under Contract DE-AC02-05CH11231, are gratefully acknowledged.

REFERENCES (1)

(2)

(3)

(4) (5)

(6)

(7)

11

Lacheisserie, E. D. T. de; Gignoux, D.; Schlenker, M. Magnetism: Fundamentals; Springer Science+Business Media: New York, 2005. Verwey, E. J. W. Electronic Conduction of Magnetite (Fe3O4) and Its Transition Point at Low Temperatures. Nature 1939, 144, 327–328. Verwey, E. J. W.; Haayman, P. W. Electronic Conductivity and Transition Point of Magnetite (“Fe3O4”). Phys. Rev. B 1941, 8, 979–987. Walz, F. The Verwey Transition-a Topical Review. J. Phys. Condens. Matter 2002, 14, R285–R340. Senn, M. S.; Wright, J. P.; Attfield, J. P. Charge Order and Three-Site Distortions in the Verwey Structure of Magnetite. Nature 2012, 481, 173–176. Chen, C. H.; Cheong, S.-W.; Hwang, H. Y. Charge-Ordered Stripes in La1−xCaxMnO3 with x>0.5. J. Appl. Phys. 1997, 81, 4326–4330. Mori, S.; Chen, C. H.; Cheong, S.-W. Pairing of ChargeOrdered Stripes in (La,Ca)MnO3. Nature 1998, 392, 473–

ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(8)

(9)

(10)

(11)

(12)

(13)

(14) (15)

(16)

(17)

(18)

(19)

(20)

(21)

(22)

(23)

(24)

(25)

(26)

12

476. Tranquada, J. M.; Sternlieb, B. J.; Axe, J. D.; Nakamura, Y.; Uchida, S. Evidence for Stripe Correlations of Spins and Holes in Copper Oxide Superconductors. Nature 1995, 375, 561–563. Andre, D.; Kim, S.-J.; Lamp, P.; Lux, S. F.; Maglia, F.; Paschos, O.; Stiaszny, B. Future Generations of Cathode Materials: An Automotive Industry Perspective. J. Mater. Chem. A 2015, 3, 6709–6732. Radin, M. D.; Hy, S.; Sina, M.; Fang, C.; Liu, H.; Vinckeviciute, J.; Zhang, M.; Whittingham, M. S.; Meng, Y. S.; Van der Ven, A. Narrowing the Gap Between Theoretical and Practical Capacities in Li-Ion Layered Oxide Cathode Materials. Adv. Energy Mater. 2017, 1602888. Demourgues, A.; Gautier, L.; Chadwick, A. V.; Delmas, C. EXAFS Study of the Jahn-Teller Distortion in Layered Nickel Oxyhydroxide. Nucl. Instruments Methods Phys. Res. Sect. B 1997, 133, 39–44. Chung, J. H.; Proffen, T.; Shamoto, S.; Ghorayeb, A. M.; Croguennec, L.; Tian, W.; Sales, B. C.; Jin, R.; Mandrus, D.; Egami, T. Local Structure of LiNiO2 Studied by Neutron Diffraction. Phys. Rev. B 2005, 71, 64410. Garc, E. Study of the Ferrodistorsive Orbital Ordering in NaNiO2 by Neutron Diffraction and Submillimeter Wave ESR. Eur. Phys. J. B 2000, 17, 615–622. Dyer, D.; Pedro, G. Alkali Metal-Nickel Oxides of the Type MNiO2. J. Am. Chem. Soc. 1954, 76, 1499. Capitaine, F.; Gravereau, P.; Delmas, C. A New Variety of LiMnO2 with a Layered Structure. Solid State Ionics 1996, 89, 197–202. Armstrong, A. R.; Bruce, P. G. Synthesis of Layered LiMnO2 as an Electrode for Rechargeable Lithium Batteries. Nature 1996, 381, 499–500. Giot, M.; Chapon, L. C.; Androulakis, J.; Green, M. A.; Radaelli, P. G.; Lappas, A. Magnetoelastic Coupling and Symmetry Breaking in the Frustrated Antiferromagnet α−NaMnO2 . Phys. Rev. Lett. 2007, 99, 247211. Zorko, A.; El Shawish, S.; Arčon, D.; Jagličić, Z.; Lappas, A.; van Tol, H.; Brunel, L. C. Magnetic Interactions in α−NaMnO2: Quantum Spin-2 System on a Spatially Anisotropic Two-Dimensional Triangular Lattice. Phys. Rev. B 2008, 77, 24412. Marianetti, C.; Morgan, D.; Ceder, G. First-Principles Investigation of the Cooperative Jahn-Teller Effect for Octahedrally Coordinated Transition-Metal Ions. Phys. Rev. B 2001, 63, 224304. Rougier, A.; Delmas, C.; Chadwick, A. V. Non-Cooperative Jahn-Teller Effect in LiNiO2: An EXAFS Study. Solid State Commun. 1995, 94, 123–127. Nakai, I.; Takahashi, K.; Shiraishi, Y.; Nakagome, T.; Nishikawa, F. Study of the Jahn-Teller Distortion in LiNiO2, a Cathode Material in a Rechargeable Lithium Battery, by in Situ X-Ray Absorption Fine Structure Analysis. J. Solid State Chem. 1998, 140, 145–148. Cao, J.; Zou, H.; Guo, C.; Chen, Z.; Pu, S. Local Trimer Orbital Ordering in LiNiO2 Studied by Quantitative Convergent Beam Electron Diffraction Technique. Solid State Ionics 2009, 180, 1209–1214. Reynaud, F.; Mertz, D.; Celestini, F.; Debierre, J.-M.; Ghorayeb, A.; Simon, P.; Stepanov, A.; Voiron, J.; Delmas, C. Orbital Frustration at the Origin of the Magnetic Behavior in LiNiO2. Phys. Rev. Lett. 2001, 86, 3638–3641. Barton, P. T.; Premchand, Y. D.; Chater, P. A.; Seshadri, R.; Rosseinsky, M. J. Chemical Inhomogeneity, Short-Range Order, and Magnetism in the LiNiO2-NiO Solid Solution. Chem. - A Eur. J. 2013, 19, 14521–14531. Chen, H.; Freeman, C. L.; Harding, J. H. Charge Disproportionation and Jahn-Teller Distortion in LiNiO2 and NaNiO2: A Density Functional Theory Study. Phys. Rev. B 2011, 84, 85108. Chen, Z.; Zou, H.; Zhu, X.; Zou, J.; Cao, J. First-Principle Investigation of Jahn-Teller Distortion and Topological

(27)

(28) (29) (30)

(31)

(32)

(33)

(34)

(35)

(36)

(37)

(38)

(39)

(40)

(41) (42)

(43)

(44)

(45)

(46)

(47) (48) (49)

Page 12 of 14

Analysis of Chemical Bonds in LiNiO2. J. Solid State Chem. 2011, 184, 1784–1790. Vernay, F.; Penc, K.; Fazekas, P.; Mila, F. Orbital Degeneracy as a Source of Frustration in LiNiO2. Phys. Rev. B - Condens. Matter Mater. Phys. 2004, 70, 1–17. Kresse, G.; Marsman, M.; Furthmüller, J. VASP the Guide. 2016. Ong, S. P.; Mo, Y.; Ceder, G. Low Hole Polaron Migration Barrier in Lithium Peroxide. Phys. Rev. B 2012, 85, 81105. Radin, M. D.; Siegel, D. J. Charge Transport in Lithium Peroxide: Relevance for Rechargeable Metal-Air Batteries. Energy Environ. Sci. 2013, 6, 2370–2379. Maxisch, T.; Zhou, F.; Ceder, G. Ab Initio Study of the Migration of Small Polarons in Olivine LixFePO4 and Their Association with Lithium Ions and Vacancies. Phys. Rev. B 2006, 73, 104301. Ong, S. P.; Chevrier, V. L.; Ceder, G. Comparison of Small Polaron Migration and Phase Separation in Olivine LiMnPO4 and LiFePO4 Using Hybrid Density Functional Theory. Phys. Rev. B 2011, 83, 075112. Johannes, M. D.; Hoang, K.; Allen, J. L.; Gaskell, K. Hole Polaron Formation and Migration in Olivine Phosphate Materials. Phys. Rev. B 2012, 85, 115106. Meredig, B.; Thompson, A.; Hansen, H. a.; Wolverton, C.; van de Walle, A. Method for Locating Low-Energy Solutions within DFT+U. Phys. Rev. B 2010, 82, 195128. Allen, J. P.; Watson, G. W. Occupation Matrix Control of Dand F-Electron Localisations Using DFT + U. Phys. Chem. Chem. Phys. 2014, 16, 21016–21031. Dorado, B.; Amadon, B.; Freyss, M.; Bertolus, M. DFT+U Calculations of the Ground State and Metastable States of Uranium Dioxide. Phys. Rev. B 2009, 79, 235125. Arroyo y de Dompablo, M. E.; Marianetti, C.; Van der Ven, A.; Ceder, G. Jahn-Teller Mediated Ordering in Layered LixMO2 Compounds. Phys. Rev. B 2001, 63, 144107. Li, X.; Ma, X.; Su, D.; Liu, L.; Chisnell, R.; Ong, S. P.; Chen, H.; Toumar, A.; Idrobo, J.; Lei, Y.; et al. Direct Visualization of the Jahn-Teller Effect Coupled to Na Ordering in Na5/8MnO2. Nat. Mater. 2014, 13, 586–592. Stock, C.; Chapon, L. C.; Adamopoulos, O.; Lappas, A.; Giot, M.; Taylor, J. W.; Green, M. A.; Brown, C. M.; Radaelli, P. G. One-Dimensional Magnetic Fluctuations in the Spin-2 Triangular Lattice α−NaMnO2. Phys. Rev. Lett. 2009, 103, 77202. Yabuuchi, N.; Kubota, K.; Dahbi, M.; Komaba, S. Research Development on Sodium-Ion Batteries. Chem. Rev. 2014, 114, 11636–11682. Ma, X.; Chen, H.; Ceder, G. Electrochemical Properties of Monoclinic NaMnO2. J. Electrochem. Soc. 2011, 158, A1307. Han, M. H.; Gonzalo, E.; Casas-Cabanas, M.; Rojo, T. Structural Evolution and Electrochemistry of Monoclinic NaNiO2 upon the First Cycling Process. J. Power Sources 2014, 258, 266–271. Clarke, S. J.; Fowkes, A. J.; Harrison, A.; Ibberson, R. M.; Rosseinsky, M. J. Synthesis, Structure, and Magnetic Properties of NaTiO2. Chem. Mater. 1998, 10, 372–384. Dhariwal, M.; Pisani, L.; Maitra, T. Competing Electronic States in High Temperature Phase of NaTi2. J. Phys. Condens. Matter 2014, 26, 205501. McQueen, T. M.; Stephens, P. W.; Huang, Q.; Klimczuk, T.; Ronning, F.; Cava, R. J. Successive Orbital Ordering Transitions in NaVO2. Phys. Rev. Lett. 2008, 101, 2–5. Aydinol, M. K.; Kohan, A. F.; Ceder, G.; Cho, K.; Joannopoulos, J. Ab Initio Study of Lithium Intercalation in Metal Oxides and Metal Dichalcogenides. Phys. Rev. B 1997, 56, 1354–1365. Bersuker, I. The Jahn-Teller Effect; Cambridge University Press: Cambridge, 2006. Kanamori, J. Crystal Distortion in Magnetic Compounds. J. Appl. Phys. 1960, 31, S14. Khomskii, D.; van den Brink, J. Anharmonic Effects on Charge and Orbital Order. Phys. Rev. Lett. 2000, 85, 3329.

ACS Paragon Plus Environment

Page 13 of 14 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(50)

(51)

(52)

(53)

(54)

(55) (56)

(57)

(58)

(59)

(60) (61)

(62)

(63)

(64) (65)

(66)

(67)

(68)

(69)

(70)

(71)

13

Chemistry of Materials Sanchez, J. M.; Ducastelle, F.; Gratias, D. Generalized Cluster Description of Multicomponent Systems. Phys. A Stat. Mech. its Appl. 1984, 128, 334–350. Ducastelle, F. Order and Phase Stability in Alloys; NorthHolland ; Sole distributors for the USA and Canada, Elsevier Science Pub. Co.: Amsterdam; New York; New York, NY, USA, 1991. Fontaine, D. De. Cluster Approach to Order-Disorder Transformations in Alloys. Solid State Phys. 1994, 47, 33– 176. Zhou, F.; Maxisch, T.; Ceder, G. Configurational Electronic Entropy and the Phase Diagram of Mixed-Valence Oxides: The Case of LixFePO4. Phys. Rev. Lett. 2006, 97, 155704. Wang, Y. C.; Chen, Z. H.; Jiang, H. The Local Projection in the Density Functional Theory plus U Approach: A Critical Assessment. J. Chem. Phys. 2016, 144, 144106. Thomas, J. C; Van der Ven, A. In preparation. CASM, v0.1.0 (2015). Available from https://github.com/prisms-center/CASMcode. doi: 10.5281/zenodo.31216 Thomas, J. C.; Van Der Ven, A. Finite-Temperature Properties of Strongly Anharmonic and Mechanically Unstable Crystal Phases from First Principles. Phys. Rev. B 2013, 88, 214111. Puchala, B.; Van Der Ven, A. Thermodynamics of the Zr-O System from First-Principles Calculations. Phys. Rev. B 2013, 88, 94108. Van der Ven, A.; Thomas, J. C.; Xu, Q.; Bhattacharya, J. Linking the Electronic Structure of Solids to Their Thermodynamic and Kinetic Properties. Math. Comput. Simul. 2010, 80, 1393–1410. Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 1993, 47, 558–561. Kresse, G.; Hafner, J. Ab Initio Molecular-Dynamics Simulation of the Liquid-Metal-Amorphous-Semiconductor Transition in Germanium. Phys. Rev. B 1994, 49, 14251– 14269. Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15–50. Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Eave Basis Set. Phys. Rev. B 1996, 54, 11169–11186. Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953. Blöchl, P.; Jepsen, O.; Andersen, O. Improved Tetrahedron Method for Brillouin-Zone Integrations. Phys. Rev. B 1994, 49, 16223–16233. Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. Dudarev, S.; Botton, G.; Savrasov, S.; Humphreys, C.; Sutton, A. Electron-Energy-Loss Spectra and the Structural Stability of Nickel Oxide: An LSDA+ U Study. Phys. Rev. B 1998, 57, 1505–1509. Wang, L.; Maxisch, T.; Ceder, G. Oxidation Energies of Transition Metal Oxides within the GGA+U Framework. Phys. Rev. B 2006, 73, 195107. Zhou, F.; Cococcioni, M.; Marianetti, C.; Morgan, D.; Ceder, G. First-Principles Prediction of Redox Potentials in Transition-Metal Compounds with LDA+U. Phys. Rev. B 2004, 70, 235121. Heyd, J.; Scuseria, G. E.; Ernzerhof, M. Hybrid Functionals Based on a Screened Coulomb Potential. J. Chem. Phys. 2003, 118, 8207. Krukau, A. V; Vydrov, O. a; Izmaylov, A. F.; Scuseria, G. E. Influence of the Exchange Screening Parameter on the Performance of Screened Hybrid Functionals. J. Chem. Phys. 2006, 125, 224106.

(72)

(73)

(74)

(75) (76)

(77) (78)

(79)

(80)

(81)

(82)

(83)

(84)

(85)

(86)

(87)

(88)

(89)

(90)

Chernova, N. A.; Nolis, G. M.; Omenya, F. O.; Zhou, H.; Li, Z.; Whittingham, M. S. What Can We Learn about Battery Materials from Their Magnetic Properties? J. Mater. Chem. 2011, 21, 9865. Chappel, E.; Núñez-Regueiro, M. D.; de Brion, S.; Chouteau, G.; Bianchi, V.; Caurant, D.; Baffier, N. Interlayer Magnetic Frustration in Quasistoichiometric Li1-xNi1+xO2. Phys. Rev. B 2002, 66, 132412. Meskine, H.; Satpathy, S. Electronic Structure and Magnetism in Sodium Nickelate: Density-Functional and Model Studies. Phys. Rev. B 2005, 72, 224423. Borgers, P. F.; Enz, U. Metamagnetism of NaNiO2. Solid State Commun. 1966, 4, 153–157. Chappel, E.; Dupont, F.; Chouteau, G.; Darie, C.; Sulpice, A. Antiferromagnetic Resonance and High Magnetic Field Properties of NaNiO2. Eur. Phys. J. B 2000, 17, 609–614. Singh, D. J. Magnetic and Electronic Properties of LiMnO2. Phys. Rev. B 1997, 55, 309–312. Mackrodt, W. C.; Williamson, E. A. First-Principles Hartree-Fock Description of the Electronic Structure of Monoclinic C2/m LixMnO2 (1 >= x >= 0). Philos. Mag. B 1998, 77, 1077–1092. Jang, Y.-I.; Chou, F. C.; Chiang, Y.-M. Magnetic Properties of Monoclinic Phase LiAl0.05Mn0.95O2. J. Phys. Chem. Solids 1999, 60, 1763–1771. Velikokhatnyi, O. I.; Chang, C.-C.; Kumta, P. N. Phase Stability and Electronic Structure of NaMnO2. J. Electrochem. Soc. 2003, 150, A1262. Abakumov, A. M.; Tsirlin, A. A.; Bakaimi, I.; Van Tendeloo, G.; Lappas, A. Multiple Twinning as a Structure Directing Mechanism in Layered Rock-Salt-Type Oxides: NaMnO2 Polymorphism, Redox Potentials, and Magnetism. Chem. Mater. 2014, 26, 3306–3315. Delmas, C.; Pérès, J. P.; Rougier, A.; Demourgues, A.; Weill, F.; Chadwick, A.; Broussely, M.; Perton, F.; Biensan, P.; Willmann, P. On the Behavior of the LixNiO2 System : An Electrochemical and Structural Overview. J. Power Sources 1997, 68, 120–125. Delmas, C.; Ménétrier, M.; Croguennec, L.; Saadoune, I.; Rougier, A.; Pouillerie, C.; Prado, G.; Grune, M.; Fournès, L. An Overview of the Li(Ni,M)O2 Systems: Syntheses , Structures and Properties. Electrochim. Acta 1999, 45, 243– 253. Radin, M. D.; Van der Ven, A. The Stability of Prismatic and Octahedral Coordination in Layered Oxides and Sulfides Intercalated with Alkali and Alkaline-Earth Metals. Chem. Mater. 2016, 28, 7898–7904. Janson, V. M.; Hoppe, R. Zur Kenntnis Der NaClStrukturfamilie Die Kristallstruktur von NaMnO2. Z. Anorg. Allg. Chem. 1973, 399, 163–169. Zorko, A.; Adamopoulos, O.; Komelj, M.; Arčon, D.; Lappas, A. Frustration-Induced Nanometre-Scale Inhomogeneity in a Triangular Antiferromagnet. Nat. Commun. 2014, 5, 3222. Yoon, W.-S.; Balasubramanian, M.; Chung, K. Y.; Yang, X.Q.; McBreen, J.; Grey, C. P.; Fischer, D. A. Investigation of the Charge Compensation Mechanism on the Electrochemically Li-Ion Deintercalated Li1−xCo1/3Ni1/3Mn1/3O2 Electrode System by Combination of Soft and Hard X-Ray Absorption Spectroscopy. J. Am. Chem. Soc. 2005, 127, 17479–17487. Hirakawa, K.; Kadowaki, H.; Ubukoshi, K. Experimental Studies of Triangular Lattice Antiferromagnets with S=1/2: NaTiO2 and LiNiO2. J. Phys. Soc. Japan 1985, 54, 3526–3536. Trease, N. M.; Seymour, I. D.; Radin, M. D.; Liu, H.; Liu, H.; Hy, S.; Chernova, N.; Parikh, P.; Devaraj, A.; Wiaderek, K. M.; et al. Identifying the Distribution of Al3+ in LiNi0.8Co0.15Al0.05O2. Chem. Mater. 2016, 28, 8170–8180. Khachaturyan, A. G. Theory of Structural Transformations in Solids.; Dover Publications: Newburyport, 2013.

ACS Paragon Plus Environment

Chemistry of Materials 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 14

TABLE OF CONTENTS GRAPHIC

14 ACS Paragon Plus Environment