Chemical Switching of the Magnetic Coupling in a MnPc Dimer by Means of Chemisorption and Axial Ligands

: We present an ab initio density functional theory study of the magnetic properties of manganese phthalocyanine dimers, where we focus on the magnetic coupling between the Mn centers and on how it is a ﬀ ected by external factors like chemisorption or atomic axial ligands. We have studied several di ﬀ erent con ﬁ gurations for the gas phase dimers, which resulted in ferromagnetic couplings of di ﬀ erent magnitudes. For the bare dimer we ﬁ nd a signi ﬁ cant ferromagnetic coupling between the Mn centers, which decreases by about 20% when a H atom is adsorbed on one of the Mn atoms and is reduced to about 7% when a Cl atom is adsorbed. The magnetic coupling is almost fully quenched when the dimer, bare or with the H ligand, is deposited on the ferromagnetic substrate Co(001). Our calculations indicate that the coupling between the two Mn atoms principally occurs via a superexchange interaction along two possible paths within a Mn − N − Mn − N four-atom loop. When these electrons get involved in chemical bonding outside the dimer itself, an appreciable alteration of the overlap between Mn and N molecular orbitals along the loop occurs, and consequently, the magnetic interaction between the Mn centers varies. We show that this is re ﬂ ected by the electronic structure of the dimer in various con ﬁ gurations and is also visible in the structure of the atomic loop. The chemical tuning of the magnetic coupling is highly relevant for the design of nanodevices like molecular spin valves, where the molecules need to be anchored to a support.


■ INTRODUCTION
A key objective of molecular electronics and spintronics (spinbased electronics) is the downscaling of the electronic components, which would bring up remarkable benefits like increase in magnetic storage density and reduction of power consumption. In this context, the possibility to use organic molecular materials of low production cost has become highly appealing. The molecules of the phthalocyanine family (Pc) have been studied quite extensively in the emerging field of organic spintronics. 1−4 Among the milestones in this field are the possibility of spin-polarized injection and transport of electrons through organic semiconductors, shown by Dediu and colleagues, 5 and the giant magnetoresistance observed for a single H 2 Pc 1 and CoPc. 6 In addition, an organic spin valve device based on immobilized layers of CuPc (spacer layer) and aligned MnPc (or NiPc) as spin-injection and spin-detection layers was fabricated and studied experimentally by Banerjee and co-workers. 7 Recently, also supermolecular nano and lowdimensional structures like molecular chains and dimers have been addressed, with an emphasis on their structures 8 and magnetic properties. 9,10 These systems are promising novel building blocks for molecular size devices as well as model systems to understand the magnetic interactions and spin transfer mechanisms between the molecules or between the molecules and substrates. An important issue in this field is to find feasible mechanisms to manipulate the molecular spin and magnetic moments in a reproducible fashion.
In the present work, we target dimers of the MnPc molecule, which is among the most studied phthalocyanines in molecular spintronics, due to the known magnetic interactions in MnPcbased bulk materials. 11,12 Interestingly, the deposition of MnPc molecules on the ferromagnetic Co(001) surface has shown to generate a so-called spinterface, 13,14 i.e., an active heterogeneous interface for spin transport. 3 Beyond that, MnPc and FePc coadsorbed in the form of a checkerboard table on Au(111) were reported to be antiferromagnetically coupled through the Ruderman−Kittel−Kasuya−Yosida (RKKY) exchange interaction via the substrate electronic states. 15 MnPc, shown in Figure 1A, consists of a single Mn center surrounded by an organic ring, representing a sort of archetype of single molecule magnet. The molecule has a spin 3/2 and a D 2h symmetry. The magnetic properties of specific phases of MnPc films were described already some decades ago. Both ferromagnetic and antiferromagnetic couplings were reported for respectively the β and α phases in bulk MnPc films, 11,12 showing that the magnetic interaction between the molecules is directly linked to their reciprocal arrangement. 12 Because of this, it is essential to unveil how the interplay between the electronic and geometric structure shapes the magnetic properties to engineer the future molecule-scale devices.
We have first studied the magnetic properties of dimers of MnPc in the gas phase, considering several possible molecular configurations by means of ab initio density functional theory (DFT) including dispersion interactions. We have then focused on the dimer with the configuration of lowest energy and studied its adsorption on a Co(001) substrate, both bare and with atomic axial ligands such as H and Cl. It has been shown in previous joint theory/experiment studies how these theoretical methods well reproduce the experimental findings for the adsorbate MnPc/Co(100) system in terms of adsorption distances and magnetic properties. 3,16 Previous theoretical studies have reported dimers of phthalocyanines in geometrical configurations that are derived from the various 3D bulk stackings, 8 where the magnetic interactions are directly associated with those in the thick films. We have here analyzed the connection between the magnetic coupling and the inter-and intramolecular chemical bondings that occur in the dimer. The Mn atoms at the center of the molecules are coupled via superexchange interaction, which unrolls via a N atom directly bonded to one of the Mn atoms and situated in the face of the Mn of the opposing molecule, via a mechanism analogous to the one described for MnPc films and bulk materials. 11,12 This type of superexchange path was also experimentally investigated by Chen et al. 17 by means of spinflip electron tunneling spectroscopy between CoPc molecules adsorbed on a Pb(111) surface, where the CoPc molecules assume a reciprocal configuration analogous to the one of the dimers in our study. We observe that when the Mn atoms take part in intermolecular chemical bonds, the magnetic coupling is markedly reduced. This happens in all cases when the dimer is chemisorbed on a Co surface via one of the MnPc rings or when a H or a Cl atom is axially adsorbed on the Mn center on one side of the dimer. In these cases, the chemical bonding scheme of the Mn atoms is changed, as can be seen in the modifications of the inter-and intramolecular distances between the Mn and N atoms involved.

■ THEORETICAL METHODS
The present study is based on ab initio DFT calculations using the VASP 18 code with the Perdew, Burke, and Ernzerhof (PBE) exchange correlation functional 19,20 and the projector augmented wave method. 21 The plane wave cutoff was 400 eV, and only the Γ-point k-mesh was used. An effective Hubbard term U eff of 4 eV was chosen to describe the Mn 3d orbitals, following Dudarev's method. 22 The value for the Hubbard term has been chosen according to the findings in our previous study. 23 This value has been proven to properly reproduce experimental data, for example, photoemission spectra. 23 The long-range dispersion forces were described with Grimme's second method (D2). 24 For comparison, some tests were performed with the (D3) 25 method showing no differences in the relevant properties. To describe the Co(001) substrate, we have used a three-layer fcc Co film as in our previous works with FePc on Co(001), 16

■ RESULTS AND DISCUSSION
We have first analyzed different possible geometries for the MnPc dimer based on possible arrangements of multilayers of metal Pc molecules. We will briefly discuss these structures which are shown in Figure S1 of the Supporting Information. Starting from the α and β polymorphs, we have built the corresponding models α× with a tilt angle of about 63°, the α+ with a tilt angle of 59°( Figure 1B), and the β type with a tilt angle of 43°, where the tilt angle is the angle between the line joining the two Mn atoms and a line orthogonal to the plain of the molecules. In the α+ and β configurations, the two molecules are shifted with respect to each other along an axis joining the Mn center to an isoindole N atom (N iso ), while in the α polymorph, the reciprocal shift occurs along a line joining the Mn center to an aza bridge N atom (N aza ) (see Figure 1A). The α+ configuration was previously proposed in a theoretical study of MnPc/F 16 CoPc dimers, 8 and the magnetic properties of MnPc and CuPc in this configuration were studied by Wu et al. 9 Furthermore, this configuration was experimentally observed in a recent STM investigation of the growth of FePc multilayer films on Au(111), as the reciprocal arrangement of the molecules between the first and second FePc layers. 27 In addition to these three models, we have studied a dimer with two MnPc lying on top of each other and rotated by 45°, a structure that was for example observed as the reciprocal geometry of two CoPc in the first two layers of films of CoPc grown on a Pb(111) substrate. 17 After a geometry relaxation of all the structures, we find that α+ is the one with lowest energy, as can be seen from Table 1 with the total energies, and it is the structure the rest of this study focuses on.
To get an estimation of the magnetic coupling between the MnPc molecules in the dimers, we computed the energy difference between the ferromagnetic (FM) and the antiferromagnetic (AFM) arrangements, ΔE FA , for all the structures given in Table 1. After geometry relaxation we obtain in all cases ferromagnetic coupling between the Mn The Journal of Physical Chemistry C pubs.acs.org/JPCC Article centers (corresponding to a positive value of ΔE FA ), but with different magnitudes in the various configurations. The top 45°d imer has the strongest coupling but has a higher total energy. In the dimers, both the α+ and the α× actually have a ferromagnetic coupling according to our calculations unlike what is reported for the α polymorphs. 11,12 ΔE FA is higher in the α+ and in the β dimers than in the α×. We attribute this to the different mechanisms that govern the magnetic coupling between the two Mn atoms which depend on their relative positions. In fact, different tilt angles in the dimers, resulting from different shifts between the planes of the two molecules, are generally expected to affect the superexchange mechanism, since the relative orientation of the molecules influences the possible overlap between different molecular orbitals. In the α+ dimer, each Mn atom is located precisely on top of a N iso of the other MnPc. From the side view of Figure 1B, interestingly, it is noticeable that the MnPcs in the dimer are indeed influenced by the presence of the opposite phthalocyanine, since they are not as flat as the single molecule in gas phase ( Figure 1A). The two central Mn atoms are 3.20 Å apart and slightly protruding from the plane of the molecule toward each other. For MnPc in gas phase, the computed bonding distances between the Mn and the N iso atoms amount to 1.97 and 1.96 Å in the two orthogonal directions because of the D 2h symmetry of the molecule induced by Jahn−Teller distortions. 23,28 When the dimer is formed, the D 2h symmetry is lifted. The distance between the Mn and the N iso that is facing the Mn of the other molecule is now 1.98 Å, slightly elongated with respect to the gas phase bond lengths, while the bonding distance with the three remaining N iso is 1.96 Å.
In the Pc films, the magnetic coupling between the metal atoms is generally explained by an intermolecular superexchange interaction that connects two successive Mn centers along the chains through the N atoms of the closest molecules. The interaction depends on the overlap of the 3d orbitals of Mn with the π orbitals centered on the N atoms in the organic shell. 12 As in the α and β phases, the relevant orbitals extending between the two molecules are the out-of-plane Mn 3d z 2 and the 3d π (i.e., the 3d xz and 3d yz ) and the N iso 2p z . One should point out that these electrons actually participate in molecular orbitals that involve also other atoms of the MnPc (C, N, and Mn). Furthermore, more than one orbital combination may contribute to the superexchange interaction, generating a more complex process than the one usually described. In a simplified scheme, the ferromagnetic superexchange coupling in the β phase is assumed to occur mainly via the overlapping of the Mn 3d z 2 with the N (π), while the antiferromagnetic coupling in the α phase via the overlapping of the Mn d π with the N (π) orbitals. In the ferromagnetic α+ configuration, the atoms involved in the superexchange path form a loop that develops along the Mn centers and two of the N iso atoms located on the opposite MnPc, while in the β configuration, shown in Figure  S1, the superexchange path develops instead along the Mn atoms and the N aza atoms. These four atoms in the α+ configuration form a rectangular loop that is highlighted in Figure 2 for three of the different structures based on the α+ dimer considered in this study and that will be discussed in the next paragraphs, namely (MnPc) 2 in gas phase ( Figure 2A) and adsorbed on the Co substrate ( Figure 2B), and (MnPc) 2 with a Cl axial ligand ( Figure 2C). The two Mn atoms in the loop are indicated as Mn1 and Mn2 and the N iso in front as N1 and N2, where the indices 1 and 2 identify the two molecules, 1 being the upper MnPc, on which the axial ligands are adsorbed, and 2 being the index of the molecule directly adsorbed on Co(001). In the bare dimer in gas phase, the superexchange loop has a regular rectangular shape with intramolecular distances between the Mn atoms and the N iso of 1.98 Å and intermolecular distance between the Mn atoms and the N iso of 2.74 Å. In this configuration, we obtain the highest ferromagnetic coupling, with a ΔE FA value of 257 meV. The atomic distances in the four-atom loop for all the configurations studied with the respective ΔE FA are reported in Table 2.
The next step was to simulate the deposition of the α+ dimer on a Co(001) surface. In our previous works regarding the adsorption of FePc on Co(001), we found a ground state corresponding to the molecule adsorbed on a top site, that is, with the Fe atom on top of a Co atom and with the Mn−N iso main axes of the molecule oriented along the [100] crystallographic direction of the surface. 16,26 The adsorption of a single MnPc on the same surface gives analogous outcomes for structures and adsorption sites as for the FePc, as reported in Table S1. Following these results, we have placed the whole (MnPc) 2 dimer on the same top Co site, and the  The Journal of Physical Chemistry C pubs.acs.org/JPCC Article relaxed structure is displayed in Figure 1C. The chemisorption casts the dimer into a heterogeneous environment with one phthalocyanine (MnPc2) in contact with the metal and the other one (MnPc1) facing vacuum. Hence, the structural symmetry between the two molecules is lifted, with an impact on the overall geometry and properties of the dimer. A sideview inspection of the adsorbed dimer ( Figure 1C) shows that MnPc2 is subject to noticeable distortions of the planar symmetry, while MnPc1 maintains a flatter geometry, akin to the gas phase molecule. The distance between the Mn atoms (Mn2 in the dimer) and the Co atom beneath is 2.43 Å for the single molecule and 2.36 Å for the adsorbed dimer, while the distance between the two Mn atoms for the adsorbed dimer is 3.32 Å, larger by 0.12 Å than in the gas phase. These distances are similar to those obtained by means of a X-rays standing wave experiment of MnPc adsorbed on Cu(001) which gave an adsorption distance of 2.240 ± 0.045 Å for the Mn to the top layer of the surface. 29 When the MnPc dimer is adsorbed on Co(001), we still obtain ferromagnetic coupling between the two molecules and between the dimer and the substrate (see Table 2). However, we find that the magnetic coupling between the two MnPc's is dramatically decreased, with ΔE FA reduced to 15 meV. In practice, we find that the deposition of the dimer on the Co substrate severely weakens the magnetic coupling between the molecules. This result may imply that in a magnetic nanodevice based on a chain or on a multilayer of MnPc the first layer would not be magnetically coupled to the remaining layers. However, the magnetic moments on the Mn centers are only very little affected by the adsorption, while for similar compounds, like for CoPc on Pb(111), 30 the first molecular layer acts as a spin insulating buffer. The quenching of the coupling is reflected by the distortion of the superexchange path presented in Figure 2. In fact, the chemical bonding formed between the Mn1 center and the Co atom beneath alters the bond lengths between the N and Mn atoms reported in Table 2. In the slightly corrugated geometry of the chemisorbed MnPc2 ( Figure 1) the N2 atom is lifted toward Mn1, resulting in a shorter Mn1N2 distance (2.46 Å) and, on the other side, inducing a longer Mn2N1 distance (2.93 Å). This can be seen from the slightly irregular loop in Figure  2B. As is also the case for the MnPc dimer with H and Cl ligands, the increase of the intermolecular Mn−N bond length between the two molecules affects the overlap of the orbitals along the loop and is correlated to the quenching of the coupling.
Experimental and theoretical studies of the axial adsorption of small molecules on several 3d metal Pc's have demonstrated how the ligands can influence the magnetic moment and spin state of the metal macrocycles. 2,31,32 Here we focus on the effect of atomic H and Cl axial ligands and on how the coupling is additionally affected by the deposition on Co(001). When a H atom is adsorbed on one side of (MnPc) 2 , the overall geometric structure of the dimer is only marginally changed with respect to the bare dimer ( Figure S2A). After geometry relaxation, the distance between Mn1 and Mn2 is 3.26 Å, which is 0.06 Å larger than in the bare dimer, and the Mn−N iso distances are between 1.95 and 1.98 Å: in MnPc1, bonded to H, the Mn−N iso distances are generally shorter by just about 0.01 Å. The two MnPc have ferromagnetic coupling with a ΔE FA of 200 meV, which is 20% lower than in the bare dimer. In comparison to the gas phase dimer, the sides of the superexchange loop are only weakly distorted, with the Mn2− N1 distance stretched from 2.74 to 2.83 Å and the Mn1−N2 distance shortened from 2.74 to 2.70 Å (see Table 2). The deposition of the dimer on the Co(001) surface on the top site ( Figure S3A) leads instead to considerable changes. We obtain a quenching of the magnetic coupling with a ΔE FA of −6 meV, where the minus sign indicates an antiferromagnetic character. In this case, we find a larger distance between the Mn centers of 3.41 Å. The Mn1−N2 bond length is 2.90 Å and the N1− Mn1 bond length is 2.03 Å, thus stretching two sides of the superexchange loop with respect to the bare dimer in the gas phase.
When a Cl ligand is adsorbed on MnPc1 ( Figure S2B), Mn1 is pulled toward the Cl atom, which lies 2.33 Å above it. The Mn1−Mn2 distance is increased in the presence of Cl to 3.62 Å in the gas phase and to 3.71 Å when adsorbed on the Co substrate ( Figure S3B); in addition, all the bond lengths between the atoms in the superexchange loop are increased as reported in Table 2. The calculations show that in this case the Cl ligand brings about a drastic reduction of the coupling: the calculated ΔE FA is 18 meV in the gas phase and 15 meV when the dimer is adsorbed on Co (see Table 2). The coupling is ferromagnetic in both cases.
The distance between the Mn1−Mn2 atoms reported in Table 2 is not a major factor in determining the strength of the magnetic interaction: for example, it only varies by 0.06 Å between (MnPc) 2 and H-(MnPc) 2 , which have ΔE FA differing by 20% from each other, and again it varies by 0.06 Å between H-(MnPc) 2 and (MnPc) 2 adsorbed on Co, where the coupling is reduced to 15 meV. Rather, the coupling depends mostly on a more complex interplay between the orbital overlaps along the superexchange path. Summarizing the results of Table 2 for all the systems studied, we can point out how the magnetic coupling is considerable for only two structures, namely, for the bare dimer and for the dimer with the H ligand in gas phase. In these cases, the bond lengths within the superexchange loop are very similar, with intramolecular Mn−N bond lengths between 1.97 and 1.98 Å and intermolecular Mn−N bond lengths between 2.70 and 2.83 Å. When the bare (MnPc) 2 is adsorbed on Co(001), ΔE FA goes down to 15 meV: here the superexchange loop is distorted on one side, with the bond length Mn1−N2 increased to 2.94 Å, which is 0.2 Å more than in gas phase. When H-(MnPc) 2 is adsorbed  Table 3 shows the computed magnetic moments on the Mn, C, and N atoms and on the ligands for the monomer and for all the examined configurations of the α+ dimers.
In (MnPc) 2 in the gas phase, in each molecule, the major part of the magnetic moment resides on the Mn atom (3.6 μ B ), while about 12% of this moment with opposite sign (−0.46 μ B ) is distributed among the N atoms and the C atoms bonded to them in the organic ring, as shown by the computed isosurface magnetization density for (MnPc) 2 in Figure 3A. This distribution of the magnetic moments is very similar to that of the single molecule in gas phase with 3.5 μ B on Mn and −0.44 μ B on the N and C atoms in the organic ring.
When MnPc is adsorbed on Co(001), both as single molecule and as a dimer, we obtain a slightly smaller moment on the Mn atom directly in contact with the surface, varying between 3.41 μ B for MnPc and 3.46 for (MnPc) 2 , as can be seen in Table 3. At the same time, we see an increase in the moments of the C atoms in the benzene rings that bond to the Co surface by about 0.1 μ B (see Figure 3B). On the other hand, the total moment of the N atoms decreases by about 1.1 μ B , since here the moment in the ring is concentrated in the N iso with negligible contributions from the N aza . These results are similar to the variation of magnetic moments for MnPc adsorbed on Cu(001) computed with the same method by Alouani and co-workers, 29 where between gas phase and adsorption on surface, a variation in the Mn moment from 3.46 to 3.34 μ B was found, with a change of the total C moment from −0.15 to −0.27 μ B and a very small change of the total N moment from −0.15 to −0.14 μ B . The small differences observed in the present work could be attributed not only to the different surfaces but also to the different adsorption sites considered in the two works, since this could affect the bonding and therefore the charge transfer of the N and C ligands to the substrate and finally their magnetic moments.
The axial ligands mainly modify the distribution of the moments in the central part of MnPc1 where they are adsorbed. In H-(MnPc) 2 , the moment on Mn1 is 2.95 μ B with −0.50 μ B on the ring (about 17% of the moment on Mn1 but with opposite sign). The H ligand brings about a decrease of the moment on Mn1 accompanied by an increase of the moment on the C atoms. The Cl ligand, on the contrary, induces an increase of the magnetic moment on Mn1 to 3.74 μ B . Now a moment of −0.17 μ B (5% of the moment on Mn1) is located on the ring and mostly concentrated on N iso , as can be seen from Figure 3C.
WhenH-(MnPc) 2 is adsorbed, the moment in Mn1 further decreases to −2.80 μ B , and it goes to 0.44 μ B on the ring. In Cl-(MnPc) 2 the moment stays at 3.74 μ B on Mn1, and it is lowered to −0.11 in the ring. The antiferromagnetic coupling observed in all these cases between the central metal and the organic ring was also been observed for several metal Pc's deposited on surfaces, 29,33 while for some Pc-based 3D materials, it was found to be ferromagnetic, for example, in FePc (CN) 2 ·2CHCl 3 crystals. 34 In addition, we have considered the magnetic interaction of the dimer to the surface, which takes place through MnPc2. We have obtained a ferromagnetic coupling of the Mn2 atom with the Co atoms as well as an antiferromagnetic coupling of the C atoms in MnPc2 to the surface. The same magnetic configuration was obtained for the single molecule and for the dimers with axial ligands. The bonding of the molecule to the surface occurs through the direct interaction of especially Mn and C atoms in the benzene rings to the underlying Co atoms, and this is coupled to variations in magnetic moment on both the MnPc and Co(001) surface. This results in a sort of  The Journal of Physical Chemistry C pubs.acs.org/JPCC Article footprint left by the molecule on the underlying surface, as can be seen in Figure 4, where the magnetic moment of the top Co atoms on the surface are shown. A profile highlighting the position of the adsorbed molecule is also sketched. The magnetic moments of the Co atoms vary between about 1.6 and 1.9 μ B , and the strongest variations are given by those Co atoms lying below the C atoms in the benzene rings, suggesting the decisive role played by these atoms in the bonding of MnPc to the surface. We had previously found the same type of patterning for the adsorption of FePc on the same substrate. 16 We have finally examined the interplay between the chemical and geometrical configurations and the electronic structure of the MnPc molecules in the dimer and how this is related to the magnetic coupling. The partial density of states (pDOS) of the single molecule and of the MnPc dimer in gas phase and on Co(001) are presented in Figure 5A−C. In Figure 5A, the pDOS of the single molecule is illustrated, with the Mn 3d orbitals, 3d z 2 , 3d π , 3d xy , and 3d x 2 −y 2 , and the N iso 2p z . For the single MnPc in gas phase we obtain a 4 E g ground state where the 3d z 2 and the 3d xy are half occupied, and the 3d π contains three electrons (see Figure 5A), in agreement with previous DFT studies 23,35 and experimental results. 36,37 MnPc's maintain the 4 E g electronic configuration in the α+ dimer and even when they are adsorbed on Co(001) and with the axial ligands ( Figure 5B,C). Because the dimer is symmetric for the two molecules, only the pDOS for one of the two is shown in Figure 5B. First of all, one can notice how the M n 3d z 2 orbital in the dimer is split into several components of lower intensity (black lines in Figure 5). This is evident when comparing the two high-intensity 3d z 2 peaks of the single MnPc at −4.3 and 1.5 eV in Figure 5A, with the several low-intensity 3d z 2 peaks in the regions between −5.4 and −3.8 eV and between 0.8 and 2.6 eV in the gas phase dimer in Figure 5B, and is a result of the interaction between the two molecules. Also, the 3d π states (gray curves) undergo some minor changes in energy position and especially in intensity, although they do not further split into several smaller components, giving a hint that the 3dπ states are less involved in the hybridization between the molecules. In all the diagrams, also the partial DOS of N iso 2p z is plotted, for the specific N iso which is located opposite to a Mn atom of the other molecule; this pDOS is very weak compared to that of the Mn 3d electrons, and its intensity is therefore multiplied by a factor 4. In (MnPc) 2 the N iso 2p z pDOS overlaps the Mn 3d orbitals in several energy intervals, like for example between −5 and −4 eV, at about 0.5 and 2 eV in Figure 5B, suggesting hybridization with the Mn 3d electrons. In this dimer as we   When the MnPc dimer is adsorbed on Co(100), the lower molecule (MnPc2) is strongly bonded to the substrate. The Mn center on a top site bonds to the Co atom above through the out of plane 3d z 2 and 3d π electrons. In this configuration, the MnPc is ferromagnetically coupled to the Co surface, in analogy to FePc on Co(100), as we had reported in our previous studies 16,26 and as was also found experimentally for example for the adsorbed systems FePc, CoPc, NiPc, and CuPc on Ag(100), FePc, CoPc, and CuPc on Co(001), 38 and MnPc on Co(001). 3 The pDOS curves of MnPc2 in Figure 5C show that both the Mn 3d z 2 and 3d π electrons are strongly affected by the adsorption, with the formation of several peaks with reduced intensity. The involvement of these electrons into the bonding to the surface corresponds to a somewhat weakening of the bonding between MnPc1 and MnPc2, which could also be deduced by the flatter geometrical structure of MnPc1, and is also responsible for the drastic decrease of the magnetic coupling between the molecules with a ΔE FA of 15 meV. Figure 6A shows the pDOS of (MnPc) 2 with adsorbed H. Although the ferromagnetic coupling is only reduced by 20% from the bare dimer, some differences in the DOS of MnPc1 are evident. For example in the Mn1 3d z 2 pDOS we observe a different splitting of the peaks and a new 3d z 2 feature at 0.57 eV. However, the pDOS of the MnPc2 and the N iso 2p z are mostly analogous to the ones of the bare dimer. Figure 6B reports the pDOS of the dimer with the Cl ligand. In this case, a significant weakening of the magnetic coupling between the Mn centers occurs, as when the dimer is adsorbed on the Co(001). The Cl is adsorbed on MnPc1, forming chemical bonds between the M n 3d z 2 and 3d π electrons to the Cl 2p (see Figure 6B between −3 and −1 eV). These new bonds weaken the interaction between the two MnPc's significantly, which are now at 3.71 Å from each other. This effect is also visible in the pDOS of MnPc2 of Figure 5C, which resembles very closely the one of the gas phase MnPc, with distinct peaks for the M n 3d z 2 and 3d π electrons. On the other side, the pDOS of MnPc1 is highly affected by the bonding to Cl.
For all these gas phase and adsorbed systems, the pDOSs provide further evidence of the role of the M n 3d z 2 and 3d π in creating the superexchange ferromagnetic interaction between the two Mn moments, since in all cases when these orbitals are involved in other chemical bondings, the superexchange path is disturbed and the coupling gets affected.

■ CONCLUSIONS
To summarize, the analysis of the interplay between the reciprocal arrangement of the molecules and the magnetic interaction in the α+ dimer confirms that the magnetic coupling between the two MnPc takes principally place via superexchange interaction, exploiting one of the two paths in a loop entangling the Mn centers and two isoindole N atoms. The path is conceptually analogous to the one pointed out for the α and β bulk phase, although in the dimers only ferromagnetic coupling is obtained. In the dimer, it is possible to generate distortions in the atomic loop, for example, by chemisorption on a substrate or by adsorption of axial ligands. We have highlighted how the coupling between the molecules can be tuned by ligand adsorption, and how it can be influenced by the adsorption of the dimer on the FM Co(001). In addition, MnPc has a ferromagnetic coupling via the Mn atom to the underlying surface, but the C and N atoms have an antiferromagnetic coupling to the same.