Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

High temperature ferromagnetism in π-conjugated two-dimensional metal–organic frameworks

Wenbin Li a, Lei Sun b, Jingshan Qi c, Pablo Jarillo-Herrero d, Mircea Dincă *b and Ju Li *c
aResearch Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
bDepartment of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA. E-mail:
cDepartment of Nuclear Science and Engineering, Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA. E-mail:
dDepartment of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA

Received 17th November 2016 , Accepted 20th January 2017

First published on 8th February 2017

We use first-principles calculations to show that the square symmetry of two-dimensional (2D) metal–organic frameworks (MOFs) made from octaamino-substituted phthalocyanines and square planar Ni2+ ions, which enable strong conjugation of π electrons, has a critical impact on the magnetic properties of the lattice. In particular, we predict the unexpected emergence of a rare high-temperature ferromagnetic half-metallic ground state in one case. Among charge neutral MOFs made from (2,3,9,10,16,17,23,24)-octaiminophthalocyanine (OIPc) metallated with divalent first-row transition metal ions (M-OIPc; M = Cr2+, Mn2+, Fe2+, Co2+, Ni2+, Cu2+, Zn2+) and connected through square planar Ni-bisphenylenediimine moieties, NiMn-OIPc exhibits a half-metallic and ferromagnetic ground state with a large exchange energy resulting from the unique strong hybridization between the d/π orbitals of Mn, the Pc ring, and the Ni-bisphenylenediimine nodes. Notably, we show that for NiMn-OIPc there is a considerable difference between the ferromagnetic ordering temperature (Tc) predicted by a 2D Ising model, which exceeds 600 K, and a Tc of 170 K predicted by our more realistic Monte Carlo simulation that includes magnetic anisotropy. Critically, our simulations adopt two spin models that incorporate magnetic anisotropy in the form of exchange anisotropy and single-ion anisotropy. We further show that in the bulk, 2D layers of NiMn-OIPc adopt a slipped-parallel stacking configuration, and exhibit interlayer magnetic coupling that is sensitive to the relative in-plane displacement between adjacent layers. These results highlight the critical role of magnetic anisotropy in modeling the properties of 2D magnetic systems. More generally, it demonstrates that strong hybridization between open-shell ions and delocalized aromatic π systems with appropriate symmetry, combined with large magnetic anisotropy, will be an effective design strategy to realize ferromagnetic 2D MOFs with high Tc.


A central focus of current research at the nanoscale is two-dimensional (2D) materials – planar sheets of materials with atomic thickness.1 2D materials such as graphene and transition metal (TM) dichalcogenides exhibit fundamentally new properties that could be essential for next-generation electronic, photonic and renewable-energy applications.2–6 Recently, 2D metal–organic frameworks (MOFs)7–16 have emerged as a new class of hybrid inorganic–organic 2D materials marked by their intrinsic structural porosity, chemical tunability and high electrical conductivity.17 These materials are synthesized through a bottom-up approach, by combining aromatic organic moieties that can ligate square-planar metal ions. The delocalization of π electrons in the 2D MOFs leads to high electrical conductivity that is uncommon in conventional MOFs. This enables a wide range of new functionalities, including chemiresistive sensing15,18 and electrochemical oxygen reduction.19 Finally, the delocalized electronic states expected for these materials may give rise to exotic new physics including topological insulator behavior,20,21 the quantum anomalous Hall effect,22 and flat-band ferromagnetism.23

The 2D MOFs synthesized so far8–15 are made from benzene- or triphenylene-derived ligands with imine, phenol, or thiophenol functionalities that bind a variety of late-transition-metal ions in a square planar environment. These 2D MOFs exhibit honeycomb lattices that stack in the third dimension, in a manner analogous to graphene and graphite. A schematic of such 2D honeycomb lattices, contrasted with 2D square lattices, is shown in Fig. 1. Although the presence of metals and the greater compositional variety of 2D MOFs allows significant modulation of band dispersion, absolute band energy, Fermi energy, and density of states, all 2D MOFs studied thus far have hexagonal symmetry and consequently exhibit electronic structures that mimic that of graphene symmetry-wise. Deviations from hexagonal symmetry should give rise to band structures that deviate entirely from those of hexagonal lattices, giving rise to electronic properties that diverge fundamentally from those expected from hexagonal 2D symmetry. It is thus highly desirable to explore 2D MOFs made from ligands and metal nodes that enforce non-hexagonal lattice geometries. In this article, we use first-principles calculations to computationally design metal-phthalocyanine (MPc)-based 2D MOFs with square lattices. Our study reveals rich magnetic behavior in this new class of 2D MOFs. In particular, we find that a charge-neutral MOF made from octaaminophthalocyanines metallated by Mn2+ ions and bound by Ni2+ ions through the phenylenediimine linkages should exhibit high-temperature ferromagnetic half-metallicity in monolayer form. In the bulk phase, MnPc based 2D MOFs exhibit strong intralayer but relatively weak interlayer magnetic coupling that is sensitive to the relative interlayer in-plane displacement. These results indicate the possibility of engineering the magnetic behavior of bulk or few-layer systems by controlling the interlayer stacking geometries.

image file: c6sc05080h-f1.tif
Fig. 1 Schematics of (a) honeycomb and (b) square lattices.

Results and discussion

Structural and magnetic ground states of monolayers

The MPc based 2D MOFs are constructed by conjugating 2,3,9,10,16,17,23,24-octaamino-substituted MPc molecules with Ni2+ ions, akin to the synthetic route of triphenylene-based 2D MOFs.11 In these materials, the amino groups are deprotonated and become imino groups. The nitrogen atoms of the imino groups in MPc molecules coordinate with Ni ions in square-planar geometry, resulting in long-range planar crystalline structures with square lattices, NiMPc, illustrated in Fig. 2. In this study, we considered materials containing metal-phthalocyanines (MPc) metallated with first row transition metals: M = Cr, Mn, Fe, Co, Ni, Cu, Zn. In addition to being thermally and chemically stable, these MPc molecules possess a wide range of intriguing electrical, optical and magnetic properties, which have been extensively investigated in the past.24–26 Despite their popularity in molecular materials, MPc moieties have not been frequently studied as components of covalent polymers.27–29 The rapid progress being made in synthesizing metal-ion mediated, π-conjugated 2D MOFs8–15 calls attention to MPc-based 2D MOFs, whose unique lattice symmetry, chemical tunability, and rich physical properties could distinguish them from other existing 2D MOFs.
image file: c6sc05080h-f2.tif
Fig. 2 Molecular models of octaamino-metal phthalocyanine (MPc) and NiMPc two-dimensional (2D) metal–organic frameworks (MOFs). M = Cr, Mn, Fe, Co, Ni, Cu and Zn.

Our first-principles calculations are based on density functional theory (DFT),30,31 as implemented in the Vienna ab initio simulation package (VASP).32,33 Exchange–correlation functional of the Perdew–Burke–Ernzerhof (PBE) form within the generalized gradient approximation (GGA)34,35 was used. Coulomb and exchange interactions of the localized d orbitals in TM elements were treated in the framework of the DFT+U method,36 using the Dudarev approach.37 The effective Coulomb (U) and exchange (J) parameters in DFT+U are U = 4 eV and J = 1 eV. This set of (U, J) parameters for TM d orbitals have been tested extensively in MPc-based systems.28,38–40 We have also tested a slightly different set of (U, J) parameters for MPc recently proposed by Brumboiu et al.41 and obtained results with quantitative similarity. Other details of the computational methods can be found in the ESI.

To determine the structural and magnetic ground states of the 2D MOFs, structural models of the 2D MOFs in unit cells were first fully relaxed in DFT+U with spin-polarized conditions. After relaxation all the systems adopt structures that can be identified as 2D square lattices in the p4m plane group. The CrPc, MnPc, FePc, CoPc and CuPc 2D MOFs exhibit unit cell magnetization of 4, 3, 2, 1, and 1 Bohr magneton (μB) respectively. The NiPc and ZnPc-based materials are non-magnetic because the TM atoms in these MPc have closed shells.26 In the systems that exhibit magnetization, the magnetic moments are predominantly localized on the TM atoms in the MPc moieties. The square-planar Ni ions connecting the MPc moieties do not carry magnetic moments that can be attributed to localized d electrons, since the Ni ions do not possess unpaired electrons in the square-planar coordination mode.

The magnetic ground states of the 2D MOFs with unit-cell magnetization are determined by comparing the energies of the systems with ferromagnetic (FM) or antiferromagnetic (AFM) coupling between the TM atoms in the MPc moieties. This is carried out in a image file: c6sc05080h-t1.tif structural supercell, which would be the magnetic unit cell if the system assumed AFM coupling, as illustrated in Fig. 3. Denoting the system energy per image file: c6sc05080h-t2.tif supercell by EAFM and EFM for AFM and FM coupling, respectively, the magnetic exchange energy per supercell is defined as EexEAFMEFM. Accordingly, the zero-temperature magnetic ground state of the system would be FM if Eex > 0, and AFM if Eex < 0. In Table 1 we list the computed values of Eex and the energy band gap corresponding to the magnetic ground states. Among the various systems investigated here, the CrPc, FePc, CoPc and CuPc-based systems are all semiconductors with weak AFM ground states. The calculated magnetic ground states of these 2D MOFs are in line with those obtained by Zhou et al. for MPc-based organometallic monolayers in the absence of the square-planar Ni linkages.28 Surprisingly, the MnPc-based system, NiMnPc, exhibits a FM metallic ground state with substantial exchange energy (Eex = 183 meV). In contrast to the rest of the NiMPc materials considered here, the exchange energy of NiMnPc, 183 meV, is significantly larger than the MnPc monolayer system studied by Zhou et al., 62 meV (40 meV using our calculation setups), when evaluated in the image file: c6sc05080h-t3.tif superlattice. In the following, we investigate the origin of the unusually strong FM coupling in the NiMnPc monolayer system.

image file: c6sc05080h-f3.tif
Fig. 3 Magnetic ground states of NiMPc 2D MOFs. Red and blue arrows denote spin-up and spin-down local magnetic moments, respectively. (a) and (b) Illustrate ferromagnetic (FM) and antiferromagnetic (AFM) coupling between metal ions in MPc moieties. The shaded regions in (a) and (b) indicate the structural unit cell and the image file: c6sc05080h-t4.tif supercell, each corresponding to the magnetic unit cell of the FM and the AFM ground states, respectively.
Table 1 Total magnetic moments per unit cell (M), exchange energies per image file: c6sc05080h-t5.tif supercell (Eex), and energy band gaps (Eg) of different MPc based 2D MOFs calculated using DFT+U
Cr Mn Fe Co Ni Cu Zn
M (μB) 4 3 2 1 0 1 0
E ex (meV) −36 183 −20 −19 −0.4
E g (eV) 0.35 Half-metal 0.28 0.35 0.35 0.30 0.30

The electronic band structure of NiMnPc near the Fermi level and the corresponding partial density of states (PDOS) are shown in Fig. 4 (the band structures of the other NiMPc 2D MOFs are shown in Fig. S1). Notably, the electronic states of NiMnPc in the vicinity of Fermi level are completely spin-polarized. The system is therefore a half-metal whose low-energy electronic phenomena would be dominated by one type of spin-polarized electrons, which could be relevant for spintronic applications. Analysis of the PDOS of bands near the Fermi level reveals that nearly all the bands are derived from the hybridization of π orbitals, namely the symmetry-equivalent dxz and dyz orbitals of Mn and Ni, as well as the pz orbitals of N and C. In combination with the metallic nature of the system, it can be concluded that the π orbitals of the system are delocalized. We note that because the electronic bands near the Fermi level are relatively dispersive, the MnPc based 2D MOFs should also have good electrical conductivity, a highly sought-after property for this class of materials.17

image file: c6sc05080h-f4.tif
Fig. 4 Electronic band structure and partial density of states (PDOS) of NiMnPc near the Fermi level. In the left panel, the spin-up bands are in red, and the spin-down bands are in blue. The states near the Fermi level are mainly derived from π orbitals, which include the dxz + dyz orbitals of Mn and Ni, and the pz orbitals of N and C. The PDOS of these orbitals are shown in the right panel.

Mechanism of magnetic exchange

To understand the nature of the unusually strong FM coupling in NiMnPc, we plotted the spin density isosurface of the system and the PDOS of Mn d orbitals in Fig. 5. As shown in Fig. 5a, positive spin density (ρρ > 0) mainly appears on the Mn atoms, and is associated with the localized d electrons. The adjacent nitrogen atoms that bond directly with Mn are oppositely polarized. Notably, the electron densities surrounding the square-planar Ni–N moieties and carbon atoms are also slightly polarized with an oscillatory character in real space. The spin polarization of these bridging units can be attributed to the polarization of the delocalized π orbitals as shown in Fig. 4. The PDOS of the Mn d orbitals in Fig. 5b shows that, except for a portion of the dxz/dyz orbitals that hybridize with the π orbitals of other atoms, the d orbitals of the Mn generally have very small dispersion, consistent with their localized nature. Given the large distance between the Mn atoms in neighboring MnPc moieties (∼18.3 Å), direct exchange between the d orbitals of Mn is rather weak, and thus cannot be responsible for the strong FM coupling between the Mn atoms. We conclude instead that the coupling occurs through the indirect exchange via the metallic π orbitals. Namely, the localized spin moments of the Mn atoms in individual MnPc moieties polarize the delocalized π orbital through direct exchange interaction. The π orbitals further mediate FM coupling of the Mn atoms at neighboring sites, as FM coupling between different MnPc moieties results in π orbitals with a consistent sign of polarization. The latter facilitates the hybridization of the π orbitals on different MnPc moieties and reduces the energy of the system. This indirect exchange through delocalized electrons is similar to the Zener d–p exchange42 or the RKKY43–45 exchange mechanism. One possible reason why this mechanism is operative in NiMnPc is that the dxz/dyz orbitals of Mn participate in the hybridization of π orbitals, leading to metallicity, stronger exchange interaction between the localized d orbitals and the π orbitals, and hence FM coupling. In contrast, in other NiMPc 2D MOFs, including CrPc, FePc, CoPc, NiPc, CuPc and ZnPc based materials, the filling of the d orbitals of TM atoms, and their energy alignment with respect to the π orbitals of the phthalocyanine ring and the linking Ni atoms, are not conducive to orbital hybridization. The localized magnetic moments of the TM atoms can therefore only interact through weak super-exchange or direct exchange, leading to AFM ground states with small exchange energies. The underlying cause of this difference in d electron hybridization can be clearly seen in the computed highest occupied molecular orbital (HOMO) of octaamino-MPc molecules, as shown in Fig. 6. Importantly, among the octaamino-MPc molecules, only the HOMO of MnPc has substantial contribution from the d orbitals of the TM atom and is fully delocalized over the entire molecule. This unique feature of strong hybridization between the metal d orbitals and the π electrons of ligand in octaamino-MnPc also occurs in amino-free MnPc molecules, as has been demonstrated both experimentally and theoretically.41,46–48
image file: c6sc05080h-f5.tif
Fig. 5 Spin density isosurface of NiMnPc and PDOS of Mn d orbitals. (a) Top and side view of the spin density isosurface at absolute spin-density value |ρρ| = 0.01 electrons per Å3. Red and blue isosurfaces correspond to positive and negative spin density, respectively. (b) PDOS of Mn dxy, dxz, dyz, dz2, dx2y2 orbitals. The PDOS of symmetry-equivalent dxz and dyz overlap with each other.

image file: c6sc05080h-f6.tif
Fig. 6 Visualization of the highest occupied molecular orbitals (HOMOs) of octaamino-MPc molecules by charge density isosurfaces. (a) Strong hybridization is evident between the dπ orbitals of Mn and the π orbitals of the ligand in octaamino-MnPc. (b) No contribution is observed from the metal d orbitals in the HOMO charge density isosurface of octaamino-MPc, where M is Cr, Fe, Co, Ni, Cu, or Zn. Both top and side views of the charge density isosurfaces (ρ = 0.001 electrons per Å3) are shown.

The square-planar Ni–N moieties in the MnPc-based 2D MOFs also play a key role in the strong FM exchange between the Mn atoms by effectively mediating the π electron conjugation in different octaamino-MPc molecules. This is consistent with the fact that 2D MOFs with some of the highest measured electrical conductivity were realized by π-conjugation through square-planar Ni–N linkage.9,11 We also find that replacing the nitrogen linker atoms in the square-planar Ni–N moieties by sulfur or oxygen does not significantly increase or decrease the magnetic exchange energy, suggesting that catechol and thiocatechol moieties may also mediate effective π electron conjugation (see also Fig. S2).

Monte-Carlo simulation of ferromagnetic phase transition

At elevated temperature, entropy favors disordered alignment of magnetic moments, causing NiMnPc to undergo a FM to paramagnetic phase transition at the Curie temperature, Tc. A critical issue affecting the theoretical determination of Tc is constructing an appropriate Hamiltonian that captures the nature of magnetic interactions in the system. The simplest model that describes isotropic interactions between local magnetic moments is the Heisenberg model image file: c6sc05080h-t6.tif, where [S with combining right harpoon above (vector)]i and [S with combining right harpoon above (vector)]j are the spin operators of magnetic moments at neighboring site i and j, and J is the exchange constant. However, it is well-known that in two dimensions, the isotropic Heisenberg model does not have spontaneous symmetry breaking at non-zero temperature due to the divergent contributions of gapless spin-waves.49 In reality, magnetic anisotropy due to spin–orbit coupling and/or magnetic dipolar interactions almost always exists, such that long-range magnetic order is usually observed with the addition of even a small amount of anisotropy.50–55 An extreme case is the 2D Ising model image file: c6sc05080h-t7.tif, which corresponds to infinitely large magnetic anisotropy, such that the spin vectors can only point in the z-direction (Si = +S or −S). Accordingly, the Ising model leads to FM or AFM ground states at non-zero temperature. Even though it has been used widely to calculate the Tc of 2D magnetic systems,28,56,57 it is not yet clear to what extent the Ising model overestimates Tc in most experimental situations, where there is a finite amount of magnetic anisotropy.

We attempt to predict more realistic Tc values for 2D FM systems by using two models that consider finite magnetic anisotropy. The first one (hereafter referred as Model 1) is the anisotropic Heisenberg model, which incorporates magnetic anisotropy by exchange anisotropy:

image file: c6sc05080h-t8.tif(1)
where [S with combining right harpoon above (vector)]i = (Six, Siy, Siz) denotes a unit vector of a classical magnetic moment at site i, the sum 〈ij〉 is over nearest neighbor pairs, J is the magnetic exchange constant, and γ is a dimensionless parameter that characterizes the strength of exchange anisotropy. For FM coupling between the spins, the values of J and (1 − γ) are positive. The magnetic “easy axis”, along which the exchange interaction is stronger, is chosen to be the z direction. The exchange anisotropy parameter γ can therefore take any values between 0 and 1. For γ = 0, the 2D isotropic Heisenberg model is obtained; γ = 1 gives an Ising-like model, that differs from the classical 2D Ising model in that the magnetic moments are still allowed to point in any arbitrary direction, not just + S or −S.

In Model 2, we assume isotropic exchange interactions, but take into consideration the presence of uniaxial single-ion anisotropy. This model enforces an energy penalty on magnetic moments pointing to directions other than the magnetic easy axis (chosen to be the z axis), irrespective of the spins at neighboring sites:

image file: c6sc05080h-t9.tif(2)

In this model, the strength of magnetic anisotropy, represented by γ, can take any positive values. A prefactor 2 is added before γ to facilitate comparison between Model 1 and Model 2, based on the energy cost of rotating the spins. The case γ = 0 in this model again corresponds to the 2D isotropic Heisenberg model, whereas the case γ = +∞ represents the classical 2D Ising model, since spin components in directions other than the magnetic easy axis are now fully penalized.

We use classical Monte Carlo simulations to determine Tc with Model 1 and Model 2 as a function of the magnetic anisotropy parameter γ. Tc is expressed as a function of J/kB, where kB is the Boltzmann constant. The Monte Carlo simulations were performed using the Metropolis algorithm, and Tc was determined by finding the intersection of Binder cumulants58 for two system sizes (40 × 40 and 80 × 80, consisting of 1600 and 6400 spins respectively). The calculated Tc as a function of γ for Model 1 and Model 2 are plotted in Fig. 7. Our simulations confirm previous studies52–55 that even a small amount of magnetic anisotropy can stabilize a FM ground state at finite temperature in the 2D Heisenberg model. When the values of γ are small (i.e. 0.001 < γ < 0.1), the Tc values determined by Model 1 and Model 2 were found to be close, as shown in the inset of Fig. 7. However, when γ > 0.1, the Tc predicted by Model 2 becomes noticeably larger. Because the value of γ in Model 1 is restricted in the range between 0 and 1, the maximum Tc of this model is reached at γ = 1.0, which gives Tc ≈ 0.9J/kB. In comparison, Model 2 predicts Tc ≈ 1.2J/kB at γ = 1.0. Importantly, in the range of γ from 0 to 1, the Tc predicted by both models are significantly lower than that of the 2D Ising model, which has an analytical solution, image file: c6sc05080h-t10.tif.59 For Model 2, the Ising limit is only asymptotically reached at large values of γ (γ > 100, see Fig. S3). Hence, caution must be exercised when using the classical Ising model to predict the Tc of 2D magnetic systems. Indeed, our results here show that in the absence of particularly strong magnetic anisotropy, the Ising model significantly overestimates Tc. Although quite dramatic, this overestimation is not surprising because in the Ising model the size of the accessible phase space for thermal excitation is significantly smaller compared to continuous spin models, low-energy excitations such as spin-waves being absent in the Ising model.

image file: c6sc05080h-f7.tif
Fig. 7 The Curie temperature (Tc) of 2D spin models in square lattices as a function of magnetic anisotropy parameter γ. The results for Model 1 (2D Heisenberg model with exchange anisotropy) and Model 2 (2D Heisenberg model with single-ion anisotropy) are plotted together. The inset shows Tc at small values of γ using log-scale horizontal axis. The Tc projected using the classical 2D Ising model is shown as a dashed line.

With a model that more comprehensively considers magnetic anisotropy in hand for magnetic phase transitions in the 2D space, we calculated the Tc of NiMnPc on the basis of the model parameters obtained from first-principles calculations. The exchange constant J in this system is related to the exchange energy Eex evaluated in the image file: c6sc05080h-t11.tif supercell by J = (1/8)Eex, which gives J = 22.9 meV. The magnetic anisotropy constant γ can be evaluated by relating it to the exchange constant J and the magnetic anisotropy energy (MAE), the latter of which can calculated from DFT+U with spin–orbit coupling.60 The value of MAE is determined as MAE ≡ EE, where E and E represent the energies of the system when the spin moment is parallel or perpendicular to the 2D plane of the MOF monolayer, respectively. A positive MAE implies that out-of-plane spin alignment is energetically more favorable than in-plane alignment. Our calculation gives MAE = 0.74 meV for a unit cell of NiMnPc monolayer. Although it is possible to specifically calculate the contribution of both exchange anisotropy and single-ion anisotropy to MAE from first principles, we have shown in Fig. 7 that the Tc values predicted by Models 1 and 2, corresponding to exchange anisotropy and single-ion anisotropy, respectively, are similar for small γ. In both models, the calculated MAE is related to γ through MAE = 2γJ. This gives γ = 0.017, which corresponds to Tc ≈ 0.66J/kB, or 170 K. In contrast, the frequently used 2D Ising model, where Tc ≈ 2.269J/kB, gives a transition temperature value higher than 600 K. Notably, when the square-planar NiN4 moiety in NiMnPc is changed to NiO4 or NiS4 (i.e., the phenylenediamine units are changed to catechols or thiocatechols) the predicted Tc are approximately 230 K and 150 K, respectively. Although the Tc values for the NiN4, NiO4, and NiS4 analogues of MnPc-based MOFs predicted using our realistic anisotropy model are still below room temperature, they sit above the liquid nitrogen temperature and thus present significant potential for practical applications (we note that we have ignored the phonon excitation and quantum fluctuation effects in calculating the Tc). Growth of these 2D MOFs as single layers on appropriate substrates, as might be expected under experimental conditions, could act as additional symmetry-breaking measures to enhance the magnetic anisotropy,61–63 further increasing Tc.

Bulk structural and magnetic properties

The discussion above pertains to single layers of 2D MOFs. The experimental realization of such single-layer structures notwithstanding, 2D MOFs synthesized thus far form bulk structures, where multiple monolayers stack to form layered bulk structures. To determine the most favorable interlayer stacking for NiMnPc, we considered two potential symmetries for bulk lattices: monoclinic and tetragonal (see Fig. S4). For the monoclinic lattice, the bulk unit cell contains one monolayer, with perpendicular a, b lattice vectors residing in the monolayer plane, and the out-of-plane c vector determining the interlayer separation and the relative (ab) displacement between neighboring layers. In the tetragonal lattice, we considered unit cells containing two monolayers with different relative (ab) displacements. In this case, the c vector is always normal to the (ab) plane and thus determines only the interlayer separation, not the relative layer displacement. Each of these unit cells gives rise to a large number of possible bulk stacked structures with different interlayer separations and relative in-plane displacements. We started from the relaxed structure of each monolayer and used DFT+U in combination with van der Waals correction in the form of DFT-D2 (ref. 64) to calculate the system energy for bulk models. In doing so, we considered both FM and AFM interlayer coupling for all structures and found that the lowest energy configuration adopts a tetragonal lattice with interlayer separation of ∼3.1 Å and relative in-plane displacement of ∼2.3 Å along both the a and b axes, as illustrated in Fig. 8a. The interlayer magnetic coupling was found to be weakly AFM, with exchange energy Eex ≈ −6 meV per unit cell. Notably, band structure calculations suggest that this lowest-energy bulk lattice is expected to be metallic (Fig. 8b); the experimental realization of this material would thus contribute a rare example of a highly electrically conductive MOF.17Fig. 8c shows the potential energy surface of NiMnPc as a function of the relative in-plane displacement between two neighboring layers in the tetragonal unit cell, with a fixed interlayer separation of 3.1 Å. The complex energy landscape results from a competition between the energy cost of orbital overlap between orbitals on neighboring layers (especially the Mn orbitals) and the energy gain from enhanced van der Waals interactions when neighboring layers have small relative in-plane displacement, as can be seen from Fig. S5. Expectedly, the interlayer magnetic coupling is itself sensitive to the relative in-plane displacement because the orbitals mediating magnetic coupling also vary with the displacement. This gives rise to oscillatory FM or AFM ground states depending on the relative displacement (Fig. 8d). Taken together, our results suggest that bulk NiMnPc should exhibit strong intralayer and relatively weak interlayer magnetic coupling that are sensitive to the stacking mode, portending the possibility that the magnetic behavior of bulk or few-layer systems can be rationally controlled by varying the interlayer stacking geometry.
image file: c6sc05080h-f8.tif
Fig. 8 Stacking mode and interlayer magnetic coupling in bulk NiMnPc. (a) Molecular structure of the lowest-energy stacking configuration as determined by DFT calculations. The lowest-energy configuration has interlayer distance ∼3.1 Å, and exhibits relative in-plane displacements of ∼2.3 Å along both a and b axes between adjacent layers. The dashed boxes in the top and side views outline the tetragonal unit cell of the bulk structure. (b) Calculated electronic band structure of the lowest-energy stacking configuration. (c) Potential energy surface generated by relative displacement of the two layers in the unit cell along a and b axes, while fixing the interlayer distance at the value of lowest-energy configuration. (d) Variation of the energy difference between interlayer AFM and FM coupling (EAFMEFM) while the two layers are displaced as in (c).


We demonstrate using first-principles calculations that metal phthalocyanine-based 2D MOFs conjugated by square-planar nickel ions exhibit rich structural and magnetic behavior. In particular, the MnPc-based materials are expected to exhibit half-metallicity and ferromagnetic ordering below 170 K, a remarkably large Curie temperature for a monolayer material. We demonstrate that such strong magnetic coupling results from the strong hybridization between the π-symmetry orbitals of Mn and the phthalocyanine ring, with the square-planar NiN4 moieties playing a key role in facilitating efficient magnetic exchange and π electron delocalization. We have also devised new theoretical means to include magnetic anisotropy in evaluating Tc in 2D magnetic systems. In particular, we show that in the absence of strong magnetic anisotropy, the widely used 2D Ising model can lead to overestimation of Tc by several times. These results suggest that π electron conjugation between magnetic ions and ligand groups, together with enhanced magnetic anisotropy through the incorporation of heavy elements with sizable spin–orbit coupling, are effective strategies to design ferromagnetic 2D MOFs with high Tc. The materials studied in this work are based on existing molecular building blocks and we are thus confident that they could be realized experimentally. Indeed, synthetic methods for octaaminophthalocyanine exist,65 and the square-planar conjugation of phenylenediamine linkers with Ni is also well documented and has already been used to make highly conductive 2D MOFs.11 Considering the tremendous amount of chemical tunability in the MPc-based 2D MOFs stemming from metal replacement24,25 and the utilization of different linking groups and ligating D4h metal ions,17 the new class of square lattices proposed here represents a rich platform for new solid-state physics phenomena in MOFs, with the promise of a wide range of potential applications in spintronics and electronics.


This work was supported by the Center for Excitonics, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Basic Energy Sciences under Award No. DE-SC0001088. M.D. thanks the Sloan Foundation, the Dreyfus Foundation, and the Research Corporation for Science Advancement for financial support. We thank Dr Dennis Sheberla, Dr Michael G. Campbell, Dr Efren Navarro-Moratalla, and Melis Tekant for helpful discussions. Figures involving atomic structures and charge density isosurfaces were rendered using VESTA.66


  1. K. S. Novoselov, D. Jiang, F. Schedin, T. J. Booth, V. V. Khotkevich, S. V. Morozov and A. K. Geim, Proc. Natl. Acad. Sci. U. S. A., 2005, 102, 10451–10453 CrossRef CAS PubMed.
  2. A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov and A. K. Geim, Rev. Mod. Phys., 2009, 81, 109–162 CrossRef CAS.
  3. Q. H. Wang, K. Kalantar-Zadeh, A. Kis, J. N. Coleman and M. S. Strano, Nat. Nanotechnol., 2012, 7, 699–712 CrossRef CAS PubMed.
  4. M. Chhowalla, H. S. Shin, G. Eda, L.-J. Li, K. P. Loh and H. Zhang, Nat. Chem., 2013, 5, 263–275 CrossRef PubMed.
  5. X. Xu, W. Yao, D. Xiao and T. F. Heinz, Nat. Phys., 2014, 10, 343–350 CrossRef CAS.
  6. F. Xia, H. Wang, D. Xiao, M. Dubey and A. Ramasubramaniam, Nat. Photonics, 2014, 8, 899–907 CrossRef CAS.
  7. M. Hmadeh, Z. Lu, Z. Liu, F. Gandara, H. Furukawa, S. Wan, V. Augustyn, R. Chang, L. Liao, F. Zhou, E. Perre, V. Ozolins, K. Suenaga, X. F. Duan, B. Dunn, Y. Yamamto, O. Terasaki and O. M. Yaghi, Chem. Mater., 2012, 24, 3511–3513 CrossRef CAS.
  8. T. Kambe, R. Sakamoto, K. Hoshiko, K. Takada, M. Miyachi, J.-H. Ryu, S. Sasaki, J. Kim, K. Nakazato, M. Takata and H. Nishihara, J. Am. Chem. Soc., 2013, 135, 2462–2465 CrossRef CAS PubMed.
  9. T. Kambe, R. Sakamoto, T. Kusamoto, T. Pal, N. Fukui, K. Hoshiko, T. Shimojima, Z. Wang, T. Hirahara, K. Ishizaka, S. Hasegawa, F. Liu and H. Nishihara, J. Am. Chem. Soc., 2014, 136, 14357–14360 CrossRef CAS PubMed.
  10. J. Cui and Z. Xu, Chem. Commun., 2014, 50, 3986–3988 RSC.
  11. D. Sheberla, L. Sun, M. A. Blood-Forsythe, S. Er, C. R. Wade, C. K. Brozek, A. Aspuru-Guzik and M. Dincă, J. Am. Chem. Soc., 2014, 136, 8859–8862 CrossRef CAS PubMed.
  12. T. Pal, T. Kambe, T. Kusamoto, M. L. Foo, R. Matsuoka, R. Sakamoto and H. Nishihara, ChemPlusChem, 2015, 80, 1255–1258 CrossRef CAS.
  13. A. J. Clough, J. W. Yoo, M. H. Mecklenburg and S. C. Marinescu, J. Am. Chem. Soc., 2015, 137, 118–121 CrossRef CAS PubMed.
  14. R. Dong, M. Pfeffermann, H. Liang, Z. Zheng, X. Zhu, J. Zhang and X. Feng, Angew. Chem., Int. Ed., 2015, 54, 12058–12063 CrossRef CAS PubMed.
  15. M. G. Campbell, D. Sheberla, S. F. Liu, T. M. Swager and M. Dincă, Angew. Chem., Int. Ed., 2015, 54, 4349–4352 CrossRef CAS PubMed.
  16. H. Maeda, R. Sakamoto and H. Nishihara, Langmuir, 2016, 32, 2527–2538 CrossRef CAS PubMed.
  17. L. Sun, M. G. Campbell and M. Dincă, Angew. Chem., Int. Ed., 2016, 55, 3566–3579 CrossRef CAS PubMed.
  18. M. G. Campbell, S. F. Liu, T. M. Swager and M. Dincă, J. Am. Chem. Soc., 2015, 137, 13780–13783 CrossRef CAS PubMed.
  19. E. M. Miner, T. Fukushima, D. Sheberla, L. Sun, Y. Surendranath and M. Dincă, Nat. Commun., 2016, 7, 10942 CrossRef CAS PubMed.
  20. Z. F. Wang, N. Su and F. Liu, Nano Lett., 2013, 13, 2842–2845 CrossRef CAS PubMed.
  21. B. Zhao, J. Zhang, W. Feng, Y. Yao and Z. Yang, Phys. Rev. B: Condens. Matter Mater. Phys., 2014, 90, 201403 CrossRef.
  22. L. Dong, Y. Kim, D. Er, A. M. Rappe and V. B. Shenoy, Phys. Rev. Lett., 2016, 116, 096601 CrossRef PubMed.
  23. M. G. Yamada, T. Soejima, N. Tsuji, D. Hirai, M. Dincă and H. Aoki, Phys. Rev. B, 2016, 94, 081102 CrossRef.
  24. N. B. McKeown, Phthalocyanine Materials: Synthesis, Structure and Function, Cambridge University Press, Cambridge, 1998 Search PubMed.
  25. K. M. Kadish, K. M. Smith and R. Guilard, The Porphyrin Handbook, Academic Press, San Diego, 2003,vol. 15–20 Search PubMed.
  26. J. Bartolome, C. Monton and I. K. Schuller, Molecular Magnets: Physics and Applications, 2014, pp. 221–245 Search PubMed.
  27. M. Abel, S. Clair, O. Ourdjini, M. Mossoyan and L. Porte, J. Am. Chem. Soc., 2010, 133, 1203–1205 CrossRef PubMed.
  28. J. Zhou and Q. Sun, J. Am. Chem. Soc., 2011, 133, 15113–15119 CrossRef CAS PubMed.
  29. E. L. Spitler and W. R. Dichtel, Nat. Chem., 2010, 2, 672–677 CrossRef CAS PubMed.
  30. P. Hohenberg and W. Kohn, Phys. Rev., 1964, 136, B864–B871 CrossRef.
  31. W. Kohn and L. J. Sham, Phys. Rev., 1965, 140, A1133–A1138 CrossRef.
  32. G. Kresse and J. Furthmüller, Comput. Mater. Sci., 1996, 6, 15–50 CrossRef CAS.
  33. G. Kresse and J. Furthmüller, Phys. Rev. B: Condens. Matter Mater. Phys., 1996, 54, 11169–11186 CrossRef CAS.
  34. J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Singh and C. Fiolhais, Phys. Rev. B: Condens. Matter Mater. Phys., 1992, 46, 6671–6687 CrossRef CAS.
  35. J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett., 1996, 77, 3865–3868 CrossRef CAS PubMed.
  36. V. I. Anisimov, F. Aryasetiawan and A. I. Lichtenstein, J. Phys.: Condens. Matter, 1997, 9, 767–808 CrossRef CAS.
  37. S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys and A. P. Sutton, Phys. Rev. B: Condens. Matter Mater. Phys., 1998, 57, 1505–1509 CrossRef CAS.
  38. A. Calzolari, A. Ferretti and M. B. Nardelli, Nanotechnology, 2007, 18, 424013 CrossRef PubMed.
  39. M. Bernien, J. Miguel, C. Weis, M. E. Ali, J. Kurde, B. Krumme, P. M. Panchmatia, B. Sanyal, M. Piantek, P. Srivastava, K. Baberschke, P. M. Oppeneer, O. Eriksson, W. Kuch and H. Wende, Phys. Rev. Lett., 2009, 102, 047202 CrossRef CAS PubMed.
  40. S. Lach, A. Altenhof, K. Tarafder, F. Schmitt, M. Ali, M. Vogel, J. Sauther, P. M. Oppeneer and C. Ziegler, Adv. Funct. Mater., 2012, 22, 989–997 CrossRef CAS.
  41. I. E. Brumboiu, S. Haldar, J. Lüder, O. Eriksson, H. C. Herper, B. Brena and B. Sanyal, J. Chem. Theory Comput., 2016, 12, 1772–1785 CrossRef CAS PubMed.
  42. C. Zener, Phys. Rev., 1951, 81, 440–444 CrossRef CAS.
  43. M. A. Ruderman and C. Kittel, Phys. Rev., 1954, 96, 99–102 CrossRef CAS.
  44. T. Kasuya, Prog. Theor. Phys., 1956, 16, 45–57 CrossRef.
  45. K. Yosida, Phys. Rev., 1957, 106, 893–898 CrossRef.
  46. M. Grobosch, C. Schmidt, R. Kraus and M. Knupfer, Org. Electron., 2010, 11, 1483–1488 CrossRef CAS.
  47. M. Grobosch, B. Mahns, C. Loose, R. Friedrich, C. Schmidt, J. Kortus and M. Knupfer, Chem. Phys. Lett., 2011, 505, 122–125 CrossRef CAS.
  48. J. Kugel, M. Karolak, J. Senkpiel, P. J. Hsu, G. Sangiovanni and M. Bode, Nano Lett., 2014, 14, 3895–3902 CrossRef PubMed.
  49. N. D. Mermin and H. Wagner, Phys. Rev. Lett., 1966, 17, 1133–1136 CrossRef CAS.
  50. C. Herring and C. Kittel, Phys. Rev., 1951, 81, 869–880 CrossRef.
  51. J. Frohlich and E. H. Lieb, Phys. Rev. Lett., 1977, 38, 440–442 CrossRef.
  52. K. Binder and D. P. Landau, Phys. Rev. B: Solid State, 1976, 13, 1140–1155 CrossRef CAS.
  53. R. P. Erickson and D. L. Mills, Phys. Rev. B: Condens. Matter Mater. Phys., 1991, 43, 11527–11530 CrossRef.
  54. P. A. Serena, N. Garcia and A. Levanyuk, Phys. Rev. B: Condens. Matter Mater. Phys., 1993, 47, 5027–5036 CrossRef.
  55. D. J. Priour, E. H. Hwang and S. Das Sarma, Phys. Rev. Lett., 2005, 95, 037201 CrossRef PubMed.
  56. H. L. Zhuang, Y. Xie, P. R. C. Kent and P. Ganesh, Phys. Rev. B: Condens. Matter Mater. Phys., 2015, 92, 035407 CrossRef.
  57. F. Wu, C. X. Huang, H. P. Wu, C. Lee, K. M. Deng, E. J. Kan and P. Jena, Nano Lett., 2015, 15, 8277–8281 CrossRef CAS PubMed.
  58. K. Binder, Z. Phys. B: Condens. Matter, 1981, 43, 119–140 CrossRef.
  59. L. Onsager, Phys. Rev., 1944, 65, 117–149 CrossRef CAS.
  60. P. Blonski and J. Hafner, Phys. Rev. B: Condens. Matter Mater. Phys., 2009, 79, 224418 CrossRef.
  61. S. Stepanow, A. Mugarza, G. Ceballos, P. Moras, J. C. Cezar, C. Carbone and P. Gambardella, Phys. Rev. B: Condens. Matter Mater. Phys., 2010, 82, 014405 CrossRef.
  62. I. G. Rau, S. Baumann, S. Rusponi, F. Donati, S. Stepanow, L. Gragnaniello, J. Dreiser, C. Piamonteze, F. Nolting, S. Gangopadhyay, O. R. Albertini, R. M. Macfarlane, C. P. Lutz, B. A. Jones, P. Gambardella, A. J. Heinrich and H. Brune, Science, 2014, 344, 988–992 CrossRef CAS PubMed.
  63. M. Gruber, F. Ibrahim, S. Boukari, H. Isshiki, L. Joly, M. Peter, M. Studniarek, V. Da Costa, H. Jabbar, V. Davesne, U. Halisdemir, J. J. Chen, J. Arabski, E. Otero, F. Choueikani, K. Chen, P. Ohresser, W. Wulfhekel, F. Scheurer, W. Weber, M. Alouani, E. Beaurepaire and M. Bowen, Nat. Mater., 2015, 14, 981–984 CrossRef CAS PubMed.
  64. S. Grimme, J. Comput. Chem., 2006, 27, 1787–1799 CrossRef CAS PubMed.
  65. F. Yuksel, S. Tuncel and V. Ahsen, J. Porphyrins Phthalocyanines, 2008, 12, 123–130 CrossRef CAS.
  66. K. Momma and F. Izumi, J. Appl. Crystallogr., 2008, 41, 653–658 CrossRef CAS.


Electronic supplementary information (ESI) available: (1) Computational methods; (2) electronic band structures of NiMPc MOF monolayers (Fig. S1); (3) calculation results for NiMnPc MOF monolayers with NiO4 and NiS4 moieties (Fig. S2); (4) ferromagnetic transition temperatures of 2D Heisenberg model with single-ion anisotropy (Fig. S3); (5) structure and energetics of bulk NiMnPc (Fig. S4 and S5); (6) atomistic coordinate data of NiMPc MOF monolayers. See DOI: 10.1039/c6sc05080h

This journal is © The Royal Society of Chemistry 2017