# Magnetically induced ferroelectricity in CuMnSnS and CuMnSnSe

###### Abstract

We investigate magnetically-induced ferroelectricity in CuMnSnS by means of Landau theory of phase transitions and of ab initio density functional theory. As expected from the Landau approach, ab initio calculations show that a non-zero ferroelectric polarization along the direction (of the order of a tenth of ) is induced by the peculiar antiferromagnetic configuration of Mn spins occurring in CuMnSnS. The comparison between , calculated either via density-functional-theory or according to Landau approach, clearly shows that ferroelectricity is mainly driven by Heisenberg-exchange terms and only to a minor extent by relativistic terms. At variance with previous examples of collinear antiferromagnets with magnetically-induced ferroelectricity (such as AFM-E HoMnO), the ionic displacements occurring upon magnetic ordering are very small, so that the exchange-striction mechanism (i.e. displacement of ions so as to minimize the magnetic coupling energy) is not effective here. Rather, the microscopic mechanism at the basis of polarization has mostly an electronic origin. In this framework, we propose the small magnetic moment at Cu sites induced by neighboring Mn magnetic moments to play a relevant role in inducing . Finally, we investigate the effect of the anion by comparing CuMnSnSe and CuMnSnS: Se-4 states, more delocalized compared to S-3 states, are able to better mediate the Mn-Mn interaction, in turn leading to a higher ferroelectric polarization in the Se-based compound.

###### pacs:

75.50.E, 77.80.-e, 75.85.+t##
I Introduction

Multiferroic materials have attracted both scientific and industrial interests, due to their profound physics and novel functionalities, respectively.Cheong () Among them, many studies have been recently performed for “improper” ferroelectric materials,Kimura (); Hur (); Picozzi (); Sergienko (); Yamauchi () where ferroelectricity is driven by non-centrosymmetric spin, charge, or orbital ordering. In this case, a stronger magnetoelectric coupling is expected, compared to the conventional proper ferroelectrics. So far, mainly manganites (such as TbMnO, TbMnO, and HoMnOKimura (); Hur (); Picozzi (); Sergienko ()) and other perovskite oxidesGianluca1 (); Sanjeev () have been carefully investigated by means of theoretical and experimental approaches, in order to better understand cross-coupling effects between spin/orbital orderings and dielectric properties.

The search for novel multiferroics other than oxides has already started: recently Nénert et al.Nenert () have suggested the ternary copper chalcogenide CuMnSnS (with nominal valences as , , , and ) as a candidate of “improper” ferroelectricity, using a symmetry-based analysis and the Landau phenomenological theory, with the ferroelectricity induced by a peculiar antiferromagnetic (AFM) Mn spin configuration. Experimentally, neutron diffraction measurements have clarified that such AFM configuration occurs under =8.8 K, with Mn-spins slightly deviating from the axis by an angle between 6 and 16 degree. Moreover, the magnetization curve has suggested the spin quantization axis to be close to the direction.Fries () Additionally, the small dependence of susceptibility upon the direction of magnetic field implied a rather small magnetic anisotropy.

Similar to the famous stannite, CuFeSnS,Caneschi () CuMnSnS crystallizes in the tetragonal sphalerite superstructure, which shows the (No.121) space group with symmetry operations {, , , , , , , }.symmetry () The crystal structure is similar to zincblende, often seen in III-V and II-VI type semiconductors, where each cation is tetrahedrally coordinated to anions. In this respect, multiferroic effects in this class of materials might prove useful in the context of semiconductor-based spintronics. Although the sphalerite crystal structure itself lacks inversion symmetry, this does not automatically imply that the crystal is polar: due to the mirror and -rotation symmetries, the ferroelectric polarization is forbidden along any direction. On the other hand, Nénert et al. have claimed that the experimentally observed AFM spin order breaks some symmetries, so that arises along the direction (i.e. the crystallographic direction) with possibly large magnetoelectric coupling.Nenert () The Landau phenomenological theory, based on the symmetry analysis of both the crystal structure and the magnetic order, is a powerful tool to investigate magnetically-induced ferroelectricity.Sergienko (); Toledano () However, in order to quantify the value of and to clarify the relation between magnetism and ferroelectricity via the electronic structure, a combination of Landau theory and ab initio density functional theory (DFT) calculations is desirable.Picozzi ()

In this paper, we investigate magnetically-induced ferroelectricity in CuMnSnS by using both Landau theory and DFT in the aim of carefully understanding the microscopic mechanism at the basis of the multiferroic behaviour. To check the effect of the anions on , we perform similar simulations by substituting S with Se, for CuMnSnSe. First, we will discuss the direction and the value of obtained by DFT, confirming the predictions based on Landau theory; afterwards, we will focus on the ferroelectric switching, in terms of different Mn-spin configurations along the adiabatic path between positive and negative ferroelectric states. Finally, the microscopic mechanism leading to ferroelectricity via spatially-polar Cu and Mn charge densities caused by Heisenberg-exchange is discussed.

## Ii Methodology and structural details

DFT calculations are performed for CuMnSnS and CuMnSnSe by using the “Vienna Ab initio Simulation Package” (VASP), based on Projector Augmented Wave (PAW) pseudopotentials.Kresse () The Perdew-Becke-Erzenhof (PBE) approach to the generalized gradient approximation (GGA) is employed for the exchange-correlation potential.Perdew () The plane wave energy cut-off is 400 eV for each atom. During structural optimization, a threshold on the atomic forces is set as 0.01 eV/. Internal atomic coordinates are optimized starting from experimental data of CuMnSnSFries () and CuMnSnSe.Sachanyuk () The Monkhost-Pack -point grid in the Brillouin zone is used.Monkhorst () Cu-, Mn-, Sn-, S-, and Se- electrons are treated as valence states. The GGA+ calculations within the Dudarev’s approachDudarev () are performed by applying a Hubbard-like potential only for Mn- states. An on-site Coulomb parameter = 4.0 eV and an exchange parameter = 0.89 eV are used.

Lattice constants are fixed to the experimental data: Å, for CuMnSnSFries () and Å, for CuMnSnSeSachanyuk () in the tetragonal structure. In order to impose the experimentally observed AFM spin configuration in CuMnSnSFries () with the propagation vector , we build the magnetic super-cell with 8 formula/units (i.e. 64 atoms/cell, see Fig.1(b) and (c)) (note that, to our knowledge, the experimental magnetic configuration for CuMnSnSe has not been reported yet). Since the observed AFM configuration breaks the symmetries (except for ), ferroelectric polarization is allowed to be magnetically-induced along the direction. To evaluate , we also perform calculations for the ferromagnetic (FM) spin configuration, which is considered as a reference non-polar state: the ferroelectric polarization is then calculated as . To connect paraelectric (PE) and ferroelectric (FE) states, calculations for noncollinear spin configurations (see below) are also performed according to Ref.Hobbs (). The Berry phase approach developed by King-Smith and VanderbiltVanderbilt () is employed to calculate the electric polarization , where we integrate over eight -point strings parallel to .

## Iii Magnetically induced ferroelectric polarization

### iii.1 Landau theory

In the Landau potential terms, the ferroelectric polarization appears as coupled with magnetic order parameters, in terms which are invariant under symmetry operations.Mostovoy (); Harris () Nenért et al.Nenert () have built up the free energy equation by considering two Mn spins ( with position (0,0,0) and with position (,,) in the nuclear unit cell) and set up two magnetic order parameters, e.g., the FM order parameter, , and the AFM order parameter, . From the free energy, they have showed that only the component is coupled with both and , and also with . Thereby, they have concluded that a spontaneous polarization arises along the direction upon the AFM configuration. Here we extend their model to fit into our AFM supercell with eight Mn spins, so that we can discuss the adiabatic switching path between the positive FE () and negative FE () states by changing the direction of Mn spins and the consequent electric polarization. (The same procedure was applied for HoMnO with E-type AFM configuration.Picozzi (); Sergienko ()) With these eight Mn spins (labelled as , as in Fig.1), two AFM order parameters are constructed as

(1) |

corresponding to the positive and negative ferroelectric domain phases with and , respectively. In order to switch between the two phases (as seen in the following section), four body-centered Mn spins (, , , ) must be flipped. Taking into account the irreducible corepresentation of the little group of with the propagation vector , one can write the possible magnetoelectric coupling terms in the second degree of magnetic order parameters and the dielectric term in the Landau potential as:

(2) | |||||

where the coefficients denote phenomenological parameters and , and are the components of the dielectric susceptibility. Since Mn magnetic moments have been experimentally observed in the plane, the terms containing and are omitted in the above equation. The first and second terms denote exchange terms, the last term is the dielectric energy and the remaining terms are relativistic terms. The minimization of Eq.(2) gives a finite polarization along :

(3) | |||||

whereas non-zero and are not obtained. This result confirms Nénert’s predictions: a spontaneous is evidently induced by the AFM coupling. However, at this point we cannot conclude whether the Heisenberg-exchange- or the relativistic-interaction terms is dominant in inducing .

### iii.2 DFT results

Figure 2 shows the partial density of states (DOS) projected onto (a) Cu-, (b) Mn- and (c) S- states in AFM CuMnSnS and onto (d) Cu-, (e) Mn- and (f) Se- states in AFM CuMnSnSe. To discuss the electronic structure, we focus on CuMnSnS. The valence bands just below the Fermi level mainly come from Cu-3 states with small contribution from Mn-3 states, while the bottom of conduction bands is mainly from Mn-3 states. We recall that, since Mn atoms are tetrahedrally surrounded by S atoms, the five-fold degenerate 3 states split into the three-fold states and two-fold states due to the crystal electric field (CEF). Additionally, the three-fold degenerated states form bonding and anti-bonding states. The majority spin states are fully occupied whereas the minority spin states are fully unoccupied, corresponding to a high-spin configuration (with 4.4 magnetic moment inside the atomic sphere with 1.3 Å radius). Cu-3 states are weakly spin polarized with a small induced magnetic moment (equal to 0.015 inside the atomic sphere with 1.3 Å radius). We note, in fact, that the distance between Mn atomic sites is rather large, so that localized Mn-spins interact via Cu and S sites. Fig.2 (a) and (b) show that there is a small hybridization between Cu- and Mn- states, so that Cu-ions are slightly spin-polarized. As we’ll show later, the small Cu spins caused by the Cu-Mn hybridization play an important role for ferroelectricity in these materials.

(eV) | SOC | ||||
---|---|---|---|---|---|

CuMnSnS | 0 | on, | 0 | 0.209 | 0 |

on, //z | 0 | 0.208 | 0 | ||

off | 0 | 0.209 | 0 | ||

4.0 | on, //x | 0 | 0.054 | 0 | |

on, //z | 0 | 0.054 | 0 | ||

off | 0 | 0.055 | 0 | ||

CuMnSnSe | 0 | on, | 0 | 0.339 | 0 |

on, //z | 0 | 0.335 | 0 | ||

off | 0 | 0.339 | 0 | ||

4.0 | on, //x | 0 | 0.092 | 0 | |

on, //z | 0 | 0.092 | 0 | ||

off | 0 | 0.095 | 0 |

Table 1 shows the DFT calculated in the optimized structure, imposing the AFM configuration. Consistent with Landau theory results, only has a finite value. Both compounds show a rather small compared to other materials with magnetically-induced ferroelectric polarization (comparing with equal to few in E-type AFM HoMnOPicozzi (), one or two orders of magnitude smaller); this is likely due to the weaker “indirect” interaction (via Mn-S(Se)-Cu-S(Se)-Mn bond).

The reason why CuMnSnSe shows a larger than CuMnSnS has to be traced back to the fact that Mn- states can hybridize more with Se- states than with (more localized) S- states, in turn mediating a stronger interaction between Mn- and Cu- states. In the same respect, we note that the difference between = 0 and 4.0 eV as far as is concerned can be equivalently understood in terms of Mn-3 and Cu-3 hybridization: the on-site Coulomb interaction pushes Mn-3 states deeper in energy (clearly shown in Fig.2) and decreases the Mn-3 and Cu-3 hybridization. Therefore, values at = 4.0 eV are smaller compared to the cases of = 0 eV.

The important finding here is that the contribution to due to spin-orbit coupling (SOC) is very small, as shown by the comparison of polarization values switching on or off the relativistic interactions (cf. Table 1). In this respect, we also note that the total energy decreases - by 0.04 meV/Mn - when the spin direction is changed from the to axis, in agreement with experiments.

We recall that our Landau theory analysis didn’t clarify which magnetic term ( Heisenberg exchange or relativistic) mainly contributes to ; on the other hand, DFT results unambiguously show the relativistic term to give only a minor contribution. Therefore we can omit the terms depending on either or in Eq. (3), but keep the term so that the equation is simplified as:

(4) | |||||

Here, depends only on the inner products of Mn spin vectors, therefore is described by the relative spin angle. terms in the above equation are neglected because these terms disapper when calculating the polarization difference between FE and PE states. Progressive rotation of (, , , ) Mn spins can switch the ferroelectric polarization between and continuously. Assuming the non-collinear (NC) spin configuration as and , Eq.(2) becomes

(5) |

Indeed, the trend is obtained when we perform DFT calculations on the NC configuration by varying between 0 and 180 with fully optimized atomic coordinates for each spin configuration. Figure 3 shows the calculated and the total energy difference along the adiabatic path as a function of the relative Mn-spin angle at CuMnSnSe. We also show , calculated with atomic coordinates fixed to the case of , so that there is no contribution to the electric polarization from atomic displacements. The similarity between polarization values for Fig.3 (b) and (c) - including or not the ionic displacements - highlights stannite as a paradigmatic example of purely electronic ferroelectricity (see also discussion below).

As expected from Landau theory and as pointed out above, the DFT-calculated is fitted by a cosine curve, with coefficients 0.0053 and 0.0015 at = 0 and 4.0 eV, respectively. The function shows convex and symmetrical behavior between positive and negative FE states. With increasing from (positive FE state), increases and reaches maximum value at (PE state) with the energy barrier 0.68 and 0.25 meV/Mn at = 0 and 4.0 eV, respectively. This energy barrier —coming from combination of term and -related term— is much lower than in the HoMnO case, where = 8 meV/Mn at = 0 eV.Picozzi () This is fully consistent with the weak Mn-Mn exchange-coupling constants involved.

## Iv Mechanism of ferroelectricity

According to the previous section, ferroelectric polarization is induced by Heisenberg-exchange coupling between Mn spins. We also noted that the Heisenberg exchange between collinear spins is more efficient in driving ferroelectric polarization than any spin-orbit Dzyaloshinskii-Moriya interactionMostovoy (): the order of magnitude of here (tenth or hundredth of ) is about one order of magnitude higher than that experimentally observed and calculated (hundredth or thousandths of ) in spiral-based systems such as TbMnO. The origin of ferroelectricity could be explained by analogy with the HoMnO case, where O ions move in order to reduce the magnetic coupling energy between Mn spins. Such “exchange-striction” mechanism has been referred to as “inverse Goodenough-Kanamori” (iGK) interaction.Picozzi (); Yamauchi2 () However, since here in CuMnSnS the Mn ions are far apart, the Mn spin interaction is too “indirect” and it is difficult (and non-intuitive) to guess how S(Se) ions would move to reduce the magnetic coupling energy. Indeed, as shown in Fig.3 (b) and (c), and have almost the same value: the ionic displacements are rather small and do not contribute much to .

Rather, to explain the “electronic” origin of the Heisenberg-driven polarization, we assume small-size spins (actually calculated as 0.015 ) on Cu sites, which interact with Mn spins via “Cu-S(Se)-Mn” bonds. Consider S(Se) sites: from the chemical point of view, all of them are equivalent, since they are tetrahedrally surrounded by one Mn, two Cu and one non-magnetic Sn atoms in the symmetry . However, from the magnetic point of view, due to the small magnetic moments at Cu sites, S(Se) ions can have two different magnetic environments at the FE state (): the magnetic moments of the surrounding Cu atoms can align either parallel or antiparallel, as shown in Fig.4 (a) and (b). The anions are therefore split into two categories, they become inequivalent upon this peculiar AFM spin ordering, breaking the symmetries and paving the way to ferroelectric polarization. Since Mn spins interact via such magnetically inequivalent S(Se) ions, the charge densities at Mn sites are modulated by Heisenberg-exchange to induce . Similarly, charge densities at Cu sites with small magnetic moments caused by Mn- and Cu- hybridization can polarize. This is clearly observed in Fig.5, where we show the difference between FE and PE states: Mn and Cu charge densities are clearly non-symmetric along , i.e. along is induced. The results are consistent with the magnetically broken point group 2 at Mn and Cu sites (the point groups at Cu and Mn sites were “originally” - i.e. without magnetic ordering - and , respectively). The “original” site symmetry at S(Se) sites is , which is already a polar point group; therefore, a magnetically induced polar charge does not appear. We note that, in the PE state (), some Cu magnetic moments point along the direction, corresponding to the NC Mn spin configuration. All S(Se) ions have then an equivalent magnetic environment, where half of the Cu magnetic moments align parallel to Mn and the other half perpendicular to it (Fig.4 (c)); in this configuration, is therefore not expected along any direction.

## V Summary

We have investigated magnetically-induced ferroelectricity in AFM CuMnSnS and CuMnSnSe by means of Landau theory and DFT calculations. The combination of both methods clearly highlights that the driving force of is not the relativistic term, but the Heisenberg exchange term. Microscopically, small localized Cu magnetic moments induced by Mn magnetic moments play an important role in inducing . The Mn-Mn spin interaction via the magnetically inequivalent S(Se) ions produces a polar modulation of the charge density on Mn and Cu sites, finally inducing a net with predominantly electronic contribution and negligible contribution from ionic displacements. We therefore show that AFM CuMnSnS and CuMnSnSe are remarkable examples of novel magnetically-induced ferroelectric materials, with improper multiferroicity branching out into non-oxide and non-octahedral based systems.

###### Acknowledgements.

The research leading to these results has received funding from the European Research Council under the European Community’s 7th Framework Programme (FP7/2007-2013) / ERC grant agreement n. 203523. Computational support by CASPUR supercomputing center (Rome) is acknowledged. TF thanks G. Nénert for helpful discussions and for his careful reading of the manuscript. The crystal structures and charge densities in this paper are plotted by using a software VESTA.VESTA ()## References

- (1) S.-W. Cheong and M. Mostvoy, Nat. Mater. 6 (2007) 13.
- (2) T. Kimura, T. Goto, H. Shintani, K. Ishizaka, T. Arima, and Y. Tokura, Nature (London) 426 (2003) 56.
- (3) N. Hur, S. Park, P. A. Sharma, J. S. Ahn, S. Guha, and S. W. Cheong, Nature (London) 429 (2004) 392.
- (4) S. Picozzi, K. Yamauchi, B. Sanyal, I. A. Segienko, and E. Dagotto, Phys. Rev. Lett. 99 (2007) 227201.
- (5) I. A. Sergienko, C. Sen, and E. Dagotto, Phys. Rev. Lett. 97 (2006) 227204.
- (6) K. Yamauchi, F. Freimuth, S. Blügel, and S. Picozzi, Phys. Rev. B 78 (2008) 014403.
- (7) G. Giovannetti, S. Kumar, D. Khomskii, S. Picozzi, and J. van den Brink, Phys. Rev. Lett. 103 (2009) 156401.
- (8) S. Kumar, G. Giovannetti, J. van den Brink, and S. Picozzi, cond-mat/9091439 (unpublished).
- (9) G. Nénert and T. T. M. Palstra, J. Phys.: Condens. Matter 21 (2009) 176002.
- (10) T. Fries, Y. Shapira, F. Palacio, M. C. Morón, G. J. McIntyre, R. Kershaw, A. Wold, and E. J. McNiff Jr., Phys. Rev. B 56 (1997) 5624.
- (11) A. Caneschi, C. Cipriani, F. Di Benedetto, and R. Sessoli, Phys. Chem. Minerals 31 (2004) 190.
- (12) and denote mirror symmetries with respect to (110) and (1-10) planes, respectively.
- (13) P. Tolédano, Ferroelectrics 161 (1994) 257.
- (14) G. Kresse and J. Furthmuller, Phys. Rev. B 54 (1996) 11169.
- (15) J. P. Perdew, K. Burke, and M. Ernzerhol, Phys. Rev. Lett. 77 (1996) 3865.
- (16) V. P. Sachanyuk, I. D. Olekseyuk, and O. V. Parasyuk, Phys. Stat. Sol. (a) 203 (2006) 459.
- (17) H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13 (1976) 5188.
- (18) S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Jumphreys, and A. P. Sutton, Phys. Rev. B 57 (1998) 1505.
- (19) D. Hobbs, G. Kresse, and J. Hafner, Phys. Rev. B 62 (2000) 11556.
- (20) R. D. King-Smith and D. Vanderbilt, Phys. Rev. B 47 (1993) 1651.
- (21) M. Mostovoy, Phys. Rev. Lett. 96 (2006) 067601.
- (22) A. B. Harris, T. Yildirim, A. Aharony, and O. Entin-Wohlman, Phys. Rev. B 73 (2006) 184433.
- (23) K. Yamauchi and S. Picozzi, J. Phys.: Condens. Matter 21 (2009) 064203.
- (24) K. Momma and F. Izumi, J. Appl. Crystallogr. 41 (2008) 653.