Magnetic Anisotropy of Homo- and Heteronuclear Terbium(III) and Dysprosium(III) Trisphthalocyaninates Derived from Paramagnetic 1H-NMR Investigation

1H-NMR spectroscopy of lanthanide complexes is a powerful tool for deriving spectral–structural correlations, which provide a clear link between the symmetry of the coordination environment of paramagnetic metal centers and their magnetic properties. In this work, we have first synthesized a series of homo- (M = M* = Dy) and heteronuclear (M ≠ M* = Dy/Y and Dy/Tb) triple-decker complexes [(BuO)8Pc]M[(BuO)8Pc]M*[(15C5)4Pc], where BuO- and 15C5- are, respectively, butoxy and 15-crown-5 substituents on phthalocyanine (Pc) ligands. We provide an algorithmic approach to assigning the 1H-NMR spectra of these complexes and extracting the axial component of the magnetic susceptibility tensor, χax . We show how this term is related to the nature of the lanthanide ion and the shape of its coordination polyhedron, providing an experimental basis for further theoretical interpretation of the revealed correlations.


Introduction
Single-molecule magnetism (SMM) is one of the most intriguing manifestations of the anisotropic coordination environment of paramagnetic lanthanides in sandwich complexes with tetrapyrrolic ligands, which was actually first described for this family of elements using the example of terbium(III) bisphthalocyaninate, Tb(Pc) 2 .Since the seminal report by Ishikawa et al. two decades ago [1], numerous examples of sandwiches containing mainly Tb 3+ and Dy 3+ metal centers have been synthesized and extensively studied in terms of magnetic relaxation dynamics, and the potential for fabrication of spintronic devices based on them has been clearly demonstrated [2][3][4].
Among such complexes, we can distinguish a prominent class of triple-decker phthalocyaninates, in which two metal ions bind three ligands, providing a wide range of combinations of different macrocycles and metal centers to achieve desired functional characteristics [5][6][7].While the synthesis and properties of a large number of heteroleptic triple-deckers have been reported [8], the family of heteronuclear derivatives containing two different rare-earth elements (REEs) is still relatively underexplored.At the same time, such a combination provides unique possibilities for studying the subtle effects of intramolecular interactions between paramagnetic ions separated by a distance of only 3.4-3.5Å.
The first vivid contribution to such studies was made by Ishikawa et al., who synthesized three series of heteroleptic triple-deckers [(BuO) 8 Pc]M*(Pc)M(Pc), where [(BuO) 8 Pc] 2− intramolecular interactions between paramagnetic ions separated by a distance of only 3.4-3.5Å.
The first vivid contribution to such studies was made by Ishikawa et al., who synthesized three series of heteroleptic triple-deckers [(BuO)8Pc]M*(Pc)M(Pc), where [(BuO)8Pc] 2− and (Pc) 2− were octa-butoxy-substituted and unsubstituted phthalocyaninato ligands, respectively [9,10].One series consisted of homonuclear complexes with M* = M = Tb 3+ , Dy 3+ , Ho 3+ , Er 3+ , Tm 3+ and Yb 3+ , and the other two families were isomeric heteronuclear complexes containing a diamagnetic Y 3+ ion at either the M* or M site, while the other site was occupied by one of the listed lanthanides.A comparison of the temperature dependence of magnetic susceptibilities between homo-and heteronuclear complexes has provided the first detection and characterization of the dipolar f-f interactions between Ln 3+ ions.
Later, the influence of these interactions on the behavior of SMMs based on tripledeckers was revealed.It was demonstrated that heteronuclear complexes containing one dia-and one paramagnetic ion, Tb 3+ and Y 3+ , acted as field-induced SMMs, while the homonuclear di-terbium complexes showed slow magnetic relaxation even in a zero dc field.Such a behavior was attributed to the f-f interactions acting as an exchange bias that suppresses the undesired quantum tunneling mechanism of magnetic relaxation [11][12][13].This effect was first demonstrated on the examples of the aforementioned [(BuO)8Pc]M*(Pc)M(Pc) complexes [14], and a deeper insight into this phenomenon was gained with the 15-crown-5-substituted sandwiches [(15C5)4Pc]M*[(15C5)4Pc]M(Pc) reported by R. Holmberg et al. [15].In both cases, M* = M = Tb or M* ≠ M = Tb/Y.
The synthesis of heteronuclear trisphthalocyaninates is relatively straightforward.It implies raise-by-one-story reactions of bisphthalocyaninates bearing one metal center with monophthalocyaninates bearing another metal center.The latter precursors can be either synthesized beforehand [16][17][18] (Figure 1a), generated in situ from metal-free ligands [9] (Figure 1b) or by transmetalation of other complexes [15,19] (Figure 1c).Summary of synthetic raise-by-one-story approaches to heteronuclear trisphthalocyaninates via addition of a pre-synthesized monophthalocyaninate to a double-decker complex (a) [16]; generation of monophthalocyaninates in situ from a metal-free ligand (b) [9]; and transmetalation of dilithium phthalocyaninate (c) [20].Summary of synthetic raise-by-one-story approaches to heteronuclear trisphthalocyaninates via addition of a pre-synthesized monophthalocyaninate to a double-decker complex (a) [16]; generation of monophthalocyaninates in situ from a metal-free ligand (b) [9]; and transmetalation of dilithium phthalocyaninate (c) [20].
In general, these approaches can be used to obtain thermodynamically and kinetically stable heteronuclear triple-deckers with any pair of middle and late lanthanide ions, including complexes with two paramagnetic ions.Such an interplay between different combinations of nonequivalent magnetically active nuclei placed in well-defined positions can be used to study f -f interactions and provides a tool for fine-tuning the SMM behavior of such complexes.A rare example of such complexes-[(n-C 8 H 13 ) 8 Pc]Tb(Pc)Dy(Pc)-has been reported by Lan et al. [20].A weak but significant interaction between the two lanthanides was clearly detected, modifying the magnetic behavior of the single lanthanide as observed in the parent mononuclear complexes.
Trisphthalocyaninates also attract attention as conformationally flexible scaffolds where the twist angle between the ligands can be controlled by various inter-and intramolecular factors, including host-guest interactions [21,22], solvation [23,24] and redox transformations [25,26].The variation in the twist angle inevitably changes the symmetry of the coordination polyhedron of the lanthanide ion sandwiched between these ligands, which has a profound effect on its magnetic properties due to the influence of this symmetry on magnetic anisotropy [27] and energetic characteristics and preferable mechanisms of magnetic relaxation [11,28,29].
Recently, we reported a series of works on conformationally flexible triple-decker complexes composed of butoxy-and crown-substituted phthalocyaninato ligands [(BuO) 8 Pc]M [(BuO) 8 Pc]M*[(15C5) 4 Pc], including homonuclear complexes M = M* = Tb 3+ and Y 3+ [23] and their heteroleptic counterparts M ̸ = M* = Tb 3+ and Y 3+ [24].We showed that the pairwise conformations of the ligands are solvation-dependent: in halogenated alkanes, the pair of BuO-and 15C5-substituted ligands adopts the staggered conformation with a twist angle of ~45 • , while the pair of BuO-substituted ligands forms the gauche conformation with a twist angle of ~22 • (Figure 2).In turn, it provides metal centers M* and M with square-antiprismatic (SAP) and distorted prismatic (DP) coordination surroundings, respectively.In contrast, solvation with aromatic solvents leads to the stabilization of staggered conformations for both ligand pairs, so that both metal centers exist in SAP environments.The stabilization of certain conformers was explained by X-ray diffraction analysis of the solvates with dichloromethane and toluene, where weak contacts with the solvent molecules were found and analyzed in terms of the quantum theory of atoms in molecules [23].Moreover, these conformations were stable not only in the crystalline state but also in solution, which was confirmed by 1 H-NMR characterization of Tb 3+ complexes-their magnetic anisotropy χ ax strongly increases upon lowering the coordination symmetry from SAP to DP, as evidenced from the structure-based analysis of lanthanide-induced shifts.
In general, these approaches can be used to obtain thermodynamically and kinetically stable heteronuclear triple-deckers with any pair of middle and late lanthanide ions, including complexes with two paramagnetic ions.Such an interplay between different combinations of nonequivalent magnetically active nuclei placed in well-defined positions can be used to study f-f interactions and provides a tool for fine-tuning the SMM behavior of such complexes.A rare example of such complexes-[(n-C8H13)8Pc]Tb(Pc)Dy(Pc)-has been reported by Lan et al. [20].A weak but significant interaction between the two lanthanides was clearly detected, modifying the magnetic behavior of the single lanthanide as observed in the parent mononuclear complexes.
Trisphthalocyaninates also attract attention as conformationally flexible scaffolds where the twist angle between the ligands can be controlled by various inter-and intramolecular factors, including host-guest interactions [21,22], solvation [23,24] and redox transformations [25,26].The variation in the twist angle inevitably changes the symmetry of the coordination polyhedron of the lanthanide ion sandwiched between these ligands, which has a profound effect on its magnetic properties due to the influence of this symmetry on magnetic anisotropy [27] and energetic characteristics and preferable mechanisms of magnetic relaxation [11,28,29].
Recently, we reported a series of works on conformationally flexible triple-decker complexes composed of butoxy-and crown-substituted phthalocyaninato ligands [(BuO)8Pc]M[(BuO)8Pc]M*[(15C5)4Pc], including homonuclear complexes M = M* = Tb 3+ and Y 3+ [23] and their heteroleptic counterparts M ≠ M* = Tb 3+ and Y 3+ [24].We showed that the pairwise conformations of the ligands are solvation-dependent: in halogenated alkanes, the pair of BuO-and 15C5-substituted ligands adopts the staggered conformation with a twist angle of ~45°, while the pair of BuO-substituted ligands forms the gauche conformation with a twist angle of ~22° (Figure 2).In turn, it provides metal centers M* and M with square-antiprismatic (SAP) and distorted prismatic (DP) coordination surroundings, respectively.In contrast, solvation with aromatic solvents leads to the stabilization of staggered conformations for both ligand pairs, so that both metal centers exist in SAP environments.The stabilization of certain conformers was explained by X-ray diffraction analysis of the solvates with dichloromethane and toluene, where weak contacts with the solvent molecules were found and analyzed in terms of the quantum theory of atoms in molecules [23].Moreover, these conformations were stable not only in the crystalline state but also in solution, which was confirmed by 1 H-NMR characterization of Tb 3+ complexes-their magnetic anisotropy  strongly increases upon lowering the coordination symmetry from SAP to DP, as evidenced from the structure-based analysis of lanthanide-induced shifts.In the present work, we extended the developed synthetic and analytical approaches to Dy 3+ complexes, including a homonuclear triple-decker with M = M* = Dy 3+ and pairs of isomeric heteronuclear analogues with M/M* = Dy 3+ /Y 3+ and Y 3+ /Dy 3+ , as well as Dy 3+ /Tb 3+ and Tb 3+ /Dy 3+ .For the sake of brevity, we will henceforth refer to butoxy-and crown-substituted ligands as [B 4 ] and [C 4 ], respectively, following our previous notation [23,24] (Figure 2).Assuming that the axial anisotropy of the magnetic susceptibility tensor χ ax deduced from NMR correlates with ligand field parameters [30][31][32][33][34], which, in turn, determine SMM behavior [26,35], these studies serve to verify the theoretical models and to design molecular magnets with improved characteristics [36].

Synthesis and Characterization
The synthesis of the [B 4 ]M[B 4 ]M*[C 4 ] complexes with the aforementioned combinations of metal centers followed the procedure reported for Y 3+ and Tb 3+ complexes (Figure 3).It started with the preparation of bis(octa-n-butoxyphthalocyaninates) M[B 4 ] 2 , M = Tb 3+ , Dy 3+ and Y 3+ by the interaction of H 2 [B 4 ] with the corresponding metal acetates in a refluxing mixture of 1-chloronaphthalene (CN) and 1-octanol (OctOH) and DBU.Previously, we demonstrated that the use of this mixture is beneficial for the fast and efficient conversion of the starting ligand into metal complexes [23].

Synthesis and Characterization
The synthesis of the [B4]M[B4]M*[C4] complexes with the aforementioned combinations of metal centers followed the procedure reported for Y 3+ and Tb 3+ complexes (Figure 3).It started with the preparation of bis(octa-n-butoxyphthalocyaninates) M[B4]2, M = Tb 3+ , Dy 3+ and Y 3+ by the interaction of H2[B4] with the corresponding metal acetates in a refluxing mixture of 1-chloronaphthalene (CN) and 1-octanol (OctOH) and DBU.Previously, we demonstrated that the use of this mixture is beneficial for the fast and efficient conversion of the starting ligand into metal complexes [23].
The resulting complexes underwent a reaction with tetra-15-crown-5-phthalocyanine, H2[C4], in the presence of acetylacetonates M*(acac)3•nH2O in a refluxing mixture of 1,2,4-trichlorobenzene (TCB) and 1-octanol.Rapid conversion of the starting reagents into the desired triple-deckers was observed using UV-vis spectroscopy, and the target com  UV-vis spectra measured in toluene and chloroform demonstrate that dysprosiumcontaining complexes exhibit pronounced solvatochromism, similarly to the previously reported Tb-and Y-containing [B4]M[B4]M*[C4] complexes [23,24] (Figure 4).In particular, the spectra of the newly synthesized complexes in toluene contain intense Q-bands with well-resolved splitting (642-643 and 695-698 nm) and less intense Soret and N-bands at 363-364 and 292-294 nm, while in the spectra in chloroform, the Q-bands are severely broadened, although the maxima of the Q-bands in the two different solvents are very UV-vis spectra measured in toluene and chloroform demonstrate that dysprosiumcontaining complexes exhibit pronounced solvatochromism, similarly to the previously reported Tb-and Y-containing [23,24] (Figure 4).In particular, the spectra of the newly synthesized complexes in toluene contain intense Q-bands with well-resolved splitting (642-643 and 695-698 nm) and less intense Soret and N-bands at 363-364 and 292-294 nm, while in the spectra in chloroform, the Q-bands are severely broadened, although the maxima of the Q-bands in the two different solvents are very close.This has been previously explained by the difference in the conformational state of molecules, as in the spectrum of the conformer with DP/SAP coordination polyhedral, more electronic transitions are symmetrically allowed in comparison with the SAP/SAP conformer [25,37].Since the positions and splitting of Q-bands in the UV-vis spectra of sandwich complexes are governed by the interligand distances, which in turn depend on the size of the metal centers [19,37,38], the faintly small difference in the spectral appearances of the synthesized complexes in each of the solvents is explained by the similarity of the radii of Y 3+ , Tb 3+ and Dy 3+ ions [39].
close.This has been previously explained by the difference in the conformational state of molecules, as in the spectrum of the conformer with DP/SAP coordination polyhedral, more electronic transitions are symmetrically allowed in comparison with the SAP/SAP conformer [25,37].Since the positions and splitting of Q-bands in the UV-vis spectra of sandwich complexes are governed by the interligand distances, which in turn depend on the size of the metal centers [19,37,38], the faintly small difference in the spectral appearances of the synthesized complexes in each of the solvents is explained by the similarity of the radii of Y 3+ , Tb 3+ and Dy 3+ ions [39].The very weak absorption of lanthanide ions in the synthesized complexes could not be detected due to the much stronger absorption of tetrapyrrolic ligands.Also, we could not detect f-luminescence of Tb 3+ and Dy 3+ ions due to reabsorption of weak lanthanidecentered emission by tetrapyrrolic ligands, which have absorption bands with extinction coefficients of 1.4−2.0× 10 4 L mol −1 cm −1 in the range of 400-600 nm, where the emission of these ions is typically observed in coordination compounds with colorless or weakly absorbing ligands [40].Nevertheless, there are reports of f-luminescence observed in the near-IR for Nd(III), Ho(III), Er(III) and Yb(III) and complexes where Pcs and related macrocycles typically do not absorb light [41].
While optical methods are not suitable to confirm the presence of lanthanides in complexes with phthalocyanine ligands, the chemical composition of the synthesized tripledeckers is unambiguously confirmed by MALDI-TOF mass spectrometry due to the good agreement between the calculated and experimentally observed isotopic distributions (Figures S1-S5).The exact arrangement of the ligands and metal centers is determined by NMR spectroscopy, as discussed in the following sections.

Analysis of Lanthanide-Induced Shifts in 1 H-NMR Spectra
The 1 H-NMR spectra of lanthanide complexes typically have strikingly different appearances from the spectra of diamagnetic organic compounds since the presence of paramagnetic metal centers causes up-or downfield shifts of resonance signals by tens or even hundreds of ppm.The signs and magnitudes of these lanthanide-induced shifts (LISs) depend on both the nature of the lanthanide ion [42][43][44][45] and the overall geometry of the complex [46,47].The very weak absorption of lanthanide ions in the synthesized complexes could not be detected due to the much stronger absorption of tetrapyrrolic ligands.Also, we could not detect f -luminescence of Tb 3+ and Dy 3+ ions due to reabsorption of weak lanthanidecentered emission by tetrapyrrolic ligands, which have absorption bands with extinction coefficients of 1.4−2.0× 10 4 L mol −1 cm −1 in the range of 400-600 nm, where the emission of these ions is typically observed in coordination compounds with colorless or weakly absorbing ligands [40].Nevertheless, there are reports of f -luminescence observed in the near-IR for Nd(III), Ho(III), Er(III) and Yb(III) and complexes where Pcs and related macrocycles typically do not absorb light [41].
While optical methods are not suitable to confirm the presence of lanthanides in complexes with phthalocyanine ligands, the chemical composition of the synthesized triple-deckers is unambiguously confirmed by MALDI-TOF mass spectrometry due to the good agreement between the calculated and experimentally observed isotopic distributions (Figures S1-S5).The exact arrangement of the ligands and metal centers is determined by NMR spectroscopy, as discussed in the following sections.

Analysis of Lanthanide-Induced Shifts in 1 H-NMR Spectra
The 1 H-NMR spectra of lanthanide complexes typically have strikingly different appearances from the spectra of diamagnetic organic compounds since the presence of paramagnetic metal centers causes up-or downfield shifts of resonance signals by tens or even hundreds of ppm.The signs and magnitudes of these lanthanide-induced shifts (LISs) depend on both the nature of the lanthanide ion [42][43][44][45] and the overall geometry of the complex [46,47].
Thus, for the k-th proton, the LIS can be expressed as the difference between its chemical shifts in isostructural para-and diamagnetic complexes (1): On the other hand, the LIS can be presented as the sum of typically negligible contact (through-bond, δ con k ) and predominant dipolar (through-space, δ dip k ) contributions: Here, G k is a geometrical parameter depending on the distance r k between the k-th proton and the lanthanide ion, and θ k is the angle between the vector (H k ; Ln 3+ ) and the main symmetry axis, D 4 in our case (Figure 5a).The spectral-structural correlation (2) renders lanthanide ions as perfect probes for the elucidation of the solution structures, for example, in structural biology [48,49].
Molecules 2024, 29, 510 6 of 18 Thus, for the k-th proton, the LIS can be expressed as the difference between its chemical shifts in isostructural para-and diamagnetic complexes (1): On the other hand, the LIS can be presented as the sum of typically negligible contact (through-bond,  ) and predominant dipolar (through-space,  ) contributions: Here,  is a geometrical parameter depending on the distance  between the k-th proton and the lanthanide ion, and  is the angle between the vector (H ; Ln ) and the main symmetry axis, D4 in our case (Figure 5a).The spectral-structural correlation (2) renders lanthanide ions as perfect probes for the elucidation of the solution structures, for example, in structural biology [48,49].4) vs. the experimental values (axes x).Grey lines show the least-squares fits between the calculated and experimental values.The complete assignment of NMR spectra for all complexes is presented in the Materials and Methods section and the Electronic Supporting Information (Figures S7-S14).
The proportionality factor  is the axial component of the magnetic susceptibility tensor.Thus, the observation of lanthanide-induced shifts in NMR spectra is a manifestation of lanthanide magnetic properties associated with crystal field parameters [42][43][44]46,47], and the value of  complements the easily affordable NMR with much more sophisticated time-and resource-consuming magnetochemical measurements [10,31,[50][51][52][53].
The dipolar approximation of the LIS using Equation (2) suggests that the ratio of LISs for the pair of k-th and l-th protons can be approximated with the ratio of their geometrical parameters,  : A combination of Equations ( 1) and ( 3) written for a pair of protons gives Equation ( 4), which can be used to estimate the positions of the resonance signals in the 1 H-NMR the experimental values (axes x).Grey lines show the least-squares fits between the calculated and experimental values.The complete assignment of NMR spectra for all complexes is presented in the Materials and Methods section and the Electronic Supporting Information (Figures S7-S14).
The proportionality factor χ ax is the axial component of the magnetic susceptibility tensor.Thus, the observation of lanthanide-induced shifts in NMR spectra is a manifestation of lanthanide magnetic properties associated with crystal field parameters [42][43][44]46,47], and the value of χ ax complements the easily affordable NMR with much more sophisticated time-and resource-consuming magnetochemical measurements [10,31,[50][51][52][53].
The dipolar approximation of the LIS using Equation (2) suggests that the ratio of LISs for the pair of k-th and l-th protons can be approximated with the ratio of their geometrical parameters, R kl : A combination of Equations ( 1) and ( 3) written for a pair of protons gives Equation ( 4), which can be used to estimate the positions of the resonance signals in the 1 H-NMR spectra of paramagnetic complexes using the spectrum of a diamagnetic complex with a well-established structure and the shift of at least one firmly assigned resonance signal in the spectrum of a paramagnetic complex: Yttrium(III) complexes are typically used as diamagnetic references, providing δ dia k values, and the ratio of geometrical parameters R kl can be taken from the structural models obtained from either X-ray characterization or DFT modeling [22].
The almost perfect agreement of the UV-vis spectra of the newly synthesized complexes with those reported previously (Figure 4) suggests an analogy in the solution structure of the entire series of which is generally expected given the close values of the ionic radii of Y 3+ (1.019 pm), Tb 3+ (1.040 pm) and Dy 3+ (1.027 pm) [39], together with the similarity of their coordination chemistry [4].This conclusion justifies the following application of the previously determined X-ray structures of solvated either with dichloromethane (DP/SAP conformer) or with toluene (SAP/SAP conformer) [23] for the analysis of 1 H-NMR spectra of paramagnetic complexes containing Tb 3+ and Dy 3+ ions using Equation (4).
It is worth noting that the X-ray structures of the solvates with CHCl Only aromatic protons and methylene protons proximal to the Pc ligands-1 i,o -CH 2 groups of butoxy-substituents and α o -CH 2 groups of crown-ether rings (black labels in Figure 5b)-were taken for further analysis as they refer to the most rigid part of the molecules, assuring the consistency between the X-ray and solution structures [23,24].Although CH 2 groups can rotate around single σ-bonds, the interconversion between the exoand endo-H protons of the considered methylenes is slow on the NMR timescale; their resonances appear as pairs of signals separated by ca. 10 ppm and coupled in 1 H-1 H COSY spectra.Such hindered rotation can be explained by numerous weak interactions of the aliphatic substituents with the substituents in the neighboring ligands [23] and solvent molecules [25].Similar splitting patterns were reported previously by Enders and Yamashita et al. [26,54], Jiang et al. [55] and Ishikawa et al. [56], together with our publications [19,23,24].
The splitting is much smaller for more distant methylene groups (grey labels in Figure 5b), and it almost vanishes for the 3 i,o -CH 2 and CH 3 groups, suggesting that these groups rotate nearly freely around the corresponding sigma bonds.Thus, these distant protons were not taken into consideration so that the deduced χ ax value is not disturbed by the conformational flexibility of the peripheral substituents.Good correlations between the experimentally observed chemical shifts and the shifts calculated from X-ray data (Figure 5c-f) validate the applied assumptions and approximations.
The multidirectional nature of LISs observed in Figure 5c-f in such complexes originates from the fact that different protons have either positive or negative G k depending on the location of these protons in the principle magnetic framework of the Dy 3+ ion (Figure 6b,c), and the largest absolute value of G k deduced from X-ray data corresponds to the aromatic proton of the inner BuO-substituted ligand bH Pc i [24].This observation allowed us to assign the resonance signals of the other protons using Equation ( 4), which was confirmed by further 1 H-1 H COSY spectra (Figures S7-S14).
Molecules 2024, 29, 510 8 of 18 The multidirectional nature of LISs observed in Figure 5c-f in such complexes originates from the fact that different protons have either positive or negative  depending on the location of these protons in the principle magnetic framework of the Dy 3+ ion (Figure 6b,c), and the largest absolute value of  deduced from X-ray data corresponds to the aromatic proton of the inner BuO-substituted ligand bHPc i [24].This observation allowed us to assign the resonance signals of the other protons using Equation ( 4), which was confirmed by further 1 H-1 H COSY spectra (Figures S7-S14).Plotting LIS from NMR spectra vs. the geometrical parameters found from X-ray data, followed by linear regression analysis, afforded axial anisotropies of Dy 3+ ions in different coordination surroundings (Figure 6a).This analysis demonstrates that the magnetic properties of this ion are sensitive to the symmetry of the coordination surrounding.Thus, switching from an SAP to a DP polyhedron in toluene and chloroform, respectively, causes a notable increase in axial anisotropy, similar to the results reported previously for Tb 3+ ions in isostructural complexes.

Analysis of Lanthanide-Induced Shifts in 1 H-NMR Spectra of the Homonuclear Complex [B4]Dy[B4]Dy[C4]
Due to the presence of two structurally nonequivalent paramagnetic centers in [B4]Dy[B4]Dy[C4], Equation (4) cannot be applied for the assignment of its 1 H-NMR spectra, as these centers might give different contributions to the total LIS value.Therefore, in the present case, the assignment is made by approximating the chemical shift of the k-th proton in the homonuclear complex with the sum of the shifts of analogous protons in the heteronuclear complexes according to Equation ( 5) [56].
The assignment is confirmed by 1 H-1 H COSY. In addition, due to the proximity of all the considered protons to the Dy 3+ ions, their signals are the most broadened, providing additional verification of their correct assignment.
A comparison of the 1 H-NMR spectra of the homonuclear complex [B4]Dy[B4]Dy[C4] measured in CDCl3 and toluene-d8 (Figure 7a,b) suggests that switching from the SAP/SAP conformer to the DP/SAP conformer also causes an increase in anisotropy, as evidenced by the broadening of the spectral range.To confirm this conclusion, we assumed that each of the Dy 3+ ions has its own  value in accordance with Equation ( 6); thus, we used two-dimensional minimization of the Wilcott factor (AF, Equation ( 7)) Plotting LIS from NMR spectra vs. the geometrical parameters found from X-ray data, followed by linear regression analysis, afforded axial anisotropies of Dy 3+ ions in different coordination surroundings (Figure 6a).This analysis demonstrates that the magnetic properties of this ion are sensitive to the symmetry of the coordination surrounding.Thus, switching from an SAP to a DP polyhedron in toluene and chloroform, respectively, causes a notable increase in axial anisotropy, similar to the results reported previously for Tb 3+ ions in isostructural complexes.4) cannot be applied for the assignment of its 1 H-NMR spectra, as these centers might give different contributions to the total LIS value.Therefore, in the present case, the assignment is made by approximating the chemical shift of the k-th proton in the homonuclear complex with the sum of the shifts of analogous protons in the heteronuclear complexes according to Equation ( 5) [56].

Analysis of Lanthanide-Induced
The assignment is confirmed by 1 H-1 H COSY. In addition, due to the proximity of all the considered protons to the Dy 3+ ions, their signals are the most broadened, providing additional verification of their correct assignment.
A comparison of the 1 H-NMR spectra of the homonuclear complex [B 4 ]Dy[B 4 ]Dy[C 4 ] measured in CDCl 3 and toluene-d 8 (Figure 7a,b) suggests that switching from the SAP/SAP conformer to the DP/SAP conformer also causes an increase in anisotropy, as evidenced by the broadening of the spectral range.To confirm this conclusion, we assumed that each of the Dy 3+ ions has its own χ ax value in accordance with Equation ( 6); thus, we used two-dimensional minimization of the Wilcott factor (AF, Equation ( 7)) across the various pairs of χ ax to find the best agreement between the calculated and experimental LIS values (Figure 7c,d).
across the various pairs of  to find the best agreement between the calculated and experimental LIS values (Figure 7c,d).
The results presented in Figure 7

complexes, [B4]Dy[B4]Tb[C4] and [B4]Tb[B4]Dy[C4]
, in CDCl3 and toluene-d8 (Figure 9) to trace the change in the axial anisotropy of both lanthanide ions simultaneously using Equations ( 5) and ( 6).The assignment was confirmed by 1 H-1 H COSY (Figures S13-S16).9) to trace the change in the axial anisotropy of both lanthanide ions simultaneously using Equations ( 5) and ( 6).The assignment was confirmed by 1 H-1 H COSY (Figures S13-S16).9) to trace the change in the axial anisotropy of both lanthanide ions simultaneously using Equations ( 5) and ( 6).The assignment was confirmed by 1 H-1 H COSY (Figures S13-S16).9a,b), which is expectedly followed by a significant decrease in its axial anisotropy-from 9.40 × 10 −31 to 7.75 × 10 −31 m 3 -in agreement with the behavior of the previously studied Tb 3+ complexes [23,24].Switching the coordination polyhedron of the Dy 3+ ion from DP to SAP also causes a decrease in axial anisotropy (Figure 9c,d), although this effect is not so strong compared to the Tb 3+ metal centers (Figure 10a,b).bels of protons are given in Figure 5b.Graphical search for values of  for corresponding lanthanide cations at the [B4]/[C4] (axes x) and [B4]/[B4] (axes y) sites corresponding to minimal values of the agreement factor, AF, in CDCl3 (e,g) and toluene-d8 (f,h).The complete assignment of NMR spectra for all complexes is presented in the Materials and Methods section and the Electronic Supporting Information (Figures S19-S26).
It can be clearly seen that the change from CDCl3 to toluene-d8 has the most pronounced influence on the 1

H-NMR spectra of [B4]Tb[B4]Dy[C4],
where the Tb 3+ ion is placed into the switchable [B4]/[B4] site (Figure 9a,b), which is expectedly followed by a significant decrease in its axial anisotropy-from 9.40 × 10 −31 to 7.75 × 10 −31 m 3 -in agreement with the behavior of the previously studied Tb 3+ complexes [23,24].Switching the coordination polyhedron of the Dy 3+ ion from DP to SAP also causes a decrease in axial anisotropy (Figure 9c,d), although this effect is not so strong compared to the Tb 3+ metal centers (Figure 10a,b).
Finally, the availability of  for each metal center allows plotting the contour maps of the overall LIS values for different conformers of  10c-f).These plots vividly explain the difference between the spectral appearances of these complexes and their homonuclear analogues.For example, it can be seen that a combination of two metal ions with essentially different anisotropies results in deformation of the zero Δ isosurface, which approaches the protons of the terminal ligands bound to the Dy 3+ ion; this is clearly seen in the experimental spectra, in which the less downfield-shifted signals refer to aromatic protons (either bHPc o or сHPc o ) and exo-protons of methylene groups (1 o ′ or α o ′, respectively).(Figure 10c-f).These plots vividly explain the difference between the spectral appearances of these complexes and their homonuclear analogues.For example, it can be seen that a combination of two metal ions with essentially different anisotropies results in deformation of the zero ∆δ isosurface, which approaches the protons of the terminal ligands bound to the Dy 3+ ion; this is clearly seen in the experimental spectra, in which the less downfield-shifted signals refer to aromatic protons (either bH Pc o or cH Pc o ) and exo-protons of methylene groups (1 o′ or α o′ , respectively).

Discussion
The presented results give several guidelines for further work on lanthanide complexes with phthalocyanine ligands:

•
The trisphthalocyanine scaffold affords the synthesis of heteronuclear complexes with a precise arrangement of rare-earth ions due to its thermodynamical and kinetic stability.
Complexes with different combinations of paramagnetic lanthanides can be efficiently obtained, which provides the basis for further investigation of intramolecular f -f interactions and the elaboration of molecular magnetic materials.In the present work, the combinations of two strongly paramagnetic Tb 3+ and Dy 3+ ions were used to obtain isomeric heteronuclear complexes, but, obviously, other combinations of middle and late lanthanides can also be used in this type of chemistry to obtain complexes with the required number of unpaired f -electrons.• The comprehensive 1 H-NMR spectroscopic characterization of strongly paramagnetic complexes is based on the appropriate structural model; therefore, this work provides algorithms for dealing with the spectra of complexes containing either one or two lanthanide ions, which are not necessarily equivalent.In this regard, the results presented demonstrate that the application of paramagnetic 1 H-NMR spectroscopy should not be limited to routine identification but can be used to extract the magnetic properties of lanthanide ions [53].In this context, our report follows the strategies applied by Enders and Yamashita, where Tb(III) and Dy(III) sandwich phthalocyaninates were comprehensively studied using NMR spectroscopy [26,35,54,57,58], and the influence of electronic and structural effects on their magnetic properties, especially χ ax , was revealed.In summary, it is expected that further magnetochemical studies of the newly synthesized lanthanide phthalocyaninates will provide more correlations between the χ ax term and the energetic properties of slow magnetic relaxation.• The addition of controllable conformational flexibility gives one more degree of freedom to control the magnetic properties of sandwich lanthanide complexes.The previously discovered correlations between the symmetry of the coordination polyhedron of the Tb 3+ ion and its magnetic properties are also valid for the Dy 3+ ion, and the effect of the conformational switching can be studied for other lanthanides to find the capabilities and limitations of theoretical models.

•
Importantly, our study evidences that the χ ax values for Dy 3+ are nearly twice smaller than those for Tb 3+ ; however, this observation does not match the expectations from Bleaney's theory, where the largest anisotropy in the lanthanide series is expected for dysprosium [43].Moreover, it contradicts the results of theoretical modeling obtained by Mironov et al. [32] for various polyhedra of lanthanide complexes, where the most pronounced influence of the surrounding coordination was anticipated for dysprosium complexes.The reason for this discrepancy may be a violation of the theory's basic assumption that the thermal energy is larger than the ligand field splitting; thus, further theoretical modeling using ab initio methods might be particularly helpful [59,60].were synthesized according to the previously reported procedures [61,62].1,2,4-Trichlorobenzene (TClB, for synthesis, 1-octanol (for synthesis), rare-earth acetylacetonates (Sigma-Aldrich, Burlington, MA, USA), and neutral alumina (50-200 µm, Macherey-Nagel, Düren, Germany) were used as received from the commercial suppliers.Chloroform (reagent grade, Ekos-1, Staraya Kupavna, Russia) was distilled over CaH 2 .

Methods
Matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectra were measured on a Bruker Daltonics Ultraflex spectrometer.Mass spectra were registered in positive ion mode using 2,5-dihydroxybenzoic acid as a matrix.UV-vis spectra in the range of 250-900 nm were measured using a Thermo Evolution 210 spectrophotometer in quartz cells with 0.5-1.0cm optical path lengths.
1 H NMR spectra were recorded at 303 K on a Bruker Avance III 600 MHz spectrometer equipped with a 5 mm Z-gradient BBO probe (zg30 pulse program from Topspin library).A total of 128 scans with a 30 • excitation pulse (90 • pulse width was 13.8 µs, DE 6.50 µs) and 1 s delay were accumulated for each system.Spectral width was taken to be 90-100 kHz, depending on the combination of lanthanides.Gaussian multiplication was used for processing.Automatic polynomial baseline correction was applied.For recording 2D COSY spectra (cosygpqf program from Topspin library), spectral width was selected to be 40-75 ppm, depending on the combination of lanthanides.The residual solvent resonances were used as internal references (δ toluene 7.09 ppm, chloroform 7.26 ppm).Typically, 5 mg of complexes was dissolved in 0.6 mL of a corresponding deuterated solvent to provide a concentration of ca.2.3 mM.The applied deuterated chloroform (99.8 atom% D, ZEOchem, Uetikon am See, Switzerland) was filtered prior to use through Pasteur pipettes filled with alumina to remove possible acidic impurities.Deuterated toluene (99.8 atom% D, ABCR, Karlsruhe, Germany) was used without additional purification.
) was added.After 10 min, the reaction mixture was cooled to ambient temperature.The reaction mixture was transferred to a chromatographic column packed with alumina in a mixture of CHCl 3 and hexane (1:1 v/v).The target complex was isolated by elution with a CHCl 3 -hexane mixture (4:1 v/v), followed by a mixture of CHCl 3 with 0→1% MeOH as a dark-blue solid ( (27.0 mg, 11.5 µmol) and H 2 [C 4 ] (18.0 mg, 14.1 µmol) in a mixture of 2.7 mL 1,2,4-trichlorobenzene and 0.3 mL n-octanol was heated to reflux at 230 • C under a stream of argon, and solid Dy(acac) 3 •nH 2 O (21.0 mg, 45.7 µmol) was added.After 8 min, the reaction mixture was cooled to ambient temperature.The reaction mixture was transferred to a chromatographic column packed with alumina in a mixture of CHCl 3 and hexane (1:1 v/v).The target complex was isolated by elution with a CHCl 3 -hexane mixture (4:1 v/v), followed by a mixture of CHCl 3 with 0→1% MeOH as a dark-blue solid (20   The target complex was isolated by elution with a CHCl 3 -hexane mixture (4:1 v/v), followed by a mixture of CHCl 3 with 0→1% MeOH as a dark-blue solid (
plexes were isolated by column chromatography on alumina.Due to the difference in polarity, [B4]M[B4]M*[C4] could be easily separated from the unreacted M[B4]2 and the side homoleptic products M*2[C4]3.

Figure 3 .
Figure 3. Synthesis of heteronuclear trisphthalocyaninates [B 4 ]M[B 4 ]M*[C 4 ].The resulting complexes underwent a reaction with tetra-15-crown-5-phthalocyanine, H 2 [C 4 ], in the presence of acetylacetonates M*(acac) 3 •nH 2 O in a refluxing mixture of 1,2,4-trichlorobenzene (TCB) and 1-octanol.Rapid conversion of the starting reagents into the desired triple-deckers was observed using UV-vis spectroscopy, and the target complexes were isolated by column chromatography on alumina.Due to the difference in polarity, [B 4 ]M[B 4 ]M*[C 4 ] could be easily separated from the unreacted M[B 4 ] 2 and the side homoleptic products M* 2 [C 4 ] 3 .UV-vis spectra measured in toluene and chloroform demonstrate that dysprosiumcontaining complexes exhibit pronounced solvatochromism, similarly to the previously reported Tb-and Y-containing[B 4 ]M[B 4 ]M*[C 4 ] complexes[23,24] (Figure4).In particular, the spectra of the newly synthesized complexes in toluene contain intense Q-bands with well-resolved splitting (642-643 and 695-698 nm) and less intense Soret and N-bands at 363-364 and 292-294 nm, while in the spectra in chloroform, the Q-bands are severely broadened, although the maxima of the Q-bands in the two different solvents are very close.This has been previously explained by the difference in the conformational state of molecules, as in the spectrum of the conformer with DP/SAP coordination polyhedral, more electronic transitions are symmetrically allowed in comparison with the SAP/SAP conformer[25,37].Since the positions and splitting of Q-bands in the UV-vis spectra of sandwich complexes are governed by the interligand distances, which in turn depend on the size of the metal centers[19,37,38], the faintly small difference in the spectral appearances

Figure 5 .
Figure 5. Explanation of values of  and  included in the geometric parameter  according to Equation (2) (a).Labels of protons: black labels indicate protons used for further analysis (b). 1 H-NMR spectra of [B4]Dy[B4]Y[C4] (c,d) and [B4]Y[B4]Dy[C4] (e,f) in CDCl3 (upper row) and toluene-d8 (lower row).Dots show the positions of chemical shifts calculated with Equation (4) vs. the experimental values (axes x).Grey lines show the least-squares fits between the calculated and experimental values.The complete assignment of NMR spectra for all complexes is presented in the Materials and Methods section and the Electronic Supporting Information (FiguresS7-S14).

Figure 5 .
Figure 5. Explanation of values of r k and θ k included in the geometric parameter G k according to Equation (2) (a).Labels of protons: black labels indicate protons used for further analysis (b). 1 H-NMR spectra of [B 4 ]Dy[B 4 ]Y[C 4 ] (c,d) and [B 4 ]Y[B 4 ]Dy[C 4 ] (e,f) in CDCl 3 (upper row) and toluene-d 8 (lower row).Dots show the positions of chemical shifts calculated with Equation (4) vs.the experimental values (axes x).Grey lines show the least-squares fits between the calculated and experimental values.The complete assignment of NMR spectra for all complexes is presented in the Materials and Methods section and the Electronic Supporting Information (FiguresS7-S14).

4 ]
3 are not yet available, although we attempted to grow single crystals by diffusion of acetonitrile, heptane or vapors of pentane into[B 4 ]Y[B 4 ]Y[C 4 ]solutions in CHCl 3 but failed to obtain material of sufficient quality to perform XRD experiments with the required precision.On the other hand, we previously demonstrated that chloroform also stabilizes staggered and gauche conformers of 15C5-and BuO-substituted triple-deckers, similarly to dichloromethane (see FigureS6for more details).Moreover, we previously showed that the spectra of[B 4 ]Tb[B 4 ]Tb[C 4 ]in CDCl 3 and CD 2 Cl 2 are almost indistinguishable[23].Altogether, this suggests that the X-ray structure of the DP/SAP conformer solvated with CH 2 Cl 2 is a reasonable structural model for the analysis of LISs measured in chloroform.2.2.1.Analysis of Lanthanide-InducedShifts in 1 H-NMR Spectra of Heteronuclear Complexes [B 4 ]Dy[B 4 ]Y[C 4 ] and [B 4 ]Y[B 4 ]Dy[C The 1 H-NMR spectra of the pair of isomeric trisphthalocyaninates [B 4 ]Dy[B 4 ]Y[C 4 ] and [B 4 ]Y[B 4 ]Dy[C 4 ] measured in CDCl 3 and toluene-d 8 evidence that the spectral range (SR) from the most upfield-shifted to the most downfield-shifted signals depends on both the position of the Dy 3+ ion and the solvent applied (Figure 5b,c).

Figure 6 .
Figure 6.Values of axial anisotropy  found by least-squares linearization of Δ vs.  according to Equation (2) for Dy 3+ ions in different surroundings (a).Contour maps of the (, ) function plotted for [B4]Dy[B4]Y[C4] (b) and [B4]Y[B4]Dy[C4] (c); the coordinates of selected protons were taken from the X-ray structures of the corresponding solvates with dichloromethane (cyan balls) and toluene (dark-blue balls).Labels of protons are given in Figure 5b.

Figure 6 .
Figure 6.Values of axial anisotropy χ ax found by least-squares linearization of ∆δ k vs. G k according to Equation (2) for Dy 3+ ions in different surroundings (a).Contour maps of the G(θ, r) function plotted for [B 4 ]Dy[B 4 ]Y[C 4 ] (b) and [B 4 ]Y[B 4 ]Dy[C 4 ] (c); the coordinates of selected protons were taken from the X-ray structures of the corresponding solvates with dichloromethane (cyan balls) and toluene (dark-blue balls).Labels of protons are given in Figure 5b.

Figure 7 .
Figure 7. 1 H-NMR spectra of [B4]Dy[B4]Dy[C4] in CDCl3 (a) and toluene-d8 (b).Dots show the chemical shifts calculated with Equation (5) (axes y) vs. the experimental values (axes x).Grey lines show the least-squares fits between the calculated and experimental values.Labels of protons are given in Figure 5c.Graphical search for values of  for Dy 3+ cations at [B4]/[C4] (axes x) and [B4]/[B4] (axes y) sites corresponding to minimal values of the agreement factor, AF, in CDCl3 (c) and toluene-d8 (d).Labels of protons are given in Figure 5b.The complete assignment of NMR spectra for all complexes is presented in the Materials and Methods section and the Electronic Supporting Information (Figures S15-S18).Plotting the contour maps of Δ for the molecules of the [B4]Dy[B4]Dy[C4] complex in CDCl3 and toluene-d8 using the derived  values explains why all signals in the spectra have shifts of the same sign-all protons fall into the region with the same sign of the net (, ) function, in contrast to the heteronuclear derivatives, where protons fall into regions with different signs of the geometrical parameter (Figure8b,c).

Figure 7 .
Figure 7. 1 H-NMR spectra of [B 4 ]Dy[B 4 ]Dy[C 4 ] in CDCl 3 (a) and toluene-d 8 (b).Dots show the chemical shifts calculated with Equation (5) (axes y) vs. the experimental values (axes x).Grey lines show the least-squares fits between the calculated and experimental values.Labels of protons are given in Figure 5c.Graphical search for values of χ ax for Dy 3+ cations at [B 4 ]/[C 4 ] (axes x) and [B 4 ]/[B 4 ] (axes y) sites corresponding to minimal values of the agreement factor, AF, in CDCl 3 (c) and toluene-d 8 (d).Labels of protons are given in Figure 5b.The complete assignment of NMR spectra for all complexes is presented in the Materials and Methods section and the Electronic Supporting Information (Figures S15-S18).The results presented in Figure 7 suggest that the change in anisotropy in the homonuclear complex [B 4 ]Dy[B 4 ]Dy[C 4 ] follows the same trends observed in heteronuclear complexes.In particular, switching between SAP/SAP (in toluene-d 8 ) and DP/SAP (in CDCl 3 ) causes an increase in the anisotropy of the ion in the conformationally flexible [B 4 ]/[B 4 ] site from 3.95 × 10 −31 to 4.28 × 10 −31 m 3 , while the χ ax value of the Dy 3+ ion in the conformationally invariant [B 4 ]/[C 4 ] site remains almost the same-3.9× 10 −31 m 3 (Figure 8a).Plotting the contour maps of ∆δ para k for the molecules of the [B 4 ]Dy[B 4 ]Dy[C 4 ] complex in CDCl 3 and toluene-d 8 using the derived χ ax values explains why all signals in the spectra have shifts of the same sign-all protons fall into the region with the same sign of the net G(θ, r) function, in contrast to the heteronuclear derivatives, where protons fall into regions with different signs of the geometrical parameter (Figure 8b,c).

Figure 10 . 10 .
Figure 10.Values of axial anisotropy  for Tb 3+ and Dy 3+ ions in heteronuclear complexes [B4]Tb[B4]Dy[C4] (a) and [B4]Dy[B4]Tb[C4] (b).Contour maps of the net LIS functions according to Figure 10.Values of axial anisotropy χ ax for Tb 3+ and Dy 3+ ions in heteronuclear complexes [B 4 ]Tb[B 4 ]Dy[C 4 ] (a) and [B 4 ]Dy[B 4 ]Tb[C 4 ] (b).Contour maps of the net LIS functions according to Equation (6) plotted for [B 4 ]Tb[B 4 ]Dy[C 4 ] (c,d) and [B 4 ]Dy[B 4 ]Tb[C 4 ] (e,f) in CDCl 3 (a,c) and toluene-d 8 (b,d).Labels of protons are given in Figure 5b.Finally, the availability of χ ax for each metal center allows plotting the contour maps of the overall LIS values for different conformers of [B 4 ]Tb[B 4 ]Dy[C 4 ] and [B 4 ]Dy[B 4 ]Tb[C 4 ](Figure10c-f).These plots vividly explain the difference between the spectral appearances of these complexes and their homonuclear analogues.For example, it can be seen that a combination of two metal ions with essentially different anisotropies results in deformation