Noncovalent dyads of lanthanide nitride cluster fullerenes Ln3N@C80 and bisphthalocyanines LnPc2: Insights from DFT calculations
Open Access
- 24 July 2025
- journal article
- research article
- Published by Springer Nature in Journal of Molecular Modeling
- Vol. 31 (8) , 1-11
Abstract
Context: Lanthanide-based systems, such as nitride cluster fullerenes Ln3N@C80 and bipthalocyanines LnPc2 (Pc = phthalocyanine ligand), are of interest for their magnetic, fluorescent and electronic properties. In this regard, we performed DFT characterization to investigate the changes in structure and electronic properties for noncovalently interacting lanthanide (Ln; where Ln = La, Ce, Gd and Lu) nitride cluster fullerenes and bisphthalocyanines to form Ln3N@C80 + LnPc2 dyads. The optimized geometries, formation and frontier orbital energies, HOMO-LUMO plots, charge and spin of Ln and N(NCF) atoms, as well as spin density plots of the dyads were analyzed in comparison with those of isolated Ln3N@C80 and LnPc2 components. In addition to LnPc2 bending distortion, the noncovalent dyad formation alters the geometry of the encapsulated Ln3N cluster, favoring more planar or pyramidal geometries, depending on the case. The HOMO and LUMO orbitals are found on bisphthalocyanines, being localized on the isoindole units, except for Ce3N@C80 + CePc2 dyad, where the LUMO was found on the central metal of CePc2. The HOMO-LUMO gap energy is lower for the dyads compared to isolated NCFs, being rather close to the gap energy of bisphthalocyanines. The changes in spin density distribution are evident in the dyads containing Ce and Gd atoms, contrary to their La and Lu-derived counterparts. The interaction of Ce3N@C80 and Gd3N@C80 with CePc2 and GdPc2, respectively, causes redistribution of the spin density, with changes in the orientation of spin-up and spin-down electrons in the encapsulated Ce3N and Gd3N clusters. Methods: The geometry optimization and electronic properties calculations based on density functional theory were performed using the DMol3 module of Material Studio 8.0 software package from Accelrys Inc. The computational parameters selected included the general gradient approximation functional PBE, combined with a long-range dispersion correction developed by Grimme (PBE-D2), the double numerical basis set (DN), equivalent to the 6-31G Pople-type basis set along with the DFT semiconductor pseudopotentials. To mitigate the self-consistent field convergence problems, the thermal smearing technique was applied, with a final very small value of 0.0001 Ha (equivalent to 31.6 K temperature), or Fermi orbital occupancy in some cases. GraphicalKeywords
This publication has 56 references indexed in Scilit:
- Magnetic behaviour of TbPc2 single-molecule magnets chemically grafted on silicon surfaceNature Communications, 2014
- Endohedral FullerenesChemical Reviews, 2013
- Sandwich Double‐Decker Lanthanide(III) “Intracavity” Complexes Based on Clamshell‐Type Phthalocyanine Ligands: Synthesis, Spectral, Electrochemical, and Spectroelectrochemical InvestigationsChemistry – A European Journal, 2012
- Supramolecular spin valvesNature Materials, 2011
- Graphene Spintronic Devices with Molecular NanomagnetsNano Letters, 2011
- Observation and electric current control of a local spin in a single-molecule magnetNature Communications, 2011
- Anchoring of Rare-Earth-Based Single-Molecule Magnets on Single-Walled Carbon NanotubesJournal of the American Chemical Society, 2009
- Electronic Structure of Surface-supported Bis(phthalocyaninato) terbium(III) Single Molecular MagnetsNano Letters, 2008
- Metal Nitride Cluster Fullerenes: Their Current State and Future ProspectsSmall, 2007
- Semiempirical GGA‐type density functional constructed with a long‐range dispersion correctionJournal of Computational Chemistry, 2006