Electronic and magnetic properties of DUT-8 ( Ni )

Please note that technical editing may introduce minor changes to the text and/or graphics, which may alter content. The journal’s standard Terms & Conditions and the Ethical guidelines still apply. In no event shall the Royal Society of Chemistry be held responsible for any errors or omissions in this Accepted Manuscript or any consequences arising from the use of any information it contains. Accepted Manuscript


Introduction
4][5] Such compounds consist of certain metal centers (metal) within their so-called secondary building unit which are connected via specific organic linkers (organic) to build a crystalline network (framework).Inside MOFs, there can be pores of different sizes.Those pores enable the MOF to be utilized for absorption due to their size and resulting large surface area.Several subclasses of frameworks, including so-called paddle wheels, can nowadays be synthesized leading to a variety of different types of MOFs.Some MOFs are flexible, 6 meaning that several stable crystalline structures can be obtained by absorption and desorption.From the diversity of available metal centers, open shell transition metals can introduce a magnetic ground state.In general, the main research on MOFs is concerned with their absorption behaviour and how to increase their surface areas and porosities.On the other hand, the magnetism in MOFs has not been studied yet, but is thought to extend the already wide field of applications for MOFs.Furthermore there is a connection to the widely studied field of binuclear complexes, [7][8][9][10][11] which are the carriers of the magnetic interaction within DUT-8(Ni) and other MOFs.A special case are Ni(II) dimer systems, which have been studied in the last few years. 7,8Within Ni(II) dimer complexes, each Ni can be viewed as a d 8 -Ni(II).As a consequence the highest possible spin is S HS = 2 while the lowest spin is S LS = 0.In all discussions low-spin accounts for an antiparallel alignment of the spins at the magnetic centers while high-spin stands for a parallel alignment.However, the magnetic ground state in DUT-8(Ni) 12,13 (Fig. 1), with one Ni dimer per magnetic unit, is still unclear.
DUT-8(Ni) is flexible corresponding to two crystalline structures with so-called paddle wheels (DUT-8(Ni) open /DUT-8(Ni) closed ), with a ratio of the unit cell volumes of V DUT-8(Ni) open E 2.6ÁV DUT-8(Ni) closed .The Ni atoms in these structures are coordinated by four oxygens of the organic linkers (2,6-naphthalenedicarboxylate C 12 O 4 H 6 ), 12 which connect the magnetic units in two dimensions to form the paddle wheels.The paddle wheels are connected via dabco units (1,4-diazabicyclo[2Á2Á2]octane C 6 N 2 H 12 ). 12This leads to an additional coordination of the Ni by a nitrogen atom in the three-dimensional periodic structure of DUT-8(Ni).In connection to the structural flexibility a very interesting question arises concerning the magnetic properties of such structures.
The main question to answer in the first place is whether the structural change, which occurs due to absorption into the DUT-8(Ni) closed to gain DUT-8(Ni) open , alters the magnetic ground state.Such an alteration could be used as an additional measure of absorption and might even be sensitive to different absorbents, as it might lead to different strengths of the coupling.To gain insight into the magnetic behaviour of this compound, DFT calculations were performed.The magnetic properties are discussed on the basis of the calculated energy difference between a low-spin and a high-spin state DE LS-HS , which describes the favoured magnetic ordering qualitatively as well as quantitatively.
For further evaluation of the strength of the magnetic interaction the Heisenberg-Dirac-Van Vleck model is applied.This model uses a Hamiltonian which describes the magnetic interaction between spins with the so-called coupling constant J.The energy values as well as the resulting magnetizations from spin-polarized DFT calculations allow the computation of this coupling constant.It scales the energy differences to the respective total spins of the high-spin state and the low-spin state (see eqn (2)).The sign of J is given by the favoured magnetic ordering (positive for high-spin and negative for low-spin) while its absolute value is a measure for the strength of such ordering.

DFT calculations
All calculations have been carried out using the DFT 14,15 based program package QUANTUM ESPRESSO (QE). 16It uses a pseudopotential approach with a plane wave 17 basis set.
The projector-augmented wave (PAW) 18 method with a GGA-PBE 19 exchange-correlation functional was used ([element].pbe-(n)-kjpaw_psl.0.1.UPF, for details see http:\\www.quantum-espresso.org/pseudopotentials/).All calculations were performed spin-polarized.The van der Waals (vdW) interaction by Grimme 20 has been used to validate the results of the crystalline structures.The results of the QE calculations have been compared with calculations using the GPAW 21 program package.This was done to compare the results obtained with different basis sets (plane wave vs. local basis).GPAW is a real-space-grid DFT code based on the projector-augmented wave (PAW) method.Pseudo wavefunctions can be described in different ways (https://wiki.fysik.dtu.dk/gpaw/algorithms.html#algorithms), of which the LCAO ansatz 22 was used.All LCAO calculations are spin-polarized and were performed using a double-z basis set with one polarization function.The set of k-points and the exchange-correlation functional are the same like in the QE calculations.
The application of the GGA-PBE approach is justified by its robustness (see below), its computational efficiency and the minimum number of required parameters.These are advantages compared to hybrid functionals (e.g.B3LYP) or DFT + U, although such approaches can correct the problem of too much delocalization of d-electrons in certain ways.Other methods, e.g.CI or CASSCF, are prohibited concerning the computational effort for large systems like MOFs.Within the research of binuclear complexes, [9][10][11] different approaches to determine the magnetic properties of small molecules have been carried out.Such approaches include wave function based methods (e.g.CI 9 or CASSCF 10 ) and recently also DFT. 11 The B3LYP hybrid functional with different amounts of exact exchange has been compared to DDCI 23 (difference dedicated configuration interaction) and to a PBE functional. 11Calculations of the coupling constant show similar results for all methodologies.Especially the results obtained with B3LYP and PBE agree even quantitatively.Considering these results and taking into account the size of DUT-8(Ni), it is reasonable to apply the GGA-PBE as a good compromise between the computational effort and accuracy.Thus all calculations in the present work were performed within the GGA-PBE approximation.
The Brillouin-zone was sampled with a 3 Â 3 Â 3 Monkhorst-Pack grid (14 k-points per spin orientation).The kinetic energy cutoff is 90 Ry for the crystalline systems and 150 Ry for the model systems.Full geometry optimizations with and without the inclusion of the vdW interaction for the crystalline structures were performed.The optimizations were carried out on a 2 Â 2 Â 2 Monkhorst-Pack grid where convergence was achieved once all forces are smaller than 10 À3 Ry/a 0 E 0.025 eV Å À1 (a 0 is the Bohr radius).
It has to be considered that within the DUT-8(Ni) open unit cell the interaction between the Ni dimers is negligible (see Section 3.2).Thus the energy difference between a low-spin and a high-spin coupling DE LS-HS has to be divided by 2. This scales DE LS-HS, DUT-8(Ni) open to one Ni dimer.With that it can be compared to DE LS-HS,DUT-8(Ni) closed .
For the calculation of DUT-8(Ni) open , no absorbents have been taken into account.Thus a comparison of the structures as shown in Fig. 1 was performed.However no major changes in the magnetic behaviour are expected once certain substances are absorbed (as long as those substances do not carry some spin polarization on their own).

Heisenberg-Dirac-Van Vleck model
A description of the magnetic behaviour, as introduced by the coupling of local spins at different magnetic centers, is given by the Heisenberg-Dirac-Van Vleck (HDVV) Hamiltonian [24][25][26] ĤHDVV ¼ À2 where J ij is the so-called coupling constant between neighbouring spins while -S i and -S j are spin operators of the magnetic centers i and j.A high-spin state, meaning a parallel alignment of the spins, is indicated by a positive value of the coupling constant J ij while a negative value refers to a low-spin state and an antiparallel alignment.The coupling constant for dimers (i = 1 and j = 2) can This journal is © the Owner Societies 2015 be related to the total energies and the total spins of two magnetic orderings Such energies and total spins can be taken from spin-polarized DFT calculations.The local spins -S i and -S j at each centers do not have to be taken into account.The denominator in eqn (2) becomes six for a Ni dimer (S HS = 2, S LS = 0), two for a Cu dimer (S HS = 1, S LS = 0) and three for a mixed Ni-Cu dimer (S HS = 3/2, S LS = 1/2).This consideration leads to the values given in Table 3.For further analysis of the magnetic interaction, calculations on supercells (Fig. 2) were performed.The supercells allow the study of different magnetic coupling configurations (Fig. 3).Thus the magnetic interaction beyond a single unit cell is analyzed, as an evaluation of such couplings is not possible with single unit cells.The energetically most favoured interaction will be either within the magnetic dimer units (local interaction) or between them (global interaction).It can be shown, e.g. for the MOF HKUST-1, 27 that the magnetic interaction can be described by molecular model systems.For an in-depth insight into the magnetic interaction within DUT-8(Ni), several molecular model systems were generated (Fig. 4 and Table 1).This allows to analyze possible alterations of the magnetic ground state which are related to changes in the chemical environment of the magnetic centers rather than the structural change in the periodic systems.Thus several ways to influence the magnetic properties have been employed.For the initial models M1 (Ni 2 (HCOO) 4 (NH 3 ) 2 , Fig. 4) the coordinates of the atoms inside the magnetic unit within the crystalline systems were used.To achieve charge neutrality the carbon of the organic linkers as well as the nitrogen of the dabco units were terminated with hydrogen.These initial structures were modified (Table 1) to study the influence of the local environment of the Ni atoms on the magnetic properties.The modifications on the chemical environment include the addition of phenyl groups, the removal of NH 3 groups and the replacement of NH 3 groups with some other functional groups.Additionally, the magnetic centers themselves have been exchanged (Cu instead of Ni) without altering the chemical environment.The influence on the magnetic behaviour caused by such alterations was investigated.The initial models (M1) contain 26 atoms while all other models (Table 1) incorporate between 18 and 26 atoms.For the molecular model systems the unit cells have been chosen to impose a vacuum of 10 Å in each direction (x,y,z) to avoid interactions between neighbouring molecules.Therefore a low-spin state is expected to be favoured.The resulting coupling constant J in such a case is negative.Further the contribution to the states around the Fermi level E F does not exclusively come from the Ni d-states.The biggest contribution next to such states are the oxygen p-states (not shown in Fig. 5).Such contributions have been found in the electronic structure of e.g.MOF-5 28 as well.

Crystalline systems
The contribution of the different spin densities for DUT-8(Ni) closed is not symmetric anymore (Fig. 5b).However, one of the Ni contributes mainly to one of the spin species while the other magnetic center contributes to the other one (Ni 1 k, Ni 2 m).The spins at the two magnetic centers are aligned antiparallel and the overall magnetic ordering should be the low-spin state as well.
The calculated energy differences between a low-spin and a high-spin coupled system for the two crystalline system give a qualitative as well as quantitative insight into the magnetic ordering (the results for DUT-   In order to check how the vdW interaction influences the energy differences, the calculations were performed again with an included vdW correction using QE.The achieved results are qualitatively as well as quantitatively the same as the ones without such a correction.For further analysis, full geometry optimization for both structures with and without vdW interaction for both high-spin and low-spin configurations were performed.The resulting energy differences between such magnetic states change only slightly Thus for no optimization there is no change in the magnetic interaction while for fully geometry optimized structures small changes in the magnetic strength can be observed.Overall the vdW interaction plays no role in the determination of the magnetic ground state.

Supercell calculations
As magnetic interactions extend beyond elementary cells, supercells were generated.This allows the calculation of different magnetic couplings between several elementary cells.Further a comparison of global magnetic interactions and local ones can be carried out.A variety of different interactions which cannot be realized within elementary cells can be studied as well.As mentioned earlier, the organic linkers within the DUT-8(Ni) closed supercell are only implemented by the periodic boundary conditions.To verify this statement the interaction along the organic linkers in DUT-8(Ni) open has been analyzed.Calculating the energy difference between a low-spin and a respective high-spin ordering along the linkers gives DE LS-HS E 0.04 meV per unit cell.Thus the interaction along the organic linkers is small and does not have to be taken explicitly into account in the DUT-8(Ni) closed supercell.In comparison to the DUT-8(Ni) open supercell, this allows the study of the interaction of one magnetic dimer unit coupled antiparallelly to all others (second last entry in Table 2).Several magnetic orderings were introduced to the Ni dimers within the supercells and the resulting energies were compared to figure out which magnetic ordering is the energetically most favoured one (Table 2).
First the results indicate that the local low-spin ordering is the energetically most favoured coupling of all possible interactions.This confirms the results achieved by the calculations on single unit cells.Second they provide the information that for local high-spin behaviour, every global low-spin ordering shows lower energies than any global high-spin ordering.The only exception occurs for a high-spin coupling along the dabco linker units (the second coupling for DUT-8(Ni) closed in Table 2).There, a very small energy difference to the global high-spin case (mmmm) is found.An explanation is the interaction between Ni and N, which has a major influence on the magnetic behaviour (see Section 3.3).For the other global low-spin couplings similar energies are found where the coupling along the dabco units is always low-spin.Thus the interactions along the dabco units determine which coupling is favoured.A local high-spin state can be gained by altering the chemical environment of the magnetic centers, as shown in the next section.

Molecular models
To further study the magnetic interaction between the Ni atoms in DUT-8(Ni), several molecular model systems were generated.With such models, the chemical environment of the magnetic centers can be modified.Additionally, the magnetic centers themselves were exchanged (Table 1).The used molecular models (M1) (Fig. 4) show a similar behaviour considering the DE LS-HS in comparison to their respective crystal structures (Table 3).These systems were modified to study influences of different chemical environments as well as an exchange of magnetic centers.If such alterations would be implemented into the crystalline structures, a similar magnetic behaviour is expected.The inclusion of the phenyl groups (M2), which describe the system with organic linkers, shows basically no change in the magnetic ground state and can be neglected for all other models.
The removal of the NH 3 groups (M3, M4) is interesting in that it leads to a high-spin state due to the missing interaction between the N to one of the Ni.An increase in the coupling strength can be observed when removing both NH 3 groups.Without the N-Ni interaction, the electron density of the Ni 'relaxes' towards to the other Ni, which changes the interaction between them.The resulting ground state has a S HS = 1 total spin.An overlap of the electron densities of the Ni due to the mentioned relaxation is assumed.Thus there are two electrons in between the Ni with a low-spin coupling (mk) while the outer electrons couple high-spin (mm).This is confirmed by calculations of the electron localization functions (ELFs) revealing paired electrons in between the Ni.As there are no paired electrons in the M1 models, the previous assumption is validated.Transferring the alterations to the crystalline systems can be thought of as a removal of a dabco unit which separates the paddle wheels within the DUT-8(Ni). 12When synthesizing this MOF, a reduced amount of dabco units should lead to a high-spin state while for a higher amount of such linker groups, a low-spin state should be present.
A different magnetic behaviour is expected when exchanging Ni with Cu (M6, M7, M8), considering the differences of their electronic configuration ([Ar]4s 2 3d 8 to [Ar]3d 10 4s 1 ).For a Ni dimer, there are two unpaired d-electrons per Ni while there is only one unpaired s-electron each in the case of Cu.Thus in Ni two 3d-states with one electron each can couple while on the other hand in Cu the magnetic interaction occurs between only one electron per atom.The high-spin energies within the models correspond to a S HS,Ni = 2 and a S HS,Cu = 1 ground state, respectively.As known from ref. 29, the magnetic ground state in Cu complexes usually results in a low-spin state (S = 0) with a rather small DE LS-HS (at around 100 meV).This is confirmed by the calculations of the M6, M7 and M8 models.The results show that the favoured ground state for both Ni and Cu is the low-spin state.The energy difference between the low-spin state and the high-spin state DE LS-HS for Ni is larger than for Cu.An explanation is the larger number of electrons which contribute to the magnetic interaction for Ni.On the other hand the DE LS-HS in the Cu dimer has been found to be higher than in the mixed Ni-Cu dimer models.This can be explained by the magnetic configurations which have to be considered in such systems.Ni provides two unpaired electrons while Cu provides only one unpaired electron.A spin of S HS = 3/2 is obtained by the highspin state while the low-spin state shows a S LS = 1/2.Unlike in the Cu dimer system, there are always two parallelly aligned electrons from Ni which interact with either another parallelly aligned (mmm) or an antiparallelly aligned electron (mmk) from Cu.The separation between the high-spin state and the low-spin state is not as distinct as in the Cu dimer, where the electrons either align parallelly (mm) or antiparallelly (mk).All those alterations can be related to a metal exchange in the crystalline systems while leaving the rest of the framework unchanged.
Further alterations were implemented to the M1 models by exchanging both NH 3 with some other functional groups, such as H 2 O, CO 2 and HCN (M9, M10, M11).In analogy to changing ligands in Cu dimers 30 the expectations arise that major changes in the magnetic ground state occur and that the Ni-N interaction plays an important role for the magnetic interaction.For M9, M10 and M11 these expectations are met.In those models, the high-spin state is favoured and obtains a S HS = 1.For the interpretation again an ELF analysis can be used.In each of them, there is a clear indication of a bond between the two Ni, thus two spins are always mk.In M9 the Ni is coordinated by yet another oxygen instead of a N.While DE LS-HS is small in comparison to the other two systems, the magnetic ground state already changes to a high-spin state.This indicates that the Ni-N interaction is vital for the magnetism in the systems.For M10 the energy difference increases by about a factor of 5 compared to M9 due to the stronger interaction of the oxygen with its attached C. In M9, the hydrogens obtain single bonds to the oxygen.In M10 there is a double bond to C, thus redistributing the valence electrons of oxygen.In M11 a further increase in DE LS-HS is observed, even though the N is neither removed nor exchanged.However, in the M1 models the N forms three bonds to H, thus obtaining a sp 3 -hybridisation.For M11 one gets NRC-H corresponding to a sp-hybridisation at the carbon atom which leads to a higher separation between the low-spin state and the high-spin state.
All discussions so far were carried out for the DUT-8(Ni) open models.In the corresponding DUT-8(Ni) closed models a similar magnetic behaviour is found (Table 3).However, the quantitative values for DE LS-HS are smaller.This can be attributed to the deviation in the structural conformation between DUT-8(Ni) open and DUT-8(Ni) closed .The biggest differences are observed for the models M4 and M6, where the DUT-8(Ni) closed models show a very small energy difference between the high-spin state  .The corresponding distance in the crystals with a dabco unit is d crystal E 6.8 Å. Trends of the energy depending on the distance were investigated while two different magnetizations were imposed to the systems.As shown earlier, the magnetic ground state for M3 has been found to be a high-spin state.Thus the two M3 molecules were generated with a local high-spin state, while the global coupling of the system was chosen to be either a high-spin (mmmm) or a low-spin one (mmkk).First, the global interaction has basically no influence on the resulting total energy.Thus, the local interaction is of major importance and only one graph is displayed in Fig. 6.Second, a shallow energy minimum is found at a distance of about 2.9 Å.Whether this minimum can be realized in the crystalline structures is however questionable.Another observation from Fig. 6 is that for a Ni-Ni distance lower than 2.5 Å the system becomes strongly unfavourable due to the formation of a bond between the two free Ni sites.This distance can be related to the metallic, crystalline fcc Ni phase, where a lattice constant of a = 3.524 Å 31

Conclusions
The electronic and magnetic properties of the flexible MOF DUT-8(Ni) have been investigated by means of DFT.For a qualitative and quantitative analysis of the magnetic behaviour the energy difference between a low-spin and a high-spin coupling has been calculated.The DFT results were used as parameters for the calculation of the coupling constant as derived from the Heisenberg-Dirac-Van Vleck Hamiltonian.It has been found that the Ni-N interaction is most important for the magnetism within the systems.In both crystalline structures (DUT-8(Ni) open /DUT-8(Ni) closed ) a low-spin ordering (S = 0) has been found to be the most favoured one.This result is confirmed by the supercell calculations, where a variety of different couplings between multiple unit cells has been analyzed, as well as from calculations including the vdW interaction.Several molecular model systems suggest a possible switch of the magnetic ground state to a high-spin character (S a 0) of the coupling by certain changes in the structure.Especially the removal of the NH 3 groups can be related to the removal of a dabco linker unit within the crystalline systems.All alterations can and should be introduced to the periodic systems both theoretically as well as experimentally to find out whether the modifications influence the magnetic ground state as indicated by the model systems.
The tetragonal unit cell (a = b E 18.679 Å and c E 9.613 Å) for DUT-8(Ni) open contains 132 atoms and two Ni dimers.The unit cell of DUT-8(Ni) closed contains 66 atoms and only one Ni dimer.It is described by a E 6.947 Å, b E 8.180 Å, c E 12.172 Å and a E 91.231, b E 103.871, g E 104.551.
The DUT-8(Ni) open supercell (Fig. 2a) has been constructed by replicating two unit cells in the c direction.For a comparable size of the DUT-8(Ni) closed supercell (Fig. 2b), four unit cells were taken into account, two in the a and two in b direction each.This leads to 264 atoms and four magnetic dimer units (I, II, III and IV) for each supercell.There is a major difference between the supercells of DUT-8(Ni) open and DUT-8(Ni) closed .This difference is connected with the organic linkers.In the DUT-8(Ni) closed supercell they are only introduced by the periodic boundary conditions.As the magnetic interaction along the organic linkers is small (see Section 3.2), an explicit consideration in the DUT-8(Ni) closed supercell is not mandatory.On the other hand those linkers are contained in the DUT-8(Ni) open supercell according to the elementary cell of DUT-8(Ni) open .

Fig. 3
Fig.3Possible couplings within the supercells leading to different global magnetic orderings.The local interaction at each magnetic dimer unit, indicated by spheres, is assumed to be high-spin.
The density of states (DOS) for DUT-8(Ni) open and DUT-8(Ni) closed around the Fermi energy E F is shown in Fig. 5.The electronic densities (n) have been decomposed into spin up (m) and spin down (k) contributions to get a first insight into the spin polarization (n m Àn k ) and the magnetic behaviour.The partial DOS (PDOS) for all Ni d-states with the respective spins are also shown next to the total DOS.Fig. 5(a), representing DUT-8(Ni) open , shows a very similar contribution of the m and k densities over the given energy range (E = E F AE 4 eV).This can be attributed to the symmetry of DUT-8(Ni) open .The PDOS for the Ni d-states shows that one Ni per dimer (dimer 1: Ni 1,2 and dimer 2: Ni 3,4 ) carries one specific spin (m) while the respective other Ni carries the opposite spin (k) at the same energy.
8(Ni) open are adjusted, meaning that the given DE LS-HS has been divided by 2 in order to account to one Ni dimer and to be comparable to DUT-8(Ni) closed ).The calculations give DE LS-HS,DUT-8(Ni) open ,QE = À199 meV DE LS-HS,DUT-8(Ni) closed ,QE = À81 meV and indicate a strong low-spin (S = 0) ordering in both crystalline structures, confirming the discussions based on the DOS.The calculated energy differences using GPAW are DE LS-HS,DUT-8(Ni) open ,GPAW = À206 meV DE LS-HS,DUT-8(Ni) closed ,GPAW = À77 meV.which is in good agreement with the results of the QE calculations.Using eqn (2) and the given DE LS-HS from the DFT

Fig. 5
Fig. 5 Density of states, splitted into spin up (m) and spin down (k) densities, for (a) DUT-8(Ni) open and (b) DUT-8(Ni) closed .States with energies around the Fermi level E = E F AE 4 eV are displayed.
and the low-spin state.This can be explained by the Ni-N distances.For DUT-8(Ni) open one has d Ni 1&2 -N E 2.11 Å for both Ni while the distances for DUT-8(Ni) closed are d Ni 1 -N E 1.91 Å and d Ni 2 -N E 2.09 Å, respectively.Thus the interaction between Ni 1 and the neighbouring N is stronger than the respective interaction for Ni 2 .As a consequence the Ni-Ni interaction is weaker than for the DUT-8(Ni) open models which leads to smaller DE LS-HS .Especially for DUT-8(Ni) closed (M4 and M6 models) the resulting values for DE LS-HS become too small to be evaluated.This shows again the importance of the Ni-N interaction on the magnetic ground state.To simulate the removal of a dabco linker unit within the crystalline structures another model system is introduced.The system consists of two M3 molecules which are facing each other at the free Ni sites.The results of this approach are displayed in Fig.6.Only the DUT-8(Ni) open models were taken into account.Similar effects should be observed for the respective DUT-8(Ni) closed models.Starting from a distance between the free Ni sites of d Ni-Ni = 2.3 Å (binding situation), it has been steadily increased by 0.1 Å steps to d Ni-Ni = 9.3 Å (two isolated systems) is found.Thus the Ni-Ni distance accounts to d Ni-Ni ¼ a ffiffi ffi 2 p ¼ 2:492 Å, in agreement with the trend seen in Fig. 6.

Table 1
Changes in the initial model systems (M1) to gain a set of different models.Abbreviations for the respective model systems are added (M2-M10)

Table 2
Results for supercell calculations onDUT-8(Ni)open and DUT-8(Ni) closed featuring different magnetic orderings.Magnetic dimer units I, II, III and IV as indicated in Fig.2.The m and k indicate the local spins at the Ni within the magnetic dimer units.M tot ¼ Ð ðn " À n # Þd 3 r is the total magnetization of the system, M abs ¼ Ð jn " À n # jd 3 r is the corresponding absolute magnetization and DE LS(local) = E LS(local) À E ref is the energy difference between the local low-spin and the considered magnetic ordering This journal is © the Owner Societies 2015 Phys.Chem.Chem.Phys., 2015, 17, 17122--17129 | 17127

Table 3
Results of spin-polarized calculations on the model systems leading to the difference in total energy between the low-spin state and the high-spin state DE = E LS À E HS and the corresponding coupling constant J