Outline
Comptes Rendus

Molecular magnets explored by density functional theory calculations
Comptes Rendus. Chimie, Corrélations de spin et de charge dans les matériaux moléculaires - Propriétés physiques, chimie et aspects « matériaux », Volume 10 (2007) no. 1-2, pp. 65-67.

Abstract

After a short introduction in the field of molecular magnetism, we focus on calculations of the magnetic anisotropy energy based on density functional theory. The results obtained make one confident in the predictive power of the formalism.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crci.2006.07.009
Keywords: Molecular magnetism, Nanomagnets, DFT magnetic anisotropy

Jens Kortus 1

1 Institut für Theoretische Physik, TU Bergakademie Freiberg, Leipziger Str. 23, 09599 Freiberg, Germany
@article{CRCHIM_2007__10_1-2_65_0,
     author = {Jens Kortus},
     title = {Molecular magnets explored by density functional theory calculations},
     journal = {Comptes Rendus. Chimie},
     pages = {65--67},
     publisher = {Elsevier},
     volume = {10},
     number = {1-2},
     year = {2007},
     doi = {10.1016/j.crci.2006.07.009},
     language = {en},
}
TY  - JOUR
AU  - Jens Kortus
TI  - Molecular magnets explored by density functional theory calculations
JO  - Comptes Rendus. Chimie
PY  - 2007
SP  - 65
EP  - 67
VL  - 10
IS  - 1-2
PB  - Elsevier
DO  - 10.1016/j.crci.2006.07.009
LA  - en
ID  - CRCHIM_2007__10_1-2_65_0
ER  - 
%0 Journal Article
%A Jens Kortus
%T Molecular magnets explored by density functional theory calculations
%J Comptes Rendus. Chimie
%D 2007
%P 65-67
%V 10
%N 1-2
%I Elsevier
%R 10.1016/j.crci.2006.07.009
%G en
%F CRCHIM_2007__10_1-2_65_0
Jens Kortus. Molecular magnets explored by density functional theory calculations. Comptes Rendus. Chimie, Corrélations de spin et de charge dans les matériaux moléculaires - Propriétés physiques, chimie et aspects « matériaux », Volume 10 (2007) no. 1-2, pp. 65-67. doi : 10.1016/j.crci.2006.07.009. https://comptes-rendus.academie-sciences.fr/chimie/articles/10.1016/j.crci.2006.07.009/

Original version of the full text

1 Introduction

“Molecular magnetism deals with magnetic properties of isolated molecules and/or assemblies of molecules” [1], as expressed by Olivier Kahn, one of the founders of this research topic. This definition is quite general and there has been recently more emphasis on the aspect of the rational design of molecular magnetic properties in the field [2]. Therefore, molecular magnetism is seen “as a discipline which conceives, realizes, studies, and uses new molecular materials bearing new but predictable magnetic (and other) physical property” [3]. Research in molecular magnetism is driven by inter-disciplinary interaction between various fields in chemistry, physics, engineering and even biology.

This field of research is very attractive for first-principles microscopic simulations, because the crystal structure of molecular magnets is well defined, reproducible and is made available rapidly. The progress in the field is clearly driven by advances in chemical synthesis of the materials and ingenious experiments. However, the combined efforts of physicists and material scientists, particularly theorists, inspire confidence that theory may become one day useful not only for explaining and postdiction, but also for guiding the synthesis of new promising materials. As for now, the theory dealing with first-principles calculations tries to keep pace with experiment, still trying to reproduce it rather than to lead.

Nevertheless, there has been a notable progress in the prediction of exchange interactions and magnetic anisotropy energies from density functional theory during the last few years. The calculation of exchange couplings seems nowadays possible and feasible even for larger systems using methods based on density functional theory (DFT). For a recent review on the topic of calculation of exchange couplings using DFT see e.g. [4]. In contrast to cases where magnetic exchange interactions follow the famous Goodenough–Kanamori rules [5–7], we still have to await similarly clean general insights derived from theory in the case of magnetic anisotropy, which is another crucial property of molecular magnets, which could revolutionize the rational design of molecular nanomagnets.

2 Mn12-acetate: the most famous example

The discovery of the molecule Mn12O12(CH3COO)16(H2O)4 with a magnetic ground state of S = 10 showing a magnetic hysteresis [8] due to the properties of the single Mn12-molecule has created large interest in the field. The observed hysteresis in molecular magnets is not due to re-orientation of domains as in conventional ferromagnets. Instead, it is an intrinsic molecular property that results from the magneto-crystalline anisotropy. This process can only be observed because the relaxation time is very large compared to the measurement time. The relaxation of the magnetization becomes indeed very slow at low temperatures (of the order of several months at 2 K). Therefore, by applying a magnetic field the material can be magnetized and will keep this magnetization for a long time if kept at very low temperatures. The so-called single-molecule magnets (SMM), which are also called molecular nanomagnets, can often be crystallized and form molecular crystals with relatively large distance between the molecular constituents. Therefore interactions between the molecular entities remain weak, so that the magnetic behavior probed by experiments is often dominated by intra-molecular effects.

Mn12-ac has been synthesized and reported first in 1980 [9]. The molecular crystal has tetragonal symmetry with space group I4¯, a single Mn12-acetate cluster in the crystal posses an S4 symmetry. There are now several modifications of Mn12-acetate known, with different crystal structures, solvent molecules and water coordination – see Ref. [10] for more information.

3 Spin-orbit coupling and magnetic anisotropy energy

Strong spin-orbit coupling would try to align the spins of electrons in the magnetic d shell with the magnetic field. In a chemical environment, the d shell experiences a crystal field splitting due to the electric field of the neighboring atoms. The energy cost or gain of aligning the electron spins in an external magnetic field will then depend on the spatial arrangement of the d shell. In two special cases, however, the magnetic anisotropy can be suppressed. If the crystal field splitting (i.e. anisotropic interaction with the ligands) is small, the spins in the d shell may freely rotate and will follow the external magnetic field, without much loss in energy. On the other hand, if the crystal field splitting is much larger than the spin-orbit coupling, the system (transition metal ion + ligands) prefers a configuration with zero (quenched) orbital moment and no magnetic anisotropy. In a semiclassical picture this orbital quenching can be interpreted, that the orbital momentum processes in the crystal field, so that its magnitude remains unchanged but all its components average to zero. However, the spin-orbit coupling cannot be completely ignored and it will mix in states with non-zero orbital momentum. This is responsible for the deviation of the g-factor from the pure spin value of 2.

Therefore, the height of the magnetic anisotropy barrier emerges by a quite subtle interplay between the strength of the spin-orbit coupling and the interaction with the crystal field due to the ligands. The temperature up to which a magnetic device will retain its preferential magnetic orientation is therefore determined by the magnetic anisotropy, hence high anisotropy is one of major priorities in the design of useful SMM devices.

Recently, Pederson and Khanna [11,12] have developed a method for accounting for second-order anisotropy energies. This method relies on a simple albeit exact method for spin-orbit coupling and a second-order perturbative treatment of the spin Hamiltonian to determine the dependence of the total energy on spin projection. It makes use of the Cartesian representation of the spin-orbit term, which is exact and also is more adaptable for multi-center systems.

4 Magnetic anisotropy in single-molecule magnets

The spin Hamiltonian describing the magnetic anisotropy up to second order can be expressed as:

H=DSz2+E(Sx2Sy2).(1)
We distinguish between axial and transverse anisotropies, with their corresponding parameters D and E.

The values of the axial anisotropy D available from a number of experiments for different SMM, and for several SMM first-principle calculations carried out with the use of the NRLMOL code are summarized in Table 1. The NRLMOL program package developed by Pederson et al. is an all-electron Gaussian-type orbital implementation of DFT [13].

Table 1

Comparison of the calculated (DFT-NRLMOL) and experimental magnetic anisotropy parameter D for selected single-molecule magnets

MoleculeSD(K)
TheoryExperiment
Mn12O12(O2CH)16(H2O)410−0.56a−0.56b
[Fe8O2(OH)12(C6H15N3)6Br6]2+10−0.53c−0.30d
[Mn10O4(2,2′-biphenoxide)4Br12]4–13−0.06e−0.05f
Co4(CH2C5H4N)4(CH3OH)4Acl46−0.64g−0.7 to −0.9h
Fe4(OCH2)6(C4H9ON)65−0.56i−0.57j
Mn9O34C32N3H3517/2−0.33−0.36l
Ni4O16C16H404−0.385−0.40l
Mn4O3Cl4(O2CCH2CH3)3(NC5H5)39/2−0.58m−0.72n

a [11,14];

b [15,16];

c [17];

d [18];

e [19];

f [20];

g [21];

h [22];

i [30];

j [23];

l [24];

m [25];

n [26].

In all the cases presented here the calculated spin ordering is in agreement with experiment. The calculated D parameters for Mn12, Mn10, Mn9, and the ferric star Fe4 are in excellent agreement with experimental values. The only remarkable discrepancy is found for Fe8, a system which seems to pose complications for the DFT treatment. Apparently the DFT may be unable to predict the ground-state density accurately enough due to important electronic correlations beyond the mean-field treatment. Similar deviations between experiment and theory have been found for a series of monomers, where D is underestimated by about a factor of two [27].

The SMM listed in Table 1 are in general characterized by a high-spin ground state. However, a high spin state does not necessarily correlate with a high anisotropy barrier. The prefactor D is also very important. In order to increase the barrier one has to understand and control D, which will be the main goal of future research in this area. In all cases where the E parameter is not zero by symmetry it has been predicted with similar accuracy as D – see relevant references for details.

The results obtained make one confident in the predictive power of the formalism. A microscopic understanding (based on the electronic structure of SMM) of the magnetic anisotropy parameters seems to be crucial for the rational design of single-molecule magnets.

Acknowledgments

The author thanks and appreciates useful discussions with A. Postnikov, Mark R. Pederson, Roberta Sessoli, Jürgen Schnack, Eliseo Ruiz, Carlo Massobrio, Marc Drillon and Stefan Blügel. This contribution is partially based on similar work in collaboration with A. V. Postnikov [28,29].


References

[1] O. Kahn Molecular Magnetism, John Wiley & Sons, Singapore, 1993

[2] M. Verdaguer Polyhedron, 20 (2001), p. 1115

[3] W. Linert, M. Verdaguer (Eds.), Molecular Magnets, Monatsh. Chem./Chem. Mon. 134 (2) (2003) (special issue).

[4] E. Ruiz, Structure and Bonding, vol. 113, Springer-Verlag, Berlin, Heidelberg, 2004

[5] J.B. Goodenough Phys. Rev., 100 (1955), p. 564

[6] J.B. Goodenough J. Phys. Chem. Solids, 6 (1958), p. 287

[7] J. Kanamori J. Phys. Chem. Solids, 10 (1959), p. 87

[8] R. Sessoli; H. Tsai; A. Schake; S.Y. Wang; J.B. Vincent; K. Folting; D. Gatteschi; G. Christou; D.N. Hendrickson J. Am. Chem. Soc., 115 (1993), p. 1804

[9] T. Lis Acta Crystallogr. Soc. B, 36 (1980), p. 2042

[10] D. Gatteschi; R. Sessoli Angew. Chem., Int. Ed. Engl., 42 (2003), p. 268

[11] M.R. Pederson; S.N. Khanna Phys. Rev. B, 59 (1999), p. 693 (R)

[12] M.R. Pederson; S.N. Khanna Phys. Rev. B, 60 (1999), p. 9566

[13] M.R. Pederson; D.V. Porezag; J. Kortus; D.C. Patton Phys. Stat. Sol. (b), 217 (2000), p. 197

[14] M.R. Pederson; S.N. Khanna Chem. Phys. Lett., 307 (1999), p. 253

[15] K.M. Mertes; Y. Suzuki; M.P. Sarachik; Y. Paltiel; H. Shtrikman; E. Zeldov; E. Rumberger; D.N. Hendrickson; G. Christou Phys. Rev. Lett., 87 (2001), p. 227205

[16] A.L. Barra; D. Gatteschi; R. Sessoli Phys. Rev. B, 56 (1997), p. 8192

[17] J. Kortus; C.S. Hellberg; M.R. Pederson; S.N. Khanna Eur. Phys. J. D, 16 (2001), p. 177

[18] M. Dressel; B. Gorshunov; K. Rajagopal; S. Vongtragool; A.A. Mukhin Phys. Rev. B, 67 (2003), p. 060405

[19] J. Kortus; T. Baruah; N. Bernstein; M.R. Pederson Phys. Rev. B, 66 (2002), p. 092403

[20] A.L. Barra; A. Caneschi; D. Gatteschi; D.P. Goldberg; R. Sessoli J. Solid-State Chem., 145 (1999), p. 484

[21] T. Baruah; M.R. Pederson Chem. Phys. Lett., 360 (2002), p. 144

[22] M. Murrie; S.J. Teat; H. Stockli-Evans; H.U. Güdel Angew. Chem., Int. Ed. Engl., 42 (2003), p. 4653

[23] S. Schromm, O. Waldmann, P. Müller, in press.

[24] S. Piligkos; G. Rajaraman; M. Soler; N. Kirchner; J. van Slageren; R. Bircher; S. Parsons; H.U. Gudel; J. Kortus; W. Wernsdorfer; G. Christou; E.K. Brechin J. Am. Chem. Soc., 127 (2005), p. 5572

[25] K. Park; M.R. Pederson; S.L. Richardson; N. Aliaga-Alcalde; G. Christou Phys. Rev. B, 68 (2003), p. 020405

[26] W. Wernsdorfer; N. Aliaga-Alcalde; D.N. Hendrickson; G. Christou Nature, 416 (2002), p. 406

[27] J. Cirera, E. Ruiz, J. Kortus, (in press).

[28] A.V. Postnikov, J. Kortus, M.R. Pederson, Physica Status solidi (b) 243 (2006) 2533.

[29] J. Kortus, A.V. Postnikov, in: M. Rieth, W. Schommers (Eds.), Handbook of Theoretical and Computational Nanotechnology, American Scientific Publishers, 2006, ISBN: 1-58883-042-X.

[30] J. Kortus; M.R. Pederson; T. Baruah; N. Bernstein; C.S. Hellberg Polyhedron, 22 (2004), p. 1871


Comments - Policy