Amrit
Sarmah
*ab and
Pavel
Hobza
ab
aInstitute of Organic Chemistry and Biochemistry of the Czech Academy of Sciences, Flemingovo nam. 2, CZ-16610 Prague 6, Czech Republic
bDepartment of Physical Chemistry, Palacký University, CZ-77146 Olomouc, Czech Republic. E-mail: amrit.sarmah@marge.uochb.cas.cz; Tel: +420 731015016
First published on 30th July 2020
The directly-linked iron–diporphyrin complexes are appealing candidates and fundamental precursors for an extended metalloporphyrin array that can potentially mimic the biological design of energy-harvesting materials. This encouraged us to appraise the layout for the modular fusion of two iron-porphyrin units. Herein, DFT-based calculations suggest that the electronic environment of diporphyrin systems can be tuned according to the topological attachment between the porphyrin units. Subsequently, a gradual increase in the electronic interaction between the constituent porphyrin units triggers a decrease in the HOMO–LUMO gap. This is essential to achieve higher electric conductivity. The spin-polarized electronic transmission is another interesting aspect of these iron–diporphyrin systems and is promising for spintronic applications. The successive theoretical interpretation of the existence of two-dimensional (2D) metalloporphyrin arrays could be the route to design a graphene analog of the covalent metal–organic framework.
The unparalleled electronic behavior along with their optical linear and non-linear properties, the π-extended porphyrinoid molecules came to prominence and opened up new horizons for technology and medical applications.11–13 In terms of sophisticated technological applications, the π-extended porphyrins are also potential precursors for the fabrication of electron- and energy-transfer devices. Certain important measures that must be tuned for the design and development of these functional porphyrin-based dyes are the nature of the bridges between aromatic components, distance and relative orientation between each chromophore, structural conformation, and the extent of steric interactions.14 It is important to note that in comparison to native porphyrins, π-extended porphyrins display substantial changes in their optical and electrochemical properties. The implications of the covalent linker between essentially coplanar chromophores improve direct conjugation, which induces split Soret and intensively red-shifted Q bands.15
The operational synthetic gateway to sculpt multiporphyrin arrays has been realized with different types of short linkers that are suitable for preparing linear or extended architectures via either meso or β position attachment.16 The fused porphyrin arrays show exciting technological advantages such as rigid and co-planar structures and low HOMO–LUMO gaps. In principle, several different approaches have been adopted to extend porphyrin chromophores without changing the actual number of pyrrole units and methine bridges. The outcome of a practical implementation of this idea is tetrabenzoporphyrin.17–19 The increment in every single benzene unit included in the macrocycle results in a bathochromic shift of about 30 nm. Subsequently, Lash et al. reported a synthetic strategy for the extension of the conjugated aromatic unit in the porphyrin array.20 Along these same lines, an important breakthrough was the preparation of a cross-shaped dye comprising five conjugated porphyrin units reported by Vicente and co-workers.21 In another interesting study, Crossley et al. reported the laboratory synthesis of porphyrin arrays containing one or more pyrazine units.22 Thompson and co-workers, based on their prior experiences on synthetic protocols of fused porphyrins, concluded that meso–meso, β–β, β–β diporphyrin-bearing pyrene moieties are appropriate starting substrates for intramolecular oxidation reactions.23
The feasibility of tetrameric metalloporphyrin sheets has recently been realized,24 indicating the possible existence of 2D infinite metalloporphyrin sheets in experimental conditions. For the first-time, Abel et al. reported the fabrication of a novel two-dimensional (2D) Fe-embedded phthalocyanine (poly-FePc) organometallic porous sheet.25 These newly categorized 2D materials exhibit promising conduction and magnetic properties; however, they are restricted to variable grain boundary conditions. In these 2D systems, the TM (transition metal) atoms are uniformly and separately distributed without clustering and display well-defined geometries and magnetic properties. Subsequently, different TM-Pc monolayers find sophisticated applications in spintronics26–29 and gas capture.30,31 Indeed, due to the unique and regular structural orientation of the TM atoms in the TM-Pc monolayers, these 2D architectures are appealing candidates for single-atom catalysts.32 In this direction, DFT-based theoretical modeling predicted that not-yet-synthesized Co-Pc and Cr-Pc monolayers can be applied as good catalysts for CO oxidation.33,34 However, to date, the electronic properties of the basic building blocks of the 2D metalloporphyrin sheets, i.e., the porphyrin-dimer units, have scarcely been addressed.
In this study, we have extended a systematic investigation to develop a better understanding of the electronic behavior of fused metalloporphyrin moieties based on density functional theory (DFT). In analogy with previous reports, significant modulation in the electronic properties of the hybrid system can be achieved depending on the fusion mode between porphyrin units. Subsequently, we also traced the relative changes in the electronic properties of the systems during the transformation from one-dimensional to two-dimensional arrays, i.e., sheet-like structures. The manuscript is organized in the following way. In Section 2, we address the computational methodologies. The results and discussion in Section 3(a) consist of an analysis of the DFT results for fused porphyrins at the molecular level along with their photophysical behaviors. The essence of electronic transport through the porphyrin dimers in different fusion modes is assayed in Section 3(b). A comprehensive theoretical examination of the electronic behavior of the 2D arrays is featured in Section 3(c). Finally, we summarize the key findings of our theoretical simulation along with some futuristic applications of the model systems in the concluding remarks in Section 4.
The periodic calculations were executed using spin-polarized density functional theory (DFT) under the generalized gradient approximation (GGA)45 in the form envisaged by Perdew, Burke, and Ernzerhof (PBE) as deployed in the Vienna ab initio simulation package (VASP).46 The projected augmented wave (PAW) method47,48 with a plane-wave basis set was used. Subsequently, the spin-polarized calculations were performed under the Vosko–Wilk–Nusair (VWN) modification49 scheme to interpolate the correlation energy. We enforced periodic boundary conditions and a vacuum space of 15 Å along the z-direction in order to avoid interactions between two FeP (iron-porphyrin) images in nearest-neighbor unit cells. All of the structures were relaxed using the conjugate gradient method without any symmetric constraints. To investigate the magnetic coupling between the TM atoms, we considered a 2 × 2 surface unit cell. A 3 × 3 × 1 Monkhorst–Pack50k-point grid centered on the Γ point was used for the calculations. The energy cutoff and convergence criteria for energy and force were set to 400 eV, 1 × 10−4 eV, and 0.01 eV Å−1, respectively.
Here, system 1 consists of a porphyrin dimer linked through a 1,3-butadiyne bridge. Similarly, triply and doubly linked attachment features were explored with system 2 and system 4. Again, a fully conjugated porphyrin dimer was attributed to system 3. The aromatic character of porphyrin, which contains 26 delocalized π electrons, is well established.52 The 1,3-butadiyne bridge between the metalloporphyrin units forms the planar Fe–PP dimeric structure. The planar structural arrangement may increase the possibility that both the porphyrin and butadiyne moieties will contribute to the extended π-electron conjugation over the entire system. Subsequently, the extended π-electron conjugation also improves the electrical conduction performance.53 Higher consistency between experimental findings and calculated values for the planar porphyrin arrays compared to their non-planar analogs is evident from a recent report,54 which encouraged us to concentrate on the planar systems. Similarly, the higher degree of extended π-electron conjugation in the other three fused Fe–PP dimer systems also opens up new avenues to tune and manipulate the characteristic electronic properties of the systems by means of structural orientation. The linked M–PP (M = Cr, Mn, Fe, Zn, Cu) dimer is the simplest building-block for one-dimensional metalloporphyrin nanowires.55 Embedding 3d TM (transition metal) atoms, i.e., Fe, into the porphyrin cores can significantly modulate the electronic properties of the system. The Fe centers are in +2 valence states51 and are coordinated to the nitrogen atoms. The combination of strong coordination interactions and relatively large distances between the two Fe centers (due to the rigid covalent bridging) effectively prevent the aggregation of these atoms. As a result, the Fe atoms can be incorporated separately into the individual porphyrin units, which can endow the Fe–PP dimer with intrinsic magnetic properties; therefore, they are promising probes in the fabrication of nanodevices for spintronics applications.
The computed total density-of-states (TDOS) plots along with the location of occupied and virtual orbitals for the four systems are reported in Fig. 2. The green and red lines indicate the occupied and virtual orbitals. The blue line in the plots represents the TDOS of the system.
As can be observed from Fig. 2, the electronic states of the first three systems have close similarities with one another except for small variations near the HOMO level. However, the fourth molecule exhibits significant changes in its electronic states around the HOMO level. A systematic analysis of the plots provides useful insights into the electronic structure. The significant electronic contribution from the p-type atomic orbitals of carbon and nitrogen atoms to the frontier molecular orbitals (FMOs) of the system is evident from the isosurface plots reported in Fig. 2. Subsequently, the corresponding PDOS plots for the model systems are included in the ESI.† The substantial contribution of the iron d-orbitals to the total density of states along with the bonding nature of the MOs constituted from their overlap with the carbon p-orbitals near the HOMO region can be distinguished from the plots. It is important to note that the electronic charge densities on the iron centers at the FMO levels are consistent for all the systems. However, there are some observed deviations from this ideal nature of the metallic centers at the HOMO of the first and the LUMO of the fourth systems. The visual graphics of the corresponding isosurface plots for the frontier molecular orbitals (FMO) (Fig. 3) further confirmed the above observations. Subsequently, we also performed calculations with a higher spin-multiplicity of the systems. The existence of unpaired spins on the central metal atom is clearly evident from the spin-density maps reported in the ESI.†
Fig. 3 Visual graphics of the frontier molecular orbitals (FMO) of the systems. The values in the boxes represent the HOMO–LUMO gaps. |
The gradual increase in the electronic interaction between the constituent porphyrin units triggers a decrease in the HOMO–LUMO gap. This is essential to achieve higher electric conductivity. Subsequently, a large number of bridging structures that can enhance the inter-porphyrin conjugation have been intensively exploited.56–58 The 1,3-butadiyne-bridged porphyrin oligomers, i.e., system 1, showed positive cooperative self-assembly (i.e. all-or-nothing assembly) as well as moderate increases in planarity and conjugation, resulting in reduction in the HOMO–LUMO band gaps. A higher degree of conjugation in system 2 facilitated greater electron delocalization, which minimizes the HOMO–LUMO gap. The preferential σ-bonding characteristics in the highest occupied MO levels of systems 1 and 2 are evident from Fig. 2. However, the π-bonding and π-antibonding natures of the HOMO and LUMO orbitals, respectively, in system 3 are associated with an increase in the HOMO–LUMO energy gap. A seemingly different FMO picture is observed for system 4. Here, the HOMO and LUMO are primarily located at the metal centers and are associated with an extremely small electronic bandgap; this indicates the profound metallic nature of the molecule, which correlates well with the DOS plot in Fig. 2.
So-called “fused diporphyrin” is an excellent and straightforward synthetic protocol to connect two porphyrin units directly through multiple covalent bonds. Within a coplanar geometrical arrangement, the systems should be quite favorable for electronic π conjugation. As we have seen, diporphyrin is the fundamental unit of an extended 1D metalloporphyrin array. These 1D structures have compelling electronic and magnetic properties and are considered as molecular wires or porphyrin-taps by virtue of their rigid shape and extended electronic conjugation.59 It is now well understood that the extraordinary electronic properties of metalloporphyrin arrays are a consequence of the modular sequence of the diporphyrin units. Therefore, it is important to develop a better understanding of the chemical nature as well as the electronic structure of the ‘linking-zone’ between two single units (Fig. 4).
In the preceding sections, we addressed the findings of different theoretical tools to account for quantitative recognition of the chemical bonding in the diporphyrin systems. It would be useful to obtain some qualitative interpretations of the results in terms of visual graphics. The localized-orbital locator (LOL) is an emergent descriptor of chemical bonding.60 According to the general convention of LOL, the gradients of the localized orbitals are maximized with the increasing extent of overlap in the localized orbitals. Additional theoretical details about LOL are included in the ESI.† Moreover, LOL has a simpler interpretation in terms of fast and slow electron regions. Qualitatively, the high/low LOL values are associated with slow/fast electrons and slow electrons are attributed to the characteristic of localized electrons, such as those in bonds or lone pairs. The light red and yellow regions in the 1,3-butadiene bridge represent two weak covalent C–C π-bonds along with the strong σ-bonding character of the C–C single bond. The bright red spot in the map depicts higher electron localization due to the presence of the C–C single bond. The projection shows spikes around the atomic zone in the 3D surface and distinguishes the blue circles in the 2D map. There are some distinct differences in the nature of chemical bonding for the fused diporphyrin systems. Apparently, the benzofused (system 3) and triple-bridge (system 2) complexes exhibit aromatic ring current due to the fast π-electron cloud at the center of the six-membered ring (the blue region). However, the α–β bridge complex (system 4) does not have this aromatic character, although it has some effective π-conjugations at the four-membered ring. Again, the benzofused and triple-bridge systems show similar bonding characteristics; however, a critical analysis of the bonding pattern on the basis of the LOL map at the fused position reveals certain differences. The regular distributions of blue and red zones for the triple-bridge system are indications of fully conjugated behavior of the π-electrons, and the regular valley-shaped 3D projection also supports this observation. On the other hand, the benzo-fused system has a lower degree of electron delocalization, as evident from the sharp spikes in the 3D projection along with the non-uniform distribution of blue and red zones in the 2D color map. The anti-aromatic character of the meso–meso bridge system was also realized from the LOL map.
System | Peak maxima (nm) | Oscillator strength | Major orbital contributions |
---|---|---|---|
System 1 | 618 (S0 → S7) | 1.41 | HOMO → LUMO (94%) |
449 (S0 → S30) | 1.06 | HOMO−1 → LUMO+4 (27%) | |
System 2 | 806 (S0 → S8) | 0.208 | HOMO−1 → LUMO+1 (83%) |
523 (S0 → S22) | 1.43 | HOMO → LUMO (76%) | |
System 3 | 548 (S0 → S12) | 0.570 | HOMO → LUMO+1 (83%) |
445 (S0 → S31) | 1.1927 | HOMO−1 → LUMO (81%) | |
System 4 | 1137 (S0 → S17) | 0.3081 | HOMO → LUMO (93%) |
560 (S0 → S40) | 1.8699 | HOMO−5 → LUMO+2 (60%) |
It is observed that in the 1,3-butadiene-linked diporphyrin complex (system 1), the low-energy Q-band is a sharper peak compared to those of the other three systems, where the Q-band appears to be a broad peak. The Q-band at λmax = 618 nm for system 1 is mainly the result of a HOMO → LUMO type of electronic transition with a significantly high oscillator strength.60 Similarly, the Q-bands for systems 2 and 4 are largely red-shifted (λmax = 806 and 1137 nm, respectively) due to the increase in delocalized π-electrons and the enlargement of the conjugate surface. However, in the benzo-fused system 3, the Q-band is blue-shifted by 70 nm relative to that of system 1; this is the consequence of a strong interporphyrin interaction mediated by the benzenoid unit.61,62 The substantial Q-band peak broadening in the fused porphyrin complexes is induced by the relative increase in the electronic interaction between the porphyrin units. The characteristic Soret band of the porphyrin is observed at λmax = 449 and 445 nm for systems 1 and 3, respectively, which is associated with their high oscillator strength. However, this particular band appears to be exceedingly red-shifted to λmax = 523 and 560 nm for systems 2 and 4, respectively. The red-shifted absorption bands associated with the fused structures suggest the presence of unusually large electronic interactions between the π orbitals on the two macrocycles. The unobtrusive and efficient intramolecular charge transfer (CT), including the CT from the metal to the ligand, the CT from the ligand to the metal, or the intramolecular CT within ligands is evident from the electron density difference map. A small amount of π–π* of electron transition can be clearly observed, which is the main reason for the large nonlinear optical (NLO) response of organic molecules.
Initially, the molecule was deposited on a Ag(111) surface, and total energy DFT calculations were performed with the VASP package to obtain the relaxed geometry of the adsorbed molecule. The adsorbed molecule was found to be located around 2.75 Å above the metal surface, and it appears to remain planar. More details about the DFT calculations are included in the ESI.† To the best of the authors’ knowledge, one-to-one correlations between the simulated images and experimental data for all the model systems included in our study are non-viable. However, we explored various experimental reports with similar porphyrin motifs to account for the relevance of our theoretical modeling.
In the reported literature, the 8-hydroxyquinoline tetramer has been extensively investigated from both experimental and theoretical aspects. For this particular system, high-resolution AFM and STM images based on experiments as well as DFT calculations are well documented.65,66 The bond contrast in the AFM images has been qualitatively compared to the bond order, where a higher local electron density leads to stronger Pauli repulsion exerted on the tip. Subsequently, it is conclusive to extend qualitative interpretations to the nature of the existing bonding patterns in the system. It is evident from the study by Zhang et al.66 that the electron localization function (ELF) is directly related to the contrast of a chemical bond appearing in the AFM image. The computed ELF maps of the four systems are depicted in Fig. 7. A significantly localized electron density donation from N to Fe in the ELF of the dimers is evident from the graphics, whereas the electrons between the β–β and meso–meso position bonding are rather delocalized. These characteristics are attributed to the fact that the bonding features in the central regions of the dimer are dative and the fused or linked positions are covalent in nature. On the basis of ELF, it can be argued that the Fe–N bond is significantly polarized and that the shared electrons are preferably localized near the N atom; thus, the AFM contrast is negligible at the metal center.
Inelastic electron tunneling spectroscopy (IETS) has re-emerged as an efficient analytical tool for understanding nanoscale and molecular junctions. Subsequently, the mapping of the intensity of the IETS feature can contribute to the understanding of the spatial distribution of the functional groups. According to our theoretical model, a CO molecule is attached to the tip, and the bonding between two atoms in the porphyrin dimer can be sensed by a particular vibrational mode of the tip CO. When we perform a gradual scanning of the adsorbed molecule with the CO-terminated tip, there should be a significant modulation in the energy and intensity of the hindered translational vibration of CO, which can be traced by inelastic electron tunneling spectroscopy (IETS) with STM.67 A systematic account of the relative variation in low-energy CO vibration can manipulate the potential energy landscape of the molecule and its surroundings. The primary advantage of the IETS-STM technique is the real-space spectroscopy and imaging of the probe. The AFM technique basically detects the local density of states (LDOS), whereas IETS based on STM has proven to be a versatile technique for the real-space imaging of molecular structures and chemical bonding of single molecules on surfaces. Comparison of IETS spectral features between experiments and simulations offers detailed information about chemical conformation, electron–vibration coupling, etc.
The comprehensive discussions of these two sophisticated analytical techniques in terms of computer simulation provide useful insights into the chemical structures of porphyrin dimers.
The transport calculations were performed based on the NEGF + DFT formalism as implemented in the ADF-BAND periodic DFT code.72,73 Comprehensive details about the calculations are available in the ESI.† This is indeed an efficient tool to investigate the electronic transport behavior of a novel single-molecular junction with metallic electrodes. In our present study, we modeled a gold–(porphyrin dimer)–gold (Au–PorD–Au) junction. The contacts were designed with Au(111) faces consisting of 3 × 3 atoms, as represented in Fig. 8. The calculations were performed at the PBE level with an SZ basis on the Au contacts and a TZP basis on the molecule, except where stated otherwise. More theoretical details are available in the ESI.†
The foremost observation from the plots reported in Fig. 8 is that there are substantial variations in the electronic transmission with the changes in the position as well as the orientation of attachment between the two porphyrin units. Another important observation is that the manifestation of spin-polarized electronic transmission, i.e., electronic transmission through the two spin-channels, is different. The extensive research work of Harmann et al. on spin-resolved transport through organometallic complexes revealed important implications of electron transport in nanoscale devices.74 In a molecular bridge, the decrease in conductance due to the appearance of destructive quantum interference (DQI) favors ferromagnetic coupling between two spin centers, while increased conductance through the bridge can favor antiferromagnetic coupling.75 It was observed that in the case of a π-conjugated system, destructive interference could block the electronic transmission, contrary to partial cancellation for the saturated system.76 Similarly, the local transmission through the complexes and their molecular orbitals shows substantial dependencies on the various components of the molecular structure that mediate electron transport.77 It is worth mentioning here that the electronic behavior of a system around the Fermi level determines the versatility of that molecule for nanoscale device fabrication. The present study accounts for a significant modulation of electronic transmission at the two spin states around the Fermi level for all the systems. From Fig. 8, it can be assumed that the origin of spin-dependent electronic transport in the system is due to the spin asymmetry of the electronic structure caused by the presence of paramagnetic metal centers in each porphyrin unit. The findings of spin-polarized calculations in the next segment suggested the presence of spin asymmetry in the model systems. This suggests that the spin-filtering capacity is inherent in the hybrid systems, not in the interface. All systems show the metallic nature of transmission associated with a sharp kink around the Fermi level. The first three systems exhibit an intense peak for the down-spin electrons close to the Fermi zone, whereas the fourth system defines the possibility of a sharp increase in the transmission through the up-spin channel at the Fermi level. The shapes and positions of these transmission peaks signify low-bias electronic conductance for these junctions. The analysis of the NEGF–DFT-based findings in terms of transmission plots uncovers certain important interpretations. Another striking feature is the quantum interference effects in the electronic transport of these systems.78–81 It appears that the electronic transmission through the spin-down channel exhibits significant quantum interference around the Fermi level. The transmission spectra exhibit a smooth regression around the Fermi energy for spin-up (red line in Fig. 8), while for spin-down (black line in Fig. 8), the transmission exhibits a sharp dip near the Fermi energy; this suggests a destructive pattern of quantum interference.82 The interference features are dramatic and in many cases dominate the characteristic transmission spectrum around the Fermi level. Here, the quantum interference patterns are different according to the topological changes in the model systems. It is worth mentioning here that the interference features arise from antiresonances (sharp dips in the transmission). The interference pattern, an antiresonance, is exclusively associated with the down-spin electronic current and affects the system transmission close to the Fermi energy. Antiresonances are simply the opposite of transmission resonances; consequently, the transmission coefficient goes to zero. One condition for the appearance of antiresonances in the transmission is the presence of nonspanning nodes in the wave function. Subsequently, these analyses also indicated some qualitative similarities in the transmission behavior of the four systems at the Fermi level, such as spin-polarized electronic transmission, distinguished broad peaks for the up-spin electrons, and sharp kinks for the down-spin electrons. This information enhances our conceptual understanding of the effects of structural orientation on the electronic transport phenomena. Further, we can argue that the electronic transmission through the junction can be tuned by varying the position as well as the orientation of attachment between the two porphyrin units, which is crucial for strategic development in the fabrication of nanoscale devices based on linked or fused porphyrin dimers.
In light of realistic application, the molecular design should be extended further to periodic models. Previous reports encouraged us to explore the possibilities of viability and the electronic structures of the 2D sheet-like orientation of fused porphyrin-based materials. Taking an analogy from earlier results, intuitively, we designed three novel 2D structures for the porphyrin oligomers. The single unit cell of these three systems is reported in Fig. 9(a). Here, we considered the unit cell to be composed of four porphyrin molecules attached in a different fashion, which should be the minimal unit to account for the relative electronic modulations. The significant differences in the electronic environments of the systems are evident from their corresponding spin-polarized band structures reported in Fig. 9(b). As can be seen from the band diagrams, although four metal centers are present in the unit cell, the systems are slightly magnetic. This is due to the strong anti-ferromagnetic interaction between the metal centers.90 The electronic behaviors of the three systems are entirely different from one another. A critical analysis of the band structure reveals that system 1 is semi-conducting in nature with a specific band-gap. Similarly, system 2 is seemingly half-metallic; one spin-channel is in the metallic region, and the other one is insulating. These characteristic electronic properties of TM-Pc monolayers affirm the promising aspects of the material for spintronics.91 Eventually, system 3 shows typical metallic behavior, with the accumulation of flat bands around the Fermi level. To obtain some better insights into the electronic levels and to trace the contribution of the Fe metal centers, we mapped the total and partial densities of states (TDOS and PDOS) for the systems, and the plots are reported in Fig. 10.
The computed TDOS shows almost symmetric distributions with some distinct deviations around the Fermi level. This correlates well with the lower magnetic moment values of the system. The contribution of electronic states of the metal centers toward the TDOS around the Fermi level (the yellow and purple shades in the plots) leads to the existence of spin-polarity in the systems. Additionally, the introduction of metal centers induced disturbance of the spin states of the p orbitals, and the plot became unsymmetrical. In the case of system 1, the electrons from iron and d and nitrogen p orbitals occupy the spin-down states just below the Fermi level, and no occupied state can be found across the Fermi level; therefore, system 1 is a semiconductor. Subsequently, system 2 shows half-metallic nature, as the p-electrons from nitrogen in the up-spin channel spread across the Fermi level and the iron d-electrons in the down-spin channel located just below the Fermi level create a gap for the down-spin channel. Meanwhile, system 3 exhibits metallic nature, with significant contributions from the iron d-electrons in the up-spin channel in the occupied states, and the down-spin channel is mostly occupied by the p-electrons from nitrogen (Fig. 11).
Fig. 11 Comparison between experimental and simulated STM images of the three systems generated based on the Tersoff–Hamann method92 within the DFT framework. (a) The relaxed structures of the three different systems, (b) relevant experimental STM images reported in previous studies, and (c) simulated STM maps; the purple box represents the molecular zone that appeared in the graphics. |
The most efficient approach to obtaining a direct correlation between theoretical atomistic models and experimental systems is scanning tunneling microscopy (STM). With the rapid progress in modern ab initio methods, it is now possible to accommodate ever-larger structures. Subsequently, the constant improvement of experimental STM leading to higher resolution images (i.e. smaller structures). Due to its simplicity, the Tersoff–Hamann (TH) method has been incorporated into most existing DFT codes. Although the TH method cannot explain the origin of atomic resolution in STM measurements, it often provides a reasonably good qualitative prediction of the experimental observations.93 On the other hand, there are certain complications with the implementation of the NEGF approach for the theoretical analysis of STM.94 The main obstacle is positioning the STM tip at an arbitrary offset with respect to the surface plane. Due to this, together with the fact that tip and sample materials are highly inhomogeneous, it is highly problematic to represent the combined substrate-tip interface within the same supercell.95 Therefore, we are currently in a position to make some sensible comparisons between experimental and calculated STM images based on the Tersoff–Hamann (TH) method. The experimental STM map has many advantages; however, it is blind to the atomic species observed in the STM images. This weak point can be covered up by comparison to the calculated STM images. If the calculated and experimentally observed STM images show close similarities, one can very reasonably assume that the underlying geometry is the same. To account for the relevance of our simulated STM images to the experimental reality, we obtained sophisticated experimental results from some recent reports.
Zhang and co-workers studied the Sonogashira cross-coupling reactions on an Au(111) surface with three different metal-coordination structures and used scanning tunneling microscopy to resolve the products formed at the successive reaction steps.96 The observed STM images of the butadiene bridge Cu–porphyrin oligomers show reasonable agreement with our simulated map. Further, the STM images reported by Nakamura et al.86 and Grill et al.97 in their respective studies manifested excellent correlation with our calculated STM surfaces for similar systems. The comprehensive electronic structure calculations in combination with the extensive analytical simulations strongly suggest the experimental feasibility of linked metalloporphyrin two-dimensional (2D) arrays and are assumed to be the porphyrin analogs of graphene.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0ma00461h |
This journal is © The Royal Society of Chemistry 2020 |