# Nonuniform STM Contrast of Self-Assembled
Tri-n-octyl-triazatriangulenium Tetrafluoroborate on
HOPG

Sergii Snegir, Yannick J Dappe, Dmytro Sysoiev, Thomas Huhn, Elke Scheer

► To cite this version:

Sergii Snegir, Yannick J Dappe, Dmytro Sysoiev, Thomas Huhn, Elke Scheer. Nonuniform STM Contrast of Self-Assembled Tri-n-octyl-triazatriangulenium Tetrafluoroborate on HOPG. ACS Omega, 2023, 8, pp.38766 - 38772. 10.1021/acsomega.3c06454. hal-04803288

HAL Id: hal-04803288

https://cnrs.hal.science/hal-04803288v1

Submitted on 25 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

# Nonuniform STM Contrast of Self-Assembled Tri-*n*-octyl-triazatriangulenum Tetrafluoroborate on HOPG

Sergii Snegir, Yannick J. Dappe, Dmytro Sysoiev, Thomas Huhn, and Elke Scheer*

<figure><img src="image_4.png" /></figure>

Cite This: https://doi.org/10.1021/acsomega.3c06454

<figure><img src="image_5.png" /></figure>

Read Online

## ACCESS

<figure><img src="image_6.png" /></figure>

Metrics & More

<figure><img src="image_7.png" /></figure>

Article Recommendations

<figure><img src="image_8.png" /></figure>

Supporting Information

**ABSTRACT:** We have assembled 4,8,12-tri-*n*-octyl-4,8,12-triazatriangulenum-tetrafluoroborate (TATA-BF₄) on highly oriented pyrolytic graphite (HOPG) and have studied the structure and tunneling properties of this self-assembled monolayer (SAM) using scanning tunneling microscopy (STM) under ambient conditions. We show that the triazatriangulenum cations TATA⁺ form hexagonally packed structures driven by the interaction between the aromatic core and the HOPG lattice, as evidenced by density functional theory (DFT) modeling. According to the DFT results, the three alkyl chains of the platform tend to follow the main crystallographic directions of HOPG, leading to a different STM appearance. The STM contrast of the SAM shows that the monolayer is formed by two types of species, namely, TATA⁺ with BF₄⁻ counterions on top and without them. The cationic TATA⁺ platform gives rise to a seemingly higher appearance than neutral TATA-BF₄, in contrast to observations made on metallic substrates. The variation of the STM tunneling parameters does not change the relative difference of contrast, revealing the stability of both species on HOPG. DFT calculations show that TATA-BF₄ on HOPG has sufficient binding energy to resist dissociation into TATA⁺ and BF₄⁻, which might occur under the action of the electric field in the tunneling gap during STM scanning.

<figure><img src="image_9.png" /></figure>

## INTRODUCTION

Triazatriangulenum (TATA) salts make up a class of very stable fluorescent dyes with a cationic organic core. TATA⁺ cations are classically synthesized by nucleophilic displacement of all methoxy groups of tris-2,6-dimethoxyphenyl methanol. This is done using a primary *n*-alkyl amine, followed by elimination of water by aid of an acid.¹ Depending on the nature of the acid, TATAs with different counterions are accessible.² Upon reaction with suitable nucleophiles, un-charged TATA species with the nucleophile covalently bound to the central carbon atom are accessible.³ Inspired by this idea, the TATA⁺ core has been functionalized with different azobenzene derivatives, and the platform was found to self-organize on Au surfaces.³,⁴ The high symmetry and π-conjugated planar structure of cationic as well as covalently bound TATAs allowed the formation of extended two-dimensional (2D) supramolecular assemblies on the surface.⁴⁻⁹ Moreover, by variation of the peripheral *n*-alkyl chain length, the platform surface density and thus the distance of the attached functional units have been adjusted.¹⁰

Despite the significant progress of using TATA⁺ as a platform, the unique properties of the cationic core, stabilized by different counterions, have not yet been well explored. Some limited number of studies showed that combining cationic TATA⁺ salts with anion receptors leads to columnar self-assembly of molecules in solution and thus to the formation of organogels consisting of submicrometer-scale

fibers.² In these structures, counterions (BF₄⁻, Cl⁻, or Br⁻) play the role of a binding agent, forming the cationic planar TATA⁺ aromatic core with a receptor anion. After deposition of TATA-BF₄ on a free-electron rich Au surface, the counterions are obscured.¹⁰ Neither X-ray photoelectron spectroscopy (XPS) nor electrochemical STM in a HClO₄ electrolyte could detect BF₄⁻ moieties on the surface.¹¹ Instead, the authors attribute this behavior to a strong interaction of the monolayer with the Au(111) substrate, leading to a screening of the positive charge of the cation and removal of the corresponding anion after adsorption of the monolayer.⁵,¹¹,¹² Recently, we succeeded in clarifying the whereabouts of the BF₄⁻ anion on Au(111) by heat-aided deposition of TATA-BF₄ from 1,2,4-trichlorobenzene solvent.¹³ The counterions could be localized only under specific tunneling conditions. We detected two types of molecular packing on the Au surface and showed that 15–20% of the TATA⁺ entities within a self-assembled monolayer (SAM) still contained a BF₄⁻ counterion on their top. Some BF₄⁻ species

Received: August 29, 2023

Accepted: September 21, 2023
were found directly on Au(111) trapped in the empty meshes of the SAM. Additionally, density functional theory (DFT) calculations showed that the loss of BF₄⁻ counterions can be explained by the decrease of the strength of the ionic bond between the TATA⁺ core and their counterions due to the interaction of the π-system with the Au lattice.¹³

In the current study, we make efforts to self-assemble TATA-BF₄ from a 1,2,4-trichlorobenzene solution on a heated HOPG substrate. HOPG was used in this study to test our hypothesis that Au would be particularly efficient in promoting the dissociation of TATA-BF₄ during the assembly of the monolayer. We assumed that HOPG would interact less strongly with the π-system of the TATA⁺ core, whereas the interaction of the *n*-alkyl chains with the underlying HOPG lattice would dominate due to the known epitaxy of *n*-alkanes with the graphite lattice.¹⁴⁻¹⁷ Under these conditions, we assumed that the number of self-assembled undissociated TATA-BF₄ entities would be significantly higher compared to that on Au(111).

## EXPERIMENTAL SECTION

Synthesis and purification of TATA-BF₄ with *n*-octyl side chains (R = $-[CH_2]_7-CH_3$) were performed following [1]. The corresponding NMR spectroscopy studies are presented in our recent paper.¹³ The powder was stored at -18 °C under a N₂ atmosphere. All STM measurements were performed with freshly prepared solutions of TATA-BF₄ in 1,2,4-trichlorobenzene. We have used the same concentration of $1.3 \times 10^{-7}$ mL⁻¹ as in our previous study, which avoided the formation of a second molecular layer.⁵,¹²

**STM Measurements.** All SAMs were studied under ambient atmospheric conditions using a commercial STM equipped with a low-current head (Bruker, Nanoscope 3A). The calibration of the STM piezoceramics was controlled in advance using a freshly cleaved HOPG surface. The STM tip was prepared by mechanical cutting of a Pt/Ir (80:20) wire. For each SAM, several STM images were recorded in the constant-current mode with the current set point ranging from 10 to 25 pA and tip bias from 0.1 to 0.4 V. Images were obtained for two different samples using four different tips to check the reproducibility and to ensure that the results are free from artifacts. All molecule-molecule distances given in the paper are averaged from different images and along different crystallographic orientations to minimize errors caused by the thermal drift of the STM tip when working under ambient conditions. Postfiltering of the noise in the obtained images was done with the commercial software package SPIP (Digital Surf).

The solution of TATA-BF₄ was deposited on freshly cleaved HOPG substrates (Micro to Nano, The Netherlands). The substrate with the deposited solution was left on a heated plate (~50 °C) under ambient atmospheric pressure in a fume cabinet before the STM measurements. After about 1 h of drying, a visual inspection of the substrate revealed complete evaporation of the solvent, and the STM images showed the SAM formation.

The statistical analysis of STM images for defining the ratio of TATA⁺ cores to TATA-BF₄ within a monitored SAM was done based on the relative height difference.

**Computational Methods.** The experimental results are complemented by periodic DFT calculations performed with the localized orbital code Fireball.¹⁸ This approach uses a self-consistent version of the Harris-Foulkes local density

approximation (LDA) functional,¹⁹,²⁰ instead of the traditional Kohn-Sham functional based on the electronic density. As a result, the total charge leading to the electrostatic potential of the system is approximated as the superposition of spherical charges around each atom. In the Fireball simulation package, self-consistency is achieved over the occupation numbers through the Harris functional.²¹ We built the optimized basis set using a linear combination of wave functions, one from the atom in its ground state and the other one from the atom in its first excited state. This combination allows us to smoothen the decay of the radial part of the wave function and optimize the overlaps between two neighboring atoms.²² The LDA exchange-correlation energy is calculated using the efficient multicenter-weighted exchange-correlation density approximation (McWEDA).²³,²⁴ In the present work, the cutoff radii for the wave functions' radial part (for *s*, *sp*, and *spd* basis sets) defining the optimized basis set are the following (in atomic units): $r_s = 4.1$ for the H atom; $r_s = 4.5, 4.2$, and $4.1$; $r_p = 4.5, 4.2$, and $4.1$ for the C, N, and F atoms, respectively; $r_s = 4.5$; $r_p = 4.9$; and $r_d = 5.2$ for the B atom. These optimized basis sets have been well tested in previous works on molecular self-assembly on metallic surfaces.¹³,²⁵,²⁶ Periodic calculations have been performed using a set of $8 \times 8 \times 1$ k-points in the surface Brillouin zone of graphene. Finally, we have used a well-established perturbation theory approach combined with the DFT formalism in order to take into account a weak repulsive Coulombic and attractive van der Waals interaction for structural optimization of the molecule adsorption on a single layer of graphene or in the TATA network.²⁷⁻²⁹ Because of van der Waals interactions between the graphene layers, the distance between the TATA and the second layer would be around 6 Å, resulting in only a very weak interaction with no significant influence on the adsorption energy.

Using this approach, we determined the adsorption energy of the TATA⁺ ion on graphene. The adsorption energy is defined in a standard manner as the difference between the total energy of the system and the energy of the isolated ion and the isolated graphene surface, both calculated independently. To this end, we have considered a 14 × 14-graphene unit cell with the TATA⁺ core on top. The core and its arms were fully optimized before setting them on the graphene surface. Then, the whole structure has been optimized until the forces went below $0.1 \text{ eV} \cdot \text{Å}^{-1}$. The interaction energy of the molecule with the surface was calculated as a function of the molecule-surface distance, taking into account weak repulsive and van der Waals interactions. From the obtained equilibrium structure, we have determined the projected density of states (PDOS) of the different species involved in this study in comparison with their PDOS when adsorbed on a gold surface. Finally, we determined the two-dimensional network of TATA platforms.

## RESULTS AND DISCUSSION

**Experimental Observations.** Figure 1 represents a 66.4 × 66.4 nm² area of the HOPG surface with two large terraces. The top left terrace with a smooth STM contrast has no indication of a molecular assembly, whereas the right bottom terrace is almost completely covered by tiny spots forming a periodic pattern. Each of these spots corresponds to the position of a single TATA (see below). The black dashed curve designates the border between the SAM and the disordered region. The observed SAM in Figure 1 forms a single domain since it has rotational symmetry observed at 120°. The
<figure><img src="image_10.png" /><figcaption>**Figure 1.** STM image (132.8 × 132.8 nm²) of a SAM of TATA-BF₄ on HOPG (STM parameters: U_t = 0.25 V and I_t = 12 pA). The black dashed curve highlights the border between the areas with and without the SAM. The red and blue ellipses are explained in the text. The height color code is used for all STM images throughout the study.</figcaption></figure>

observed domain contains several STM contrast peculiarities that are different in nature. Two elongated broad lines with a bright (red markup) and dark (blue markup) STM contrast highlight the surface area where the STM tunneling conditions are not uniform due to a lattice mismatch of the underlying graphite sheets. This stems from either a structural defect or a slight misorientation of a lower lying HOPG layer with respect to the top layer. The resultant wave functions in this area adopt either a maximum or a minimum, leading to the formation of a bright or dark STM contrast, respectively. Despite this effect, the structure of the SAM in these regions remains uniform.

The analysis of the STM contrast of Figure 1 revealed that entities within a SAM may have a different STM appearance. Figure 2a shows an enlarged area of Figure 1 free from distinct defects of the HOPG layers. The image shows a part of a SAM with irregular STM contrast. The cross sections A–B and C–D along the arbitrarily selected rows of several neighboring bright spots show ~0.4 Å larger relative height compared to others (Figure 2b), revealing higher local electron tunneling probability. It was counted that among 529 entities imaged in Figure 2a, 128 are showing a relative height between 0.2 and 0.4 Å, suggesting that about 24% of arbitrarily located species show a brighter contrast. We explain the slight variation of the relative height of TATA⁺ (brightness of spots in Figure 2) by disorder in the HOPG substrate. This disorder leads to a small variation of the interaction between the two topmost graphene layers of HOPG and hence a variation of the effective charge of TATA⁺, which is crucial for the STM contrast as we explain in the Discussion section. Scanning the same surface area using different STM tunneling parameters did not change the positions of the bright species. However, upon scanning, local distortions in the vicinity of these species became apparent under certain tunneling parameters.

By variation of the tunneling parameters, different molecular orbitals are probed. The resultant STM contrast of the species underwent a significant change from a large uniform spot (U_t = 0.25 V and I_t = 12 pA, Figure 2a) to clearly distinguishable ring-like structures (U_t = 0.3 V and I_t = 15 pA, Figure 2c) and finally to a combination of interconnected elongated objects (U_t = 0.3 V and I_t = 14 pA, Figure 2d). Cross-sectional analysis of all the observed 2D structures in Figure 2 under the above STM tunneling conditions allowed for reliable identification of a hexagonal arrangement of TATAs within a SAM with a

lattice constant of *a* = 2.2 ± 0.1 nm. The minor variation of the current between Figure 2c,d reflects the shape of some molecular orbitals (MOs) which contribute to the current. Although several MOs can be energetically close to each other, their shapes may differ. By a slight decrease in the tunneling current from 0.15 to 0.14 pA, the distance between the STM tip and the sample is increased. This alters the overlap between the wave functions of the metal tip and the MOs of the sample.

**Theoretical Modeling.** In order to determine the adsorption energy of TATA on HOPG, we placed the TATA⁺ ion in different positions on a graphene plane, considering the alignment of one alkyl chain along the [110] direction.<sup>14–17</sup> The graphene plane was used as the model of the top layer of HOPG. As such, we initially have considered a configuration of the TATA⁺ core being rotated with respect to the graphene hexagonal rings by 30°, then a rotation of TATA⁺ by additional 30° to align the core with the hexagonal rings of the graphene layer, and finally a translation to give an AB-like stacking with respect to the graphene layer. The latter one turns out to be the most favorable position and is represented in Figure 3a. This configuration has been calculated after full DFT optimization and taking into account van der Waals interactions, as described above. The corresponding adsorption energy curves as a function of the molecule–graphene distance are represented in Figure 3b. The obtained energy lies around

<figure><img src="image_11.png" /><figcaption>**Figure 2.** (a) STM image (50.5 × 50.5 nm²) of a SAM of TATA-BF₄ deposited on HOPG. The corresponding STM parameters are U_t = 0.25 V and I_t = 12 pA. (b) Cross sections along the corresponding directions in (a). (c,d) STM images (11.3 × 11.3 nm²) of an area of (a) recorded under U_t = 0.3 V, I_t = 15 pA and U_t = 0.3 V, I_t = 14 pA, respectively. The white shaded diamonds in (c,d) indicate the hexagonal unit cell with lattice parameter *a*. The height color code is the same as that in Figure 1.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 3.** (a) Top and side view of the DFT-optimized configuration of the TATA⁺ ion adsorbed on graphene. The color code is dark gray for carbon, blue for nitrogen, white for hydrogen, and light gray for carbon atoms in graphene. (b) Evolution of the adsorption energy for different adsorption configurations of TATA⁺ on graphene. (c) Evolution of the molecule-molecule interaction energy in the hexagonal molecular network as a function of the intermolecular distance.</figcaption></figure>

3 eV per TATA⁺ ion for an equilibrium distance of 1.6 Å between the lowest hydrogen atom of the ion and the graphene plane, i.e. 3.32 Å, between the aromatic core of the ion and the graphene plane (indicated in the bottom of Figure 3a). This corresponds to the usual equilibrium distance between graphitic materials.³⁰ In order to characterize the experimentally observed TATA molecular network on HOPG and to compare the molecule-graphene and molecule-molecule interactions, we have also theoretically modeled the potential network by constructing a hexagonal unit cell of 7 TATA⁺ ions, as represented in Figure 3c. We have then calculated the interaction energy between the central ion and its next neighbors in the cell, for different intermolecular distances, using the previous formalism. As a result, we find an equilibrium distance of 21.5 Å between two neighboring ions, corresponding to the van der Waals equilibrium distance. This is in good agreement with the experimental observations. The corresponding interaction energy variation is represented in Figure 3c.

**Discussion.** Figure 4 shows a small part (8.9 × 8.9 nm²) of Figure 2d with the superimposed chemical structure of n-octyl-TATA⁺. According to DFT calculations, the structure of the molecule has a distorted star-like geometry because it allows gaining a minimum of energy on graphene (Figure 3a). BF₄⁻ was intentionally removed from the structure to simplify the figure. The aromatic cores of the ion structures were placed directly above the surface areas that under the monitored tunneling conditions appeared as ring-like objects (Figure 2c,d). A similar ring-like contrast was observed for TATA on Au(111), reflecting positions of the molecular orbital of the

<figure><img src="image_13.png" /><figcaption>**Figure 4.** STM image (8.9 × 8.9 nm²) with superimposed DFT-optimized configuration of TATA. The STM parameters are Uₜ = 0.3 V and Iₜ = 14 pA.</figcaption></figure>

central cores of the TATA entities.¹³ However, the STM contrast of the molecular core on HOPG has peculiarities that are observed under Uₜ = 0.3 V and Iₜ = 14 pA. We attribute these peculiarities to a strong epitaxy of the aromatic TATA⁺ core with the underlying HOPG lattice (Figure 3a). Our DFT calculations strongly favor an orientation of the six hexagonal rings of the TATA⁺ core with respect to the graphene lattice. This is similar to the case of two adjacent graphene planes in a consecutive AB pattern by translation along one C–C bond (see above).¹⁴⁻¹⁶ Hence, the carbon atoms of the TATA⁺ core are located in an alternating fashion above the center of a graphene hexagon and above the carbon atom of the
underlying graphene. The latter coorientation creates specific tunneling conditions, allowing STM visualization of overlapping atoms of TATA⁺ and graphene (Figure 3a). Three of these atoms in the TATA⁺ core are nitrogen atoms, which might lead to the observed variation in STM contrast compared with the situation with two overlapping carbon atoms. Therefore, this effect can explain the presence of the observed bright spots in the STM contrast (white rings in Figure 4). By application of this hypothesis, the exact position of the aromatic core of TATA⁺ on the STM image (Figure 4) can be explained as well as the exact positions of the alkyl chains. The alkyl chains lie right above the elongated 2D objects in the recorded STM image (Figure 4). At first sight, it seems counterintuitive that the alkyl chains that are known to be “invisible” for STM on the Au(111) surface have distinguishable STM contrast on HOPG. However, it is also known that alkanes tend to orient themselves along [110] crystallographic directions and hence tend to form ideal epitaxial structures on a HOPG lattice.<sup>14–17</sup> Such epitaxy leads to the appearance of periodically repeated protrusions due to created resonant tunneling conditions, allowing visualization of the entire length of the alkyl chains.<sup>14–16</sup> This situation holds true only for two of the three alkyl chains, whereas the third chain remains invisible due to a slight misfit with the HOPG lattice. From DFT calculations (Figure 3a), it becomes apparent that two alkyl chains at least partially lie above the carbon atoms along [110], while the third alkyl chain cannot accommodate a position along or even close to the [110] direction. Apparently, the interaction of the aromatic core of TATA⁺ with graphene outweighs the simultaneous epitaxy of all the three alkyl chains on graphene. From these results, we deduce that the adsorption configuration of TATA⁺ on HOPG is imposed mainly by the aromatic core of TATA⁺, whereas the alkyl chains are occupying favorable positions on the lattice by tilting their N–C–C angle and/or by bending the chain (see the DFT-optimized structure in Figure 3a).

At the same time, the packing density of the SAM is defined exclusively by the length of the alkyl chains. TATA platforms cannot be packed any closer because of the steric hindrance between neighboring entities (Figure 3c). This leads to a minimum intermolecular interaction energy of 0.1 eV/mol at a separation distance of 21.5 Å ≈ 9 × T<sub>HOPG</sub> (T<sub>HOPG</sub> = 2.46 Å). The analysis of STM images recorded under different tunneling conditions (given in Figure 2a–d) shows that the same self-assembly mechanism is active for both species with normal and brighter STM contrast.

One open question so far is whether the bright species [with ~0.4 Å bigger relative height for the tunneling conditions of Figure 1 compared to “dark” molecules] correspond to TATAs that still bear their BF₄⁻ counterions on top or to naked TATA⁺ ions instead.

To address this question, we calculated the DOS for TATA⁺ and TATA-BF₄ on graphene in the configuration shown in Figure 5a,b respectively. The DOS for the TATA⁺ ion shows a pronounced resonance at the Fermi level E<sub>F</sub> corresponding to the highest occupied molecular orbital (HOMO) of the ion. It is located close to E<sub>F</sub> because of the weaker van der Waals interaction of TATA⁺ with the substrate, which prevents the ion from being fully neutralized from graphene. In contrast, the HOMO of the neutral molecule TATA-BF₄ is located more than 1 eV below E<sub>F</sub>. We argue that the enhanced DOS at E<sub>F</sub> of TATA⁺ may result in an enhanced electronic transmission probability to the STM tip and, hence, in a brighter STM

<figure><img src="image_14.png" /><figcaption>**Figure 5.** Geometrical configuration of (a) TATA⁺, (b) TATA-BF₄ on graphene, (c) TATA⁺ on Au(111), and (d) TATA-BF₄ on Au(111). (e) Calculated PDOS for the four configurations shown in panels (a–d).</figcaption></figure>

appearance. The absolute value of the transmission critically depends on the position of the molecule on the substrate as well as the distance between the TATA⁺ and the BF₄⁻ counterion and the distance to the STM tip. These parameters are not known in the experiment, hindering us from performing meaningful transmission calculations.

We hence posit that the bright protrusions in the STM figure do NOT show the undissociated TATA-BF₄ molecules but instead the naked TATA⁺ ions. Counting these latter species, we find a fraction of roughly 24% bright molecules, meaning that 76% of the molecules stay undissociated. This interpretation is additionally supported by the slightly varying STM appearance of the bright spots when applying larger voltages. As can be seen in Figure 5e, the full width at half-maximum (fwhm) of the resonance at E<sub>F</sub> of TATA⁺ is about 0.4 eV. Increasing the tunneling voltage, U<sub>v</sub> from 0.25 to 0.3 V slightly changed the appearance of the bright species, while the dark species remain unaffected. In our earlier study of TATA-BF₄ on Au(111),<sup>13</sup> we observed the opposite effect in both aspects: undissociated TATA-BF₄ appeared as bright spots and with a relative abundance of 15–20%. An important difference between the two systems is that the bright spots reflect different aspects of HOPG and Au(111). While on Au(111) they can be seen as topographic contrast (higher elevation due to the presence of the counterions), they signal electronic effects in the case of HOPG. For completeness, we show in Figure 5e also the DOS for TATA⁺ and TATA-BF₄ on Au(111), where in both cases, the HOMO is well below E<sub>F</sub> and broadened, indicative of the partial charge transfer due to the strong interaction with the metal substrate.

As mentioned in the Introduction section, in our earlier study of TATA-BF₄ on Au(111), we were able to observe the split-off BF₄⁻ counterions in the meshes of the TATA-SAMs.<sup>13</sup> On HOPG, we do not detect free counterions, most likely because of their smaller amount and weak interaction with
HOPG, which might enhance their mobility. Moreover, in the current study, all bright spots stay constant during several STM scans and during application of different tunneling bias and currents, which was not the case for TATA-BF₄ on a Au surface. Such observations support the idea that the bright species correspond to the naked TATA⁺ platforms.

Since about a quarter of all species lose their counterions, one may ask at which moment in the sample preparation they are cleaved off the platform. One possibility would be that part of the molecules predissociate in the solution and can adsorb on HOPG in both predissociated (TATA⁺) and undissociated (TATA-BF₄) forms. The final formation of the SAM leads to the stabilization of either TATA-BF₄ or TATA⁺ species on HOPG. Obviously, more integral methods such as XPS might bring additional insights to answer the question. However, the self-assembly of TATA-BF₄ does not lead to the formation of sufficiently large domains for a distinct XPS study. Also, particular approaches applied to promote the SAM formation were insufficient to detect BF₄⁻ on the surface by XPS.¹¹

## CONCLUSIONS

A SAM of TATA-BF₄ on HOPG is formed by mainly van der Waals interaction between the aromatic core of TATA⁺ and the substrate. The 6-fold rotational symmetry of graphene (top layer of HOPG) partially imposes the alkyl chains to orient with respect to the graphene lattice, forming a hexagonal molecular packing, revealed by a star-like structure of the entities on the surface. The presence of BF₄⁻ counterions on the majority of the TATA⁺ platforms does not impair the substrate-molecule interaction and therefore does not distort the hexagonal packing of the SAM. Undissociated TATA-BF₄ can be easily detected on HOPG using different STM tunneling conditions. The main reason for this behavior is that the electric field in the tunneling gap is insufficient to overcome the dissociation energy of TATA-BF₄.

## ASSOCIATED CONTENT

### Supporting Information

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acsomega.3c06454.

Details on theoretical modeling of TATA⁺ and TATA-BF₄ on HOPG including the result for different stacking geometries, information on noise filtering of some of the STM images, details of the calculation of the projected density of states, and structural data of the adsorbed molecules (PDF)

## AUTHOR INFORMATION

### Corresponding Author

Elke Scheer – Department of Physics, University of Konstanz, Konstanz 78464, Germany; orcid.org/0000-0003-3788-6979; Email: elke.scheer@uni-konstanz.de

### Authors

Sergii Snegir – Department of Physics, University of Konstanz, Konstanz 78464, Germany; orcid.org/0000-0002-4024-5956

Yannick J. Dappe – SPEC, CEA, CNRS, Université Paris-Saclay, Gif-sur-Yvette Cédex 91191, France; orcid.org/0000-0002-1358-3474

Dmytro Sysoiev – Department of Chemistry, University of Konstanz, Konstanz 78464, Germany

Thomas Huhn – Department of Chemistry, University of Konstanz, Konstanz 78464, Germany; orcid.org/0000-0001-6292-4215

Complete contact information is available at: https://pubs.acs.org/10.1021/acsomega.3c06454

### Notes

The authors declare no competing financial interest.

## ACKNOWLEDGMENTS

The authors appreciate the financial support of the Deutsche Forschungsgemeinschaft (DFG) through the project Plasmo-chrom (no. 406778771).

## REFERENCES

(1) Laursen, B. W.; Krebs, F. C. Synthesis of a Triazatriangulenum Salt. Angew. Chem., Int. Ed. 2000, 39, 3432–3434.

(2) Haketa, Y.; Sasaki, S.; Ohta, N.; Masunaga, H.; Ogawa, H.; Mizuno, N.; Araoka, F.; Takezoe, H.; Maeda, H. Oriented Salts: Dimension-Controlled Charge-by-Charge Assemblies from Planar Receptor-Anion Complexes. Angew. Chem., Int. Ed. 2010, 49, 10079–10083.

(3) Kubitschke, J.; Näther, C.; Herges, R. Synthesis of Functionalized Triazatriangulenes for Application in Photo-Switchable Self-Assembled Monolayers. Eur. J. Org Chem. 2010, 2010, 5041–5055.

(4) Jung, U.; Kuhn, S.; Cornelissen, U.; Tuczek, F.; Strunskus, T.; Zaporojtchenko, V.; Kubitschke, J.; Herges, R.; Magnussen, O. Azobenzene-Containing Triazatriangulenum Adlayers on Au(111): Structural and Spectroscopic Characterization. Langmuir 2011, 27, 5899–5908.

(5) Kuhn, S.; Baisch, B.; Jung, U.; Johannsen, T.; Kubitschke, J.; Herges, R.; Magnussen, O. Self-assembly of triazatriangulenum-based functional adlayers on Au(111) surfaces. Phys. Chem. Chem. Phys. 2010, 12, 4481–4487.

(6) Kuhn, S.; Jung, U.; Ulrich, S.; Herges, R.; Magnussen, O. Adlayers based on molecular platforms of trioxatriangulenum. Chem. Commun. 2011, 47, 8880–8882.

(7) Jung, U.; Schütt, C.; Filinova, O.; Kubitschke, J.; Herges, R.; Magnussen, O. Photoswitching of Azobenzene-Functionalized Molecular Platforms on Au Surfaces. J. Phys. Chem. C 2012, 116, 25943–25948.

(8) Otte, F. L.; Lemke, S.; Schütt, C.; Krekiehn, N. R.; Jung, U.; Magnussen, O. M.; Herges, R. Ordered Monolayers of Free-Standing Porphyrins on Gold. J. Am. Chem. Soc. 2014, 136, 11248–11251.

(9) Hammerich, M.; Herges, R. Laterally Mounted Azobenzenes on Platforms. J. Org. Chem. 2015, 80, 11233–11236.

(10) Lemke, S.; Ulrich, S.; Claußen, F.; Bloedorn, A.; Jung, U.; Herges, R.; Magnussen, O. M. Triazatriangulenum adlayers on Au(111): Superstructure as a function of alkyl side chain length. Surf. Sci. 2015, 632, 71–76.

(11) Ulrich, S.; Jung, U.; Strunskus, T.; Schütt, C.; Bloedorn, A.; Lemke, S.; Ludwig, E.; Kipp, L.; Faupel, F.; Magnussen, O.; Herges, R. X-ray spectroscopy characterization of azobenzene-functionalized triazatriangulenum adlayers on Au(111) surfaces. Phys. Chem. Chem. Phys. 2015, 17, 17053–17062.

(12) Baisch, B.; Raffa, D.; Jung, U.; Magnussen, O. M.; Nicolas, C.; Lacour, J.; Kubitschke, J.; Herges, R. Mounting Freestanding Molecular Functions onto Surfaces: The Platform Approach. J. Am. Chem. Soc. 2009, 131, 442–443.

(13) Snegir, S.; Dappe, Y. J.; Sysoiev, D.; Pluchery, O.; Huhn, T.; Scheer, E. Where do the counterions go? Tip-induced dissociation of self-assembled triazatriangulenum-based molecules on Au(111). Phys. Chem. Chem. Phys. 2021, 23, 9930–9937.

(14) Leunissen, M. E.; Graswinckel, W. S.; van Enckevort, W. J. P.; Vlieg, E. Epitaxial Nucleation and Growth of n-Alkane Crystals on Graphite (0001). Cryst. Growth Des. 2004, 4, 361–367.
(15) Groszek, A. J. Selective adsorption at graphite/hydrocarbon interfaces. *Proc. R. Soc. London, Ser. A* **1970**, *314*, 473-498.

(16) Claypool, C. L.; Faglioni, F.; Goddard, W. A.; Gray, H. B.; Lewis, N. S.; Marcus, R. A. Source of image contrast in STM images of functionalized alkanes on graphite: A systematic functional group approach. *J. Phys. Chem. B* **1997**, *101*, 5978-5995.

(17) Faglioni, F.; Claypool, C. L.; Lewis, N. S.; Goddard, W. A. Theoretical Description of the STM Images of Alkanes and Substituted Alkanes Adsorbed on Graphite. *J. Phys. Chem. B* **1997**, *101*, 5996-6020.

(18) Lewis, J. P.; Jelínek, P.; Ortega, J.; Demkov, A. A.; Trabada, D. G.; Haycock, B.; Wang, H.; Adams, G.; Tomfohr, J. K.; Abad, E.; Wang, H.; Drabold, D. A. Advances and applications in the FIREBALL ab initio tight-binding molecular-dynamics formalism. *Phys. Status Solidi B* **2011**, 248, 1989-2007.

(19) Harris, J. Simplified method for calculating the energy of weakly interacting fragments. *Phys. Rev. B* **1985**, *31*, 1770-1779.

(20) Foulkes, W. M. C.; Haydock, R. Tight-binding models and density-functional theory. *Phys. Rev. B* **1989**, *39*, 12520-12536.

(21) Demkov, A. A.; Ortega, J.; Sankey, O. F.; Grumbach, M. P. Electronic structure approach for complex silicas. *Phys. Rev. B* **1995**, *52*, 1618-1630.

(22) Basanta, M. A.; Dappe, Y. J.; Jelínek, P.; Ortega, J. Optimized atomic-like orbitals for first-principles tight-binding molecular dynamics. *Comput. Mater. Sci.* **2007**, *39*, 759-766.

(23) Jelínek, P.; Wang, H.; Lewis, J. P.; Sankey, O. F.; Ortega, J. Multicenter approach to the exchange-correlation interactions in ab initio tight-binding methods. *Phys. Rev. B* **2005**, *71*, 235101.

(24) Sankey, O. F.; Niklewski, D. J. Ab initio multicenter tight-binding model for molecular-dynamics simulations and other applications in covalent systems. *Phys. Rev. B* **1989**, *40*, 3979-3995.

(25) Snegir, S.; Dappe, Y. J.; Kapitanchuk, O. L.; Coursault, D.; Lacaze, E. Kinked row-induced chirality driven by molecule-substrate interactions. *Phys. Chem. Chem. Phys.* **2020**, *22*, 7259-7267.

(26) Sleczkowski, P.; Dappe, Y. J.; Croset, B.; Shimizu, Y.; Tanaka, D.; Minobe, R.; Uchida, K.; Lacaze, E. Two-Dimensional Self-Assembly Monitored by Hydrogen Bonds for Triphenylene-Based Molecules: New Role for Azobenzenes. *J. Phys. Chem. C* **2016**, *120*, 22388-22397.

(27) Basanta, M. A.; Dappe, Y. J.; Ortega, J.; Flores, F. Van der Waals forces in the local-orbital Density Functional Theory. *Europhys. Lett.* **2005**, *70*, 355-361.

(28) Dappe, Y. J.; Martínez, J. I. Effect of van der Waals forces on the stacking of coronenes encapsulated in a single-wall carbon nanotube and many-body excitation spectrum. *Carbon* **2013**, *54*, 113-123.

(29) Gonzalez, C.; Abad, E.; Dappe, Y. J.; Cuevas, J. C. Theoretical study of carbon-based tips for scanning tunnelling microscopy. *Nanotechnology* **2016**, *27*, 105201.

(30) Esfandiarpoor, S.; Fazli, M.; Ganji, M. D. Reactive molecular dynamic simulations on the gas separation performance of porous graphene membrane. *Sci. Rep.* **2017**, *7*, 16561.