Tetrathiafulvalene-supported triple-decker phthalocyaninato dysprosium(III) complex: synthesis, properties and surface assembly.

Self-assembly of functional compounds into a prerequisite nanostructure with desirable dimension and morphology by controlling and optimizing intermolecular interaction attracts an extensive research interest for chemists and material scientist. In this work, a new triple-decker sandwich-type lanthanide complex with phthalocyanine and redox-active Schiff base ligand including tetrathiafulvalene (TTF) units has been synthesized, and characterized by single crystal X-ray diffraction analysis, absorption spectra, electrochemical and magnetic measurements. Interestingly, the non-centrosymmetric target complex displays a bias dependent selective adsorption on a solid surface, as observed by scanning tunneling microscopy (STM) at the single molecule level. Density function theory (DFT) calculations are utilized to reveal the formation mechanism of the molecular assemblies, and show that such electrical field dependent selective adsorption is regulated by the interaction between the external electric field and intrinsic molecular properties. Our results suggest that this type of multi-decker complex involving TTF units shows intriguing multifunctional properties from the viewpoint of structure, electric and magnetic behaviors, and fabrication through self-assembly.

With the development of science and technology, the investigation and application of new functional materials have been the main task for production and life. Single-molecule magnets (SMMs) as molecular species exhibiting quantum tunneling and slow magnetic-relaxation processes, have been investigated from the viewpoint of applied science for use in next-generation devices, such as memory storage, quantum computers and spin-based molecular electronics12345. As a representative, sandwich-type phthalocyaninato and/or porphyrinato lanthanide complexes possess intriguing molecular structures, unique electronic and optical properties due to intramolecular interactions and intrinsic nature of the metal ion67891011, rendering them potentially useful in material science, including molecular spintronics, chemical sensors and field-effect transistors1213141516. Self-assembly of these functional molecules into a prerequisite nanostructure with desirable dimension and morphology via controlling and optimizing intermolecular interaction has been a subject of extensive research interest for chemists and material scientist17181920.

Recently, the need for new materials with more diversified and more sophisticated properties is continuously increasing. One of the goals is to prepare materials that possess not only one expected property or function but also combine two or more of them in one system. In particular, owing to the important applications in molecular spintronics, the synthesis of novel materials possessing synergy or interplay between electrical conductivity with magnetism has attracted attention in the last few years2122232425. Tetrathiafulvalene (TTF) and its derivatives, well-known sulfur-rich organic molecules, are ideally suited as components in these types of systems because they can act as reversible and stable electron donors. Their cationic radical salts constitute an important class of conducting molecular materials26272829, and have been incorporated as units on molecular machines3031. Absorbing redox-active TTF organic units on surfaces, and assembling SMMs systems comprised of these molecules with the paramagnetic metal ions may lead to new multifunctional materials with interesting structures and properties26323334353637383940.

Herein, we report the synthesis and surface assembly of the first sandwich-type triple-decker lanthanide complex based on a Schiff base ligand containing electrochemically active TTF unit H2L (H2L = 2,2′-((2-(4,5-bis(methylthio)-1,3-dithiol-2-ylidene)-1,3-benzodithiole-5,6-diyl)bis(nitrilomethylidyne)bis(4-chlorophenol)), which extends the classical system of sandwich-type tetrapyrrole oligomers limited to homoleptic or heterolepic phthalocyanine and/or porphyrin lanthanide multi-decker complexes67891011. For future applications in molecular electronics or spintronics, such as molecular field-effect transistors or spin field-effect transistors, the introduction of redox-active TTF derivatives will offer opportunities for new multifunctional molecular materials to fabricate single-molecule devices.

Results

Reaction of phthalocyanine, dysprosium acetylacetonate and TTF-supported Schiff base ligand in a 1:2:2 molar ratio in 1,2-dichloroethane and methanol resulted into the formation of target complex [Dy2(Pc)L2(CH3OH)]·5(ClCH2CH2Cl) (complex 1, Supplementary Fig. S1), which is stable in both solution and solid state, and corresponding characterization has been accomplished by single crystal X-ray diffraction analysis, absorption spectra, electrochemical and magnetic measurements. More interestingly, the assembling behavior of such novel complex was studied by scanning tunneling microscopy (STM) on a highly oriented pyrolytic graphite (HOPG) surface, for probing the surface patterns and revealing the interface information at the single-molecule level.

Molecule structural description

Complex 1 crystallizes in the triclinic space group Pī with the unique non-centrosymmetric tripler-decker sandwich-type molecular structure. As shown in Fig. 1, the Dy1 ion is eight-coordinated to six N atoms (four from the Pc and two from the H2L) and two phenolic O atoms, while the seven-coordinated Dy2 ion is surrounded by two N atoms of another TTF-supported Schiff base ligand and five O atoms (four from the H2L and one from a methanol molecule). Employing Continuous Shape Measures (CShMs) method41 allows us to determine the degree of distortion of the Dy(III) coordination sphere (Supplementary Table S1). The results show that the coordination geometry of Dy1 ion can be described as a similar square antiprism (SAPR) with S value of 0.928 (The S values indicate the proximity to the ideal polyhedron, S = 0 corresponds to the non-distorted polyhedron), while seven-coordinated Dy2 adopts a conformation close to a capped trigonal prism (CTPR). The two Dy3+ ions are bridged by the O1 and O2 atoms with Dy1···Dy2 distance of 3.887(5) Å. The dihedral angles between the mean planes of N1-N3-N5-N7 and O1-O2-N9-N10, O1-O2-N9-N10 and O3-O4-N11-N12, N1-N3-N5-N7 and O3-O4-N11-N12, are 1.10(15)°, 17.09(17)°, and 16.02(19)°, respectively. Two TTF cores adopt boat-like conformation, and the average deviation from a least-squares plane for one TTF core containing S1, S2, S3 and S4 atoms is 0.365 Å, while another TTF core containing S7, S8, S9 and S10 atoms is a deviation of 0.157 Å, indicating the distortion of TTF-supported Schiff base ligand in the middle layer is larger than that in the bottom layer. In the crystal packing (Supplementary Fig. S2), a wave-shaped conformation can be observed due to the different distortion of the TTF units, and the shortest intermolecular distance between the two Dy3+ ions is 8.541(6) Å without showing significant intermolecular interactions and shorter S···S contacts.

Figure 1

Molecular structure for complex 1 with thermal ellipsoids at 30% probability.

Hydrogen atoms and solvent molecules are omitted for clarity.

Electrochemical and spectroscopic properties

Electrochemical properties were investigated by cyclic voltammetry (Supplementary Fig. S3). The H2L ligand and complex 1 show similar two-step reversible oxidation processes, suggesting the presence of redox activities in present system. Additionally, UV-Vis absorption spectrum for complex 1 shows typical feature of phthalocyaninate metal complexes with a strong Q band around 686 nm and an additional weak absorption at about 615 nm (Supplementary Fig. S4), originating from the π-π* transitions of the phthalocyaninate ligand42.

Static and dynamic magnetic properties

The static magnetic behaviors for complex 1 were investigated under 100 Oe direct current (dc) magnetic field. As shown in Fig. 2a, the χ value of 27.55 cm3 K mol−1 at 300 K is slightly smaller than the theoretical value for two uncoupled Dy(III) ions (28.34 cm3 K mol−1, H, S = /, L = 5, J = /, g = /). When the temperature is lowered, the χ value drops to a minimum value of 23.52 cm3 K mol−1 at about 6.0 K, mainly ascribed to the progressive depopulation of excited Stark sublevels of Dy(III)1043. Upon further cooling, the curve increases abruptly again and reaches a maximum value of 24.58 cm3 K mol−1 at 1.8 K, likely due to the weak ferromagnetic Dy(III)-Dy(III) interactions44. Field-dependent magnetization for complex 1 increases rapidly at low fields (Fig. 2a, Inset), and then achieves a maximum value of 11.96 Nβ at 70 kOe without reaching the theoretical saturation value (20 Nβ for two isolated Dy(III) ions), indicating the presence of the crystal-field effects and magnetic anisotropy for Dy(III) ion in this sandwich system45.

Figure 2

(a) Temperature-dependent χMT values for complex 1 under 100 Oe dc field. Inset: Field-dependent magnetization plots for 1 at 1.8 K. (b) Temperature dependence of the in-phase and out-of-phase ac magnetic susceptibilities for 1 under 3 kOe dc field in the frequency range 1–1000 Hz.

Alternating current (ac) susceptibility measurements were carried out to study the dynamics of the magnetization for 1 (Supplementary Figs. S5 and S6). The temperature dependence of ac magnetic susceptibilities exhibit obvious frequency-dependent character in the in-phase (χ) and out-of-phase (χ) signals under H = 0 Oe and H = 2 Oe. Nevertheless, no maxima of χ signals were observed above 1.8 K owing to the fast relaxation associated with the quantum tunneling at zero magnetic field. Such behavior is attributed to the degeneracy of the ground Kramers states of each single Dy(III) ion, which is consistent with some similar Dy-SMMs1044. Performing an applied dc field of 3 kOe, the peaks of χ″ signals in the temperature-dependent ac susceptibilities clearly appear at frequencies as low as 16 Hz (Fig. 2b), indicating field induced SMM behavior with the suppression of the quantum tunneling of the magnetization (QTM) between sublevels. Furthermore, there are probably associated with the interesting dynamic magnetic behavior due to the slightly distinct anisotropic Dy(III) centers situating at the different coordination environment4446. According to frequency-dependent ac susceptibilities, two different χ″ signals can be observed as the temperature is decreased, which demonstrates the possible presence of two relaxation processes corresponding to the high-frequency signals (fast relaxation phase) and the low-frequency signals (slow relaxation phase) (Supplementary Fig. S7). As shown in the Supplementary Figs. S8 and S9, the Cole-Cole plots show double-separated relaxation phases at low temperature and similar semicircles at high temperature, respectively. By fitting the data with a generalized Debye model47 above 4.0 K, the α values were obtainted in the range of 0.25–0.56, and the energy barrier can be further estimated using the Arrhenius law, τ = τ0 exp(Ueff/kT), giving the parameters Ueff = 29.0 K and τ0 = 3.6 × 10−6 s (R = 0.9965) (Supplementary Fig. S10 and Table S2). While the data of Cole-Cole plots can not be successfully fitted below 4.0 K for double relaxation process with an extended Debye model7 due to the fact that the low frequency peaks of χ″ signals were not clearly shown within the available frequency range of the ac measurements.

Scanning tunneling microscopy (STM) investigations

As complex 1 has the interesting non-centrosymmetric structure, it may have intrinsic molecular dipole and respond to the external electronic field on surface. To better explore the conformation and assembling of this novel complex, STM investigations were performed since such technique has been successfully used to probe the adlayers of several triple-decker sandwich compounds141548495051. In our system, the sole triple-decker complex 1 cannot be stably adsorbed on HOPG surface, as characterized by STM. Nevertheless, with the help of a supramolecular template formed by 1,3,5-tris(10-carboxydecyloxy)benzene (TCDB), complex 1 can be immobilized into the formed nano cavity. On HOPG, through the hydrogen bonding interactions between molecules, TCDB can form rectangular networks with the inner cavity size of 2.8 nm × 1.7 nm (Supplementary Fig. S7). Besides, the shape and cavity size can be tunable based on the guest molecules because of the flexibility of long alkyl chains in TCDB molecule. In our previous work, this supramolecular network has been widely used as molecular templates to host alien ensembles and as spatial confiners to control the chemical reactions52. For our present system, as shown in Fig. 3, at a mixing molar ratio of 1:2 (TCDB: complex 1), one TCDB cavity can accommodate one complex 1 molecule, leading to form a TCDB/complex 1 binary structure. When the molar ratio is increased to 1:4 (TCDB: complex 1), each TCDB network cavity can hold two triple-decker molecules, forming the TCDB/2 complex 1 assembly (Fig. 4).

Figure 3

(a, c) High resolution STM images of TCDB/complex 1 scanned at 700 mV and –500 mV, respectively. Tunneling parameter: Iset = 225.1 pA. (b, d) Suggested molecular models for the scanned structures in (a) and (c), respectively. (e, f) Tentative scheme models for illustrating the conformation of the complex 1 on HOPG surface at an applied positive bias and negative bias, respectively. To better illustrate the two different sides of the complex adsorbed on HOPG, we label the Pc part in purple and one H2L unit in green (all H atoms are omitted). The linear arrows in the models represent the direction of dipolar and the electric property for the triple-decker complex 1. The curved arrow displays the rotation of the top H2L in 1-phenyloctane solution.

Figure 4

(a, c) High resolution STM images for TCDB/2 complex 1 binary structure with the mixing molar ratio of 1:4. (a) Scanned at positive bias. Iset = 232.5 pA, Vbias = 743.1 mV. (c) Scanned at a negative bias. Iset = 232.5 pA, Vbias = −624.3 mV. (b) Suggested molecular models for the scanned structures in (a). (d) Suggested molecular models for the scanned structures in (c). (e, f) Tentative scheme models for illustrating the conformation of two complex 1 molecules on HOPG surface at an applied positive bias and negative bias, respectively. The linear arrows in the models represent the direction of dipolar and the electric property for the triple-decker complex 1. The curved arrow displays the rotation of the top H2L in 1-phenyloctane solution.

Interestingly, the complex 1 can perform the selective adsorption on HOPG, depending on the applied electrical field. At a positive bias scanning condition (Figs. 3a and 4a), close-packed triple layers with square-shaped bright spots can be found on HOPG. On the basis of the attribution of tunneling current, the bright spots should be attributed to the triple-decker. In the present case, since the size (m × n = (1.4 ± 0.1) nm × (1.4 ± 0.1) nm) of the top part is similar to that of the individual Pc molecule53, we suggest that the Pc units point toward the solution while the H2L toward the HOPG surface, that is to say, the conformation of the triple-decker complex 1 at positive bias is the top Pc and the bottom H2L. If the STM measurement is performed at the negative bias, there appears another assembly situation. Under the negative bias, we could observe that the assembled structure of the binary TCDB/complex 1 (Fig. 3c) was disordered, and particularly existing some defects marked by the white arrows in the TCDB/2 complex 1 structure (Fig. 4c). Furthermore, the shape of the bright spots has changed into the rectangular strips, with the measured size of m′ × n′ = (1.6 ± 0.1) nm × (1.3 ± 0.1) nm, which is the same as the assembled structures of TTF54. Hence, it may be highly possible that the H2L locates at the top and the Pc at the bottom, when performed at negative bias.

As displayed in Figs. 3e and 4e, when the STM tip was applied at a positive bias, it is rational that the H2L part is stably adsorbed toward HOPG surface, with the dipole direction along the external applied field. Besides, because of the perfect symmetry of Pc molecule, the top conformation in this case is regularly assembled. By contrary, at the negative bias, the triple-decker complex 1 is apt to be adsorbed with the Pc molecule toward the HOPG surface (Figs. 3f and 4f). Importantly, since the top H2L is a non-symmetric structure, it may allow free rotation in the upper space. Hence, for the binary structures, the top bright spots observed at a negative bias are irregularly arranged. It should be noted that we did not observe any defect for the top conformation in the TCDB/complex 1 structure, which may result from the enough space for the free rotation of each H2L of complex 1 in TCDB cavity. But for the TCDB/2 complex 1 structure, the rotation of the top H2L of the two complex molecules may be different, resulting in the contact and collision of top H2L of the two complex molecules in a confined network during the rotation process, and thus losing some top H2L parts.

Discussion

Based on these observed phenomena, the density functional theory (DFT) calculations have been performed to further investigate the target-molecule self-assembled structures on HOPG. The calculated lattice parameters for the TCDB/complex 1 self-assemblies are summarized in Table 1. It is clear that the calculated parameters of the theoretical models agree well with the experimental values. Meanwhile, we calculated the total energy and total energy per unit area for the TCDB/complex 1 system as shown in Table 2. The total energy includes the interaction between adsorbates (TCDB and complex 1), and the interaction between adsorbates and graphite. Therefore, we could compare thermodynamic stability of the different networks by the total energy per unit area. Evidently, the total energies per unit area of the four kinds of assemblies are nearly equal. It means that the TCDB/complex 1 assembly and TCDB/2 complex 1 assembly, from the view point of thermodynamic, could be co-adsorbed under the same condition, which is consistent with the STM observation.

Table 1

Experimental (Expt.) and calculated (Cal.) lattice parameters for the TCDB/complex 1 self-assemblies

		Unit cell parameters
Assembled system	Adsorbed type	a (nm)		b (nm)	α (deg)
TCDB/complex 1	Pc up	Expt.	4.9 ± 0.1	3.2 ± 0.1	74 ± 2
		Cal.	4.75	3.15	74
	H2L up	Expt.	4.8 ± 0.2	2.8 ± 0.1	73 ± 2
		Cal.	4.50	2.90	74
TCDB/2 complex 1	Pc up	Expt.	5.0 ± 0.1	1.9 ± 0.1	82 ± 2
		Cal.	4.92	2.00	85
	H2L up	Expt.	4.7 ± 0.1	2.0 ± 0.1	83 ± 2
		Cal.	4.80	2.00	85

Table 2

Total energies (Etotal), and total energy per unit area for the observed TCDB/complex 1 self-assemblies. The total energy includes the interaction between adsorbates (TCDB and complex 1) and the interaction between adsorbates and graphite. Here, the more negative energy means the system is more stable

Assembled system	Adsorbed type	Total energies (kcal mol−1)	Total energy per unit area (kcal mol−1 nm−2)
TCDB/complex 1	Pc up	−635.52	−44.19
	H2L up	−642.13	−51.19
TCDB/2 complex 1	Pc up	−494.73	−50.47
	H2L up	−499.84	−52.27

Furthermore, DFT calculations are performed to understand why the different patterns of complex 1 were recorded under the positive/negative STM bias though total energies per unit area of TCDB/complex 1 assemblies nearly equals that of TCDB/2 complex 1 assemblies. To explore the origin of such electrical field dependent selective adsorption, we investigate the interaction between the external electric field and intrinsic molecular properties. DFT results show that the triple-decker complex 1 is with intrinsic molecular dipole P (4.94 Debye), and the dipole moment directs from the Pc to the H2L. When extra electrical field was applied, the extra stability energy Eenergy could be calculated with the equation Eenergy = −P ·Eexternal, which is caused by the interaction between the intrinsic molecular dipole P and the applied electric field Eexternal. Considering that the distance between the STM tip and the substrate surface is 1 nm, the interaction between the molecular dipole and the external electric field could be calculated to be around 0.06 to 0.15 eV (with an applied bias ranging from ±0.60 V to ±1.50 V). In comparison with the thermal energy kT at room temperature (about 0.03 eV), this energy is significant and is considered sufficient to induce selective adsorption of the molecule in the electric field.

In summary, the unique triple-decker dinuclear lanthanide sandwich-type complex based on phthalocyanine and Schiff base ligand containing redox-active TTF unit is successfully synthesized, which has an interestingly non-centrosymmetric structure with two central Dy(III) ions situating in two different coordination geometries. The spectroscopic, electrochemical, and magnetic properties are fully investigated. More importantly, the assembling structures of this triple-decker compound on HOPG surface have been studied by STM, showing a bias dependent selective adsorption. By density function theory calculations, we reveal the formation mechanism of the molecular assemblies, and indicate that such electrical field dependent selective adsorption may be regulated by the interaction between the external electric field and intrinsic molecular properties including molecular dipole, electronics, and symmetry. The synthetic strategy, reported in this paper of introducing TTF unit into the sandwich-type Ln-Pc system, provides a possible route for preparing new multifunctional molecular materials. Moreover, it can be envisioned that such electrical-field-induced selective assembly will be of interest for preparing electrical molecular devices.

Methods

Details of the synthesis and characterization methods of compound 1 was described in the Supplementary Information. The crystallographic data and the bond lengths and bond angles are listed in Supplementary Tables S3 and S4 in the Supplementary Information. CCDC reference number 964624 (1).

TCDB and complex 1 were dissolved in 1-phenyloctane with concentration less than 10−4 M. A droplet (0.4 μL) of the 1-phenyloctane solution containing TCDB/complex 1 mixture (1:2 or 1:4) was deposited onto a freshly cleaved surface (5 mm × 5 mm) of HOPG, respectively. A few minutes later, the sample was studied by STM at a positive bias. After imaging the ordered structures, the bias was changed to a negative direction, and the STM images were captured again. It should be noted that it takes a few minutes to obtain stable assembled structures when changing the bias, this may result from the adsorption/desorption process of the complex 1 at surface.

The STM measurements were performed on a Nano IIIa scanning probe microscope system (Bruker, USA) under ambient conditions. All STM images were recorded in constant current mode using a mechanically cut Pt/Ir (80/20) tip.

Periodic theoretical calculations are performed using density functional theory (DFT) provided by the DMol3 code55. The periodic boundary conditions (PBC) is used to describe the 2D periodic structure on the graphite in this work. The Perdew and Wang parameterization of the local exchange correlation energy are applied in local spin density approximation (LSDA) to describe exchange and correlation56. We expand the all-electron spin-unrestricted Kohn–Sham wave functions in a local atomic orbital basis. For the large system, the numerical basis set is applied. All calculations are all-electron ones, and performed with the medium mesh. Self-consistent field procedure is done with a convergence criterion of 10−5 a.u. on the energy and electron density.

Author Contributions

J.-L.Z. and Q.-D.Z. contributed to the conception and design of the experiments, analysis of the data and revised the paper. F.G. and L.C. carried out the synthetic experiments and characterization. K.D. performed the theoretical computation and analysis. X.-M.Z. designed and carried out the STM experiments, analyzed the data. F.G. and X.-M.Z. wrote the paper.