Magneto-orbital helices as a route to coupling magnetism and ferroelectricity in multiferroic CaMn7O12

Abstract

Orbital physics drives a rich phenomenology in transition-metal oxides, providing the microscopic underpinning for effects such as Colossal Magnetoresistance. In particular, magnetic and lattice degrees of freedom are coupled through orbital ordering, and it has long been hoped that this coupling could be exploited to create high-temperature multiferroics with large values of the electrical polarization. Here we report an unprecedented magneto-orbital texture in multiferroic CaMn₇ O₁₂ , found to give rise to the largest magnetically induced ferroelectric polarization measured to date. X-ray diffraction characterization of the structural modulation in these 'magneto-orbital helices', and analysis of magnetic exchange shows that orbital order is crucial in stabilising a chiral magnetic structure, thus allowing for electric polarization. Additionally, the presence of a global structural rotation enables the coupling between this polarization and magnetic helicity required for multiferroicity. These novel principles open up the possibility of discovering new multiferroics with even larger polarization and higher transition temperatures.

ARTICLE

Received 27 Jul 2012 | Accepted 14 Nov 2012 | Published 11 Dec 2012

Orbital physics drives a rich phenomenology in transition-metal oxides, providing the microscopic underpinning for effects such as Colossal Magnetoresistance. In particular, magnetic and lattice degrees of freedom are coupled through orbital ordering, and it has long been hoped that this coupling could be exploited to create high-temperature multiferroics with large values of the electrical polarization. Here we report an unprecedented magneto-orbital texture in multiferroic CaMn7O12, found to give rise to the largest magnetically induced ferroelectric polarization measured to date. X-ray diffraction characterization of the structural modulation in these magneto-orbital helices, and analysis of magnetic exchange shows that orbital order is crucial in stabilising a chiral magnetic structure, thus allowing for electric polarization. Additionally, the presence of a global structural rotation enables the coupling between this polarization and magnetic helicity required for multiferroicity. These novel principles open up the possibility of discovering new multiferroics with even larger polarization and higher transition temperatures.

DOI: 10.1038/ncomms2294

Magneto-orbital helices as a route to coupling magnetism and ferroelectricity in multiferroic CaMn7O12

N.J. Perks1, R.D. Johnson1,2, C. Martin3, L.C. Chapon4 & P.G. Radaelli1

1 Department of Physics, Clarendon Laboratory, University of Oxford, Oxford OX1 3PU, UK. 2 ISIS facility, Rutherford Appleton Laboratory-STFC, Chilton, Didcot OX11 0QX, UK. 3 Laboratoire CRISMAT, ENSICAEN, UMR F-6508 CNRS, 6 Boulevard du Marechal Juin, F-14050 Caen, France.

4 Institut Laue-Langevin, BP 156X, 38042 Grenoble, France. Correspondence and requests for materials should be addressed to N.J.P. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 3:1277 | DOI: 10.1038/ncomms2294 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms2294

Orbital ordering (OO) is one of the most striking manifestations of electronic correlations in a variety of materials, most notably in 3d transition-metal oxides1.

In these systems, below a certain temperature TOO, a distinct

anisotropic pattern of orbital occupation emerges from a nearly isotropic high-temperature state. As OO is strongly connected with both structural distortions (the cooperative Jahn Teller effect) and with magnetism, OO has long been regarded as a promising ingredient to allow the cross-coupling of electric properties. Recently, it has been proposed that some charge and orbitally ordered states could be electrically polar, so that OO could promote high-temperature, magnetically controllable ferroelectricity2. Here we show that a unique form of incommensurate OO in the complex manganese oxide CaMn7O12, appearing below 250 K, is responsible for the largest observed magnetically induced ferroelectric polarization. Unlike all other known systems, orbital occupation in CaMn7O12 does

not select xed preferential directions, but rotates around the threefold crystallographic axis, preserving the overall point-group symmetry. The resulting pattern of magnetic exchange interactions produces a magnetic helix having twice the period of the orbital modulation below 90 K.

Upon cooling below B400 K, CaMn7O12 adopts a trigonal structure (space group R 3), shown in Fig. 1a. It can be described as a perovskite superstructure, with 3/4 of the A sites occupied by Mn3 (hereafter labelled Mn1), and the remaining occupied by

Ca. The pseudocubic B sites form three inequivalent triangular

layers stacked along the hexagonal c axis (shown in Fig. 1a, and coloured green, blue and red throughout this report). Within each layer, 1/4 of the B sites retain 3 symmetry (Mn3, black-outlined circles in Fig. 1a), whereas the other 3/4 have 1 symmetry (Mn2). Charge ordering develops gradually on further cooling, with Mn2 and Mn3 tending towards Mn3 and Mn4 , respectively. An incommensurate structural modulation along the hexagonal c axis appears below TOO 250 K, with a propagation vector qc (0, 0,

2.0770.001) at 150 K, producing a continuous variation in Mn-O bond lengths3.

The magnetic structure below TN1 90 K was determined

using neutron powder diffraction4, and is shown in Fig. 1b. Magnetic anisotropy results in the spins lying within the hexagonal ab plane. Within each of the triangular layers, Mn2 spins are parallel to each other, while Mn3 spins are approximately antiparallel. Adjacent layers are rotated by 1241 in the ab plane with respect to each other, leading to an incommensurate magnetic modulation with propagation vector qm (0, 0, 1.037). As observed in many magnetically ordered

crystals (see, for example ref. 5), below TN1, qc 2qm. Here the

structural modulation is also present in the paramagnetic phase, and the magnetic modulation locks in to it as it develops at TN1.

Below TN2 50 K, another magnetic transition occurs, with the

appearance of a second propagation vector4. S"awiski et al.6 have shown that the average propagation vector is constant across TN2,

and described the low-temperature magnetic phase as a beating of the high-temperature phasea model that provides a good t

c

b

a

c

b

a

Mn2 layer A Mn2 layer B Mn2 layer C

Mn3 layer A Mn3 layer B Mn3 layer C

Figure 1 | The crystal and magnetic structures and associated orbital orientations of CaMn7O12. (a) (left panel) The room temperature crystal structure in the hexagonal setting, highlighting the three triangular layers (A, blue, B, red and C, green) stacked along the hexagonal c axis (cubic [111] direction). The structure can also be described within the pseudocubic perovskite setting (right panel), with Mn3 sites (black-outlines spheres, with nominal valence 4 )

occupying the apex sites along the [111] diagonal. (b) The magnetic structure at 90 K. The triangular layers are again shown as blue, red and green, with a 1241 spin rotation between adjacent layers. (c) A schematic showing the B sites in the pseudocubic setting, dening the x, y and z bond directions, which are used to describe orbital orientation.

Mn2 x bond Mn2 y bond Mn2 z bond

2 NATURE COMMUNICATIONS | 3:1277 | DOI: 10.1038/ncomms2294 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2012 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms2294 ARTICLE

to their neutron data. Although the origin of such complex low-temperature structure remains obscure, its symmetry properties (in particular its chirality) are very similar to the high-temperature phase, and the transition at TN2 has clearly very little effect on the electrical properties4. In this paper, we discuss primarily the simpler magnetic phase between TN1 and TN2. In this phase, the magnetic structure induced by the orbital texture breaks inversion symmetry and couples to the crystal structure via the relativistic antisymmetric exchange to yield the large switchable electrical polarization.

ResultsThe structural modulation. We now focus on the incommensurate structural modulation appearing below TOO, and demon

strate that this corresponds to the onset of a novel form of OO. Single-crystal X-ray diffraction data were measured at 150 K. Many different crystals from different batches were tested, consistently yielding sharp superstructure peaks, indicative of a long-range-ordered modulated structure. The modulated structure was rened using JANA20067. Consistent with previous reports3, we found that the modulated crystal structure is described by the centrosymmetric four-dimensional space group R 3(00g)0 with propagation vector qc (0, 0, 2.077). In Fig. 2a, we plot the Mn-O

bond lengths for adjacent unit cells along the c axis around Mn2 (for variation in the Mn1-O and Mn3-O bond lengths see Supplementary Fig. S1). This structural modulation has strong implications for the nature of orbital occupation as follows. We dene orthogonal x, y and z bond directions, schematically shown in Fig. 1c. In the unmodulated, average structure, the Mn3 O6 octahedra around Mn2 are compressed, giving two short and four long bonds, implying an unusual occupation of the x2 y2

orbital. In reality, the structural modulation leads to an incommensurate rotation of a single pair of long bonds in the xy plane. We can characterize this orbital rotation by employing a formalism rst proposed by Goodenough8,9, which quanties orbital occupation according to a mixing angle, tan y

3

p x z/2y x z. Following the conventions in

Wu et al.9, we need only to consider the angular sector 2p/3ryr4p/3 (inset to Fig. 2). As shown in Fig. 2b, the orbital occupation rotates periodically between the 3x2 r2 and

3y2 r2 orbitals, passing through the x2 y2 only at the nodal

points. Fig. 2c shows the Mn2 charge modulation accompanying the orbital modulation, calculated using the Bond Valence Sum method10. Valence uctuations are surprisingly small (1%see Supplementary Fig. S2 and Supplementary Discussion), and may even be entirely an artefact of the model that does not account for the modulation of thermal parameters. This indicates that the character of the modulation is predominantly orbital.

The magneto-orbital helix. Although OO is a widespread phenomenon11,12, an incommensurate modulation where the orbital occupation rotates perpendicular to the propagation direction is unprecedented. The resulting orbital arrangement has a profound effect on magnetic ordering below TN 90 K. This can be

seen by applying the GoodenoughKanamoriAnderson (GKA) rules1316, which predict the nature of exchange interactions. In the absence of the structural modulation, the pseudo-perovskite structure is heavily distorted, and antiferromagnetic interactions are weak due to the large departure of the Mn-O-Mn bond angles from 1801. The alternation of lled and empty x2 y2 orbitals on

Mn2 and Mn3 promotes ferromagnetic exchange, and a net ferromagnetic structure is predicted. This is conrmed by rst-principle calculations17, which also show that a rotating magnetic structure would require spinorbit coupling. The situation changes drastically in the presence of the orbital modulation,

which enables the helical magnetic structure to be stabilised by Heisenberg exchange and to lock to the superstructure. This is because the magnetic interaction between adjacent A, B and C layers is predominantly mediated by the orbitals that are most affected by the modulation, which introduces magnetic frustration. This is shown schematically in Fig. 3a, which depicts a region of the crystal structure where the modulation amplitude is largest. Mn3 in the B layer is connected to the Mn2 sites in the A layer by lled 3x2 r2 orbitals, yielding a

ferromagnetic interaction J1, and to the Mn2 sites in the C layer by empty orbitals, yielding an antiferromagnetic interaction J2 (crucially, J1 and J2 would be equal in the absence of OO). This produces an effective ferromagnetic interaction between spins in each of the A, B and C layers, as all J1(J2)s are equivalent by symmetry (shown in Fig. 3b). Across layers there is an effective next-nearest-neighbour antiferromagnetic interaction J12, which

is strongest at the crest of the modulation, but remains antiferromagnetic everywhere except possibly at the nodal point. The overall nearest-neighbour interaction between layers, J3, results from the competition between full-empty (x z)

orbitals and empty-empty (y z) orbitals, which are ferro- and

antiferromagnetic, respectively. Although the sign of J3 cannot be predicted a priori, it is experimentally found to be positive (antiferromagnetic). The magnetic helix, therefore, results from the competition between antiferromagnetic nearest-neighbour and next-nearest-neighbour interactions. We can test this hypothesis with a very simple model, which captures the effect of the average interactions J1, J2 and J3, assuming a magnetic propagation vector (0,0,qz) and identical spins S on Mn2 and Mn3 sites. The magnetic Heisenberg exchange energy per Mn3 site can be written as:

E S2 3J1 cos f 3J2 cosf 2c 6J3 cos c

1 where c 2/3pqz and f are the angles between Mn2-Mn2 and

Mn3-Mn2 spins in adjacent layers, respectively. Minimizing the energy in equation 1 yields the phase diagram shown in Fig. 3c. The wide incommensurate phase eld is bound on either side by two distinct commensurate collinear phases, both with qz 1.5.

The white line indicates the range of parameters for which the most stable propagation vector is that determined experimentally. The red dot corresponds precisely to the experimental magnetic structure (f 311), and is found for J2/J1 0.85 and

J3/J1 0.63, in good agreement with our qualitative

considerations based on the GKA rules. The coupling between Mn1 and adjacent Mn2 sites is expected (and observed) to be ferromagnetic (Fig. 3a) as the sites are linked via oxygen anions with highly-distorted bonds at less than the GKA critical angle. The most important ingredient not included in our simple model is the qc 2qm lock-in term, which intertwines the magnetic

and orbital modulations to form a magneto-orbital helix. Phenomenologically, different forms of the lock-in term are allowed by symmetry, with the most likely candidate involving a third-harmonic squaring up of the magnetic modulation.

DiscussionBelow TN1, CaMn7O12 acquires the largest measured magnetically induced electrical polarization along the c axis4, which can be explained by the existence of the magneto-orbital helix. CaMn7O12 has point symmetry 3one of the ferroaxial groups.

Ferroaxiality is a ferroic property like ferroelectricity and ferromagnetism, characterized by the presence of domains having a well-dened rotation between different elements of the crystal structure. In CaMn7O12, the six oxygen atoms surrounding an individual Mn3 are all rotated in the same direction (clockwise in the domain illustrated in Fig. 4a). The electrical polarization can arise through coupling between the

NATURE COMMUNICATIONS | 3:1277 | DOI: 10.1038/ncomms2294 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms2294

2.15

2.1

2.05

Mn2-0 bond length ()

3

c z

x

2

1.95

225

200

1

2 3

y

2

[afii9835] (Degrees)

Mn2 bond valence sum (+e)

175

150

125

100

3.16

3.14

3.12

3.1

3.08

3.06

3y 2r 2 >

y 2z 2 >

120

x 2y 2 >

3x 2r 2 > z 2x 2 >

[afii9835]

0 3z 2r 2 >

c

0

2

4 6 8

240

10 12 14 16 18

1

R z

Figure 2 | The incommensurate structural modulation coupled to orbital rotation and valence uctuations. (a) Variation in Mn2-O bond lengths along x, y and z (coloured according to the schematic in Fig. 1c) as a function of Rz, plotted across 19 unit cells along the c axis. (b) Orbital mixing angle y for the Mn2 octahedra as a function of Rz. y correlates directly with orbital occupation, as shown by the circular inset. The orbital rotation is shown, along with the coupled helical spin rotation. (c) The Mn2 charge modulation, which accompanies the structural modulation along the c axis. The incommensurate orbital modulation for half a period of the structural modulation is depicted at the side of the gure, with positions 1, 2 and 3 corresponding to the labelled points in (a).

chiral magneto-orbital helix and the global rotation, as described by the phenomenological invariant asA P, where a is a coupling

constant, s is the magnetic helicity, P is the electrical polarization and A is an axial vector representing the structural rotation. It is important to stress that the helicity of this particular magnetic structure can be reversed by a global rotation of all the spins, so that the symmetric-exchange energy of two inversion-related domains must be the same. This implies that the spinorbit interaction is essential to couple magnetic helicity and polarization, although striction effects can be important in achieving large values of P17. The relevant microscopic mechanism is the antisymmetric DzyaloshinskiiMoriya (DM) exchange, where energy can be gained by distorting the crystal structure and/or the electronic density in the presence of noncollinear spin congurations. Although both structural and electronic distortions are allowed in the present case, we will employ the former as an illustration. When two metal sites

carrying non-collinear spins are joined by common ligand atoms (in this case Mn-O-Mn), energy can be gained by displacing the ligand through a vector u so that DE

lu [r12 (S1 S2)] lD (S1 S2), where S1 and S2 are the

spins on the two sites, r12 is the position vector connecting them and l is a coupling constant. D u r12 is the familiar DM

vector18,19. Energy is, therefore, minimized by a pattern of local ligand (oxygen) displacements u associated with pairs of spins and either parallel or antiparallel to the vectors r12 (S1 S2),

depending on the sign of l. Figure 4b shows the pattern of r12 (S1 S2) vectors (black arrows) for clusters of Mn2 around

a single Mn3. Assuming lo0, the DM interaction will favour oxygen displacements parallel to these arrows, making the Mn3 -O-Mn4 bonds between the red and the blue layer atter (that is, the bond angles closer to 1801) and Mn3-O

Mn4 bonds between the red and the green layer more acute (that is, the bond angles farther from 1801). It is easy to see

4 NATURE COMMUNICATIONS | 3:1277 | DOI: 10.1038/ncomms2294 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2012 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms2294 ARTICLE

J1

J12

J2

Mn2 layer A Mn2 layer B

Mn2 layer C

Mn3 layer A Mn3 layer B Mn1 (A-site)

J3

c

ab

3.0

CM2

J1

CM1

J2

J2

J1

J 3/J 1

1.50

1.25

1.00

0.75

k z

J1

2.0

1.0

0.0

ICM

J2

b

J3

a

0.0 1.0 2.0 3.0 J2/J1

Figure 3 | The effect of OO on magnetic exchange interactions promotes a chiral magnetic structure. (a) Three pseudocubic unit cells illustrating the coupling between the incommensurate orbital modulation and the magnetic structure. The occupied Mn2 (Mn3 ) orbitals are shown in grey, giving rise to ferromagnetic exchange interactions (J1), shown by a dashed line, and antiferromagnetic exchange interactions (J2), shown by a dotted line. The nearest-neighbour and effective next-nearest-neighbour antiferromagnetic exchanges between layers are labelled J3 and J12, respectively. The ferromagnetic exchange along the Mn1Mn2 chains is shown by a solid black line. (b) Figure 3a viewed down the c axis, showing the symmetry equivalent exchange interactions. (c) Phase diagram of the CaMn7O12 magnetic structures, obtained by minimizing the classical Heisenberg exchange energy in equation 1 as a function of the angles c and f for different exchange constants J1, J2 and J3. Spins on Mn2 and Mn3 sites were taken to be equal (different spins would be equivalent to renormalising the exchange constants). CM1 and CM2 are two commensurate collinear phases, both with propagation vector (0, 0, 1.5), in which the Mn2 and Mn3 in each layer are antiparallel (for CM1) and parallel (for CM2). ICM is an incommensurate phase with rotating spins (magnetic helix). The shading indicates the value of the z-component of the propagation vector. The white line indicates the locus of points with the experimental propagation vector (0, 0, 1.037). The red dot corresponds to the experimental magnetic structure (f 311).

that this pattern of displacements will result in the central Mn4 moving upwards (yellow arrow), creating a local polarization along the positive c direction (identical considerations for l40 lead to a polarization along the negative c direction). As the magnetic structure is globally chiral, these local polarizations add up to give a net electrical polarization. It is important to stress that such a simple relation between local u and global P vectors is only possible in the presence of a global structural rotation, which is allowed in point-group 3 (no net polarization would arise for 1801 Mn3 -O-Mn4 bonds). Ferroaxiality is thus an essential ingredient for the development of electrical polarization in

CaMn7O12. Moreover, the polarization is constrained to lie along the high-symmetry axial vector, enforcing a simple switchable 2-domain structure, where the magnetic chirality determines the direction of P (the ferroaxiality itself does not switch under applied electric eld). The large value of P is readily explained by considering the angle of rotation f between the spins S1 and S2, as P should be roughly proportional to sinf. The average angle of

rotation is /fS 124/21 621, whereas for the prototypical

magnetic multiferroic TbMnO3, f 2001, so that sin/fS in

CaMn7O12 is 2.5 times larger than in TbMnO3. It is also clear to see that the small deviation from commensurability is not crucial in determining the orbital physics and the large electrical polarization in CaMn7O12: P would be almost as large if the magnetic propagation vector was exactly (0,0,1), as what matters is the large rotation angle between spins in adjacent A, B and C layers. This model should also be applicable to the complex low-temperature magnetic structure below TN2, as the modulated

structure is unaffected by the transition and the magnetic structure is still chiral.

We have shown that the magnetic and structural modulations in CaMn7O12 are intertwined to form an incommensurate magneto-orbital helix. This unprecedented texture represents a novel manifestation of the deep connection between orbital physics and magnetism in transition-metal oxides. The structural modulation is crucial in stabilising a chiral magnetic structure,

NATURE COMMUNICATIONS | 3:1277 | DOI: 10.1038/ncomms2294 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2012 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms2294

A

b

a

Mn 3+

Mn 4+

Mn 3+

P

Figure 4 | The coupling between magnetic helicity and ferroelectricity leading to a giant polarization. (a) The global structural rotation described by the axial vector A along the c axis is evident in the pattern of oxygen atoms (red spheres), which couples to the clockwise (anticlockwise) spin rotation between Mn4 and green (blue) Mn3 layers. (b) A perspective view on the Mn3 Mn4 Mn3 conguration. The inverse

DzyaloshinskiiMoriya interaction favours a displacement of the oxygen anions in the direction of the black arrows. This results in an additional displacement of the Mn4 cation, indicated by the small yellow arrow, and thus a net polarization along the positive c axis, shown by the large yellow arrow.

which breaks inversion symmetry, producing the largest magnetically induced electrical polarization ever observed in any system. We have also shown how the presence of a global structural rotation (ferroaxiality) promotes the coupling between magneto-orbital helicity and ferroelectricity. These observations open up the possibility of nding a new class of multiferroic materials in which giant ferroelectric polarizations may arise, underpinning the development of multiferroic technology.

Methods

Note on magnetic propagation vector. qm (0, 0, 1.037) is equivalent to the

previously reported magnetic propagation vector (0, 1, 0.963) through addition of the reciprocal lattice vector (0,1,2) to qm. In this study, qm is chosen to lie within

the rst Brillouin zone of the hexagonal unit cell.

Experimental details. High-quality single crystals of CaMn7O12, with approximate dimensions 60 60 60 mm were grown by the ux method. A 5 g mixture of

CaCl2:MnO2 was heated for 25 h in an alumina crucible in air at 850 1C and then cooled at a rate of 5 1C per hour. X-ray diffraction and bulk property

measurements showed that the crystals were not twinned at room temperature. Single-crystal X-ray diffraction at 150 K was performed using a SuperNova X-ray diffraction system from Agilent Technologies, tted with an Oxford Cryosystems L N2 cryostream. Crysalis Pro software was used for data collection and processing, enabling concurrent data reduction and identication of additional diffraction peaks that originate in the structural modulation, observed as satellites of Bragg reections from the average structure. A powder sample was prepared by grinding single crystals, and sieving through a 35-mm mesh. This sample was used to obtain accurate lattice parameters at 150 K through X-ray powder diffraction using a Siemens D5000 diffractometer. The incommensurately modulated structure was rened using Jana (2006) software7 with lattice parameters obtained from the X-ray powder diffraction measurement, the propagation vector obtained from single-crystal X-ray diffraction, and an initial structure taken from S"awiski et al.3

References

1. Hotta, T. Orbital ordering phenomena in d- and f-electron systems. Rep. Prog. Phys. 69, 20612150 (2006).

2. Van den Brink, J. & Khomskii, D. I. Multiferroicity due to charge ordering.J. Phys. Condens. Matter 20, 434217 (2008).3. S"awiski, W. et al. Modulation of atomic positions in CaCuxMn7 xO12

(xr0.1). Acta. Cryst. 65, 535542 (2009).4. Johnson, R. D. et al. Giant improper ferroelectricity in the ferroaxial magnet CaMn7O12. Phys. Rev. Lett. 108, 067201 (2012).

5. Noda, Y. et al. Magnetic and ferroelectric properties of multiferroic RMn2O5.

J. Phys.: Condens. Matter 20, 434206 (2008).6. S"awiski, W., Przenios, R., Sosnowska, I. & Chrobak, A. Beats in the magnetic modulation of multiferroic CaMn7O12. J. Phys. Soc. Jpn 81, 094708 (2012).

7. Petricek, V., Dusek, M. & Palatinus, L. Jana. Technical Report. Institute of Physics, Prague, Czech Republic (2006).

8. Goodenough, J. B. & Cooper, S. L. Localized to intinerant: electronic transition in perovskite oxide (Springer, 2001).

9. Wu, H. et al. Orbital order in La0.5Sr1.5MnO4. Phys. Rev. B 84, 155126 (2011).

10. Brown, I. D. & Altermatt, D. Bond-valence parameters obtained from a systematic analysis of the inorganic crystal-structure database. Acta. Cryst. B41, 244247 (1985).

11. Radaelli, P. G., Cox, D. E., Marezio, M. & Cheong, S. W. Charge, orbital and magnetic ordering in La0.5Ca0.5MnO3. Phys. Rev. B 55, 30153023 (1997).

12. Feiner, F. & Oles, A. M. Electronic origin of magnetic and orbital ordering in insulating LaMnO3. Phys. Rev. B 59, 32953298 (1999).

13. Anderson, P. W. Antiferromagnetism - theory of superexchange interaction. Phys. Rev. 79, 350356 (1950).

14. Kanamori, J. Theory of the magnetic properties of ferrous and cobaltous oxides.1. Prog. Theoret. Phys. 17, 177196 (1957).15. Kanamori, J. Theory of the magnetic properties of ferrous and cobaltous oxides.2. Prog. Theoret. Phys. 17, 197222 (1957).16. Goodenough, J. B. An interpretation of the magnetic properties of the perovskite-type mixed crystals La1 xSrxCoO3. J. Phys. Chem. Solids 6, 287297 (1958).

17. Lu, X. Z., Whangbo, M. H., Dong, S., Gong, G. & Xiang, J. Giant ferroelectric polarisation of CaMn7O12 induced by a combined effect of Dzyaloshinskii-

Moriya interaction and exchange striction. Phys. Rev. Lett. 108, 187204 (2012).18. Dzyaloshinskii, I. Theory of helicoidal structures in antiferromagnets. 1. nonmetals. 1960-1971. Sov. Phys. JETP 19, 960971 (1964).

19. Moriya, T. Anisotropic superexchange interactions and weak ferromagnetism. Phys. Rev. 120, 9198 (1960).

Acknowledgements

The work done at the University of Oxford was funded by an EPSRC grant number EP/J003557/1, entitled New Concepts in Multiferroics and Magnetoelectrics.

Author contributions

N.P. conducted the X-ray diffraction experiment and performed the data analysis, samples were grown by C.M., and the manuscript was written by N.P., R.D.J, and P.G.R. The project direction and ideas underpinning the results were conceptualized and developed by R.D.J., L.C.C. and P.G.R.

Additional information

Supplementary Information accompanies this paper on http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Perks, N.J., et al. Magneto-orbital helices as a route to coupling magnetism and ferroelectricity in multiferroic CaMn7O12. Nat. Commun.

3:1277 doi: 10.1038/ncomms2294 (2012).

6 NATURE COMMUNICATIONS | 3:1277 | DOI: 10.1038/ncomms2294 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2012 Macmillan Publishers Limited. All rights reserved.