4Nanosystem Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki 305-8568, Japan
5Departamento de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, Spain
6Department of Physics, University of Konstanz, D-78457 Konstanz, Germany
Guest Editor: P. J. Skabara
Beilstein J. Org. Chem. 2015, 11, 1068–1078. doi:10.3762/bjoc.11.120
Received 11 Mar 2015, Accepted 22 May 2015, Published 24 Jun 2015
We describe the synthesis and single-molecule electrical transport properties of a molecular wire containing a π-extended tetrathiafulvalene (exTTF) group and its charge-transfer complex with F4TCNQ. We form single-molecule junctions using the in situ break junction technique using a homebuilt scanning tunneling microscope with a range of conductance between 10 G0 down to 10−7 G0. Within this range we do not observe a clear conductance signature of the neutral parent molecule, suggesting either that its conductance is too low or that it does not form a stable junction. Conversely, we do find a clear conductance signature in the experiments carried out on the charge-transfer complex. Due to the fact we expected this species to have a higher conductance than the neutral molecule, we believe this supports the idea that the conductance of the neutral molecule is very low, below our measurement sensitivity. This idea is further supported by theoretical calculations. To the best of our knowledge, these are the first reported single-molecule conductance measurements on a molecular charge-transfer species.
Keywords: break junction measurements; charge-transfer complex; DFT-based transport; molecular electronics; tetrathiafulvalene
The development of molecular electronics is a current challenge in nanoscience. The ultimate goal is to fabricate different electronic devices based on a variety of active elements such as wires, transistors, diodes or switches (to name just a few), where each is built from individual, suitably functionalized molecules. Although the controlled handling of molecules to form reliable molecule-based circuits remains a demanding task, modern organic synthesis allows the design and preparation of nearly any challenging molecule .
Among the aforementioned electroactive elements, the study of molecular wires has received great attention and, in this regard, a great variety of molecules of different nature involving single, double and triple C–C bonds (conjugated or not) have been extensively studied [2-4]. Furthermore, most of these systems have been covalently connected to a great variety of different anchor groups in order to improve the connection to various metal electrodes. However, despite the huge number of molecular systems that find application as wires, the use of electroactive molecules exhibiting different oxidation states that can modify/control the conductance through the wire have been considerably less studied.
In the realm of organic chemistry there are a great variety of organic compounds able to show different redox states that are very appealing candidates to be used as wires and/or switches in molecular electronics. In this regard, molecules having redox centers such as viologen [5,6], aniline [7,8], thiophene , anthraquinone  and ferrocene  have been previously studied. However, a particularly suitable redox-active molecule for molecular electronics is the well-known electron donor tetrathiafulvalene (TTF) molecule. Pristine TTF, as well as the tetraselenafulvalene analogue (TSF), have been previously reported. In this study, the authors hypothesized that in the Au–TTF–Au junctions, the molecule is connected to the electrodes in a face-to-face overlapping configuration . In contrast, since the first report on a suitably functionalized TTF as a molecular wire using two thioacetate anchoring groups , most of the TTF derivatives synthesized for this purpose have been functionalized with sulfur atoms as alligator clips. These belong either to a fused ring on the TTF , or to a chain covalently connected to the TTF core bearing a thiol group at the end [15,16]. Furthermore, extended TTF cruciform molecules, formed by two orthogonally placed, π-conjugated moieties bearing the 1,3-dithiole rings at the ends, have also been used for single-molecule measurements .
A singular TTF analogue is the so-called π-extended TTF (exTTF, (9,10-bis(1,3-dithiol-2-ylidene)-9,10-dihydroanthracene) which, in contrast to pristine TTF that exhibits two oxidation peaks to form the radical cation and dication species, shows only one oxidation peak involving a two electron process to form the dication state. Furthermore, the geometrical properties of exTTF are quite different from pristine TTF. Thus, whereas TTF is mostly planar in the neutral form, exTTF is a highly distorted, out-of-plane molecule with a butterfly shape in its neutral state. It undergoes a dramatic gain of planarity and aromaticity upon oxidation (Figure 1) [18,19]. This gain of stability upon oxidation has been skillfully used in a variety of D–π–A systems, namely exTTF–π-bridge–C60 derivatives, for determining the attenuation factor of the molecular wire (oligomer) acting as the π-bridge [20-22].
In this paper, we describe the synthesis of a new exTTF derivative, suitably functionalized with two (p-acetylthio)phenylethynyl substituents at positions 2 and 6 of the anthracene central core. We then describe single-molecule conductance measurements on this new derivative, along with measurements of the charge-transfer complex formed with F4TCNQ. Finally, we present theoretical calculations to understand its electrical transport properties.
The synthesis of the target molecule is shown in Scheme 1 and starts from the previously reported 2,6-diiodo-exTTF 1 . Reaction of 1 under Sonogashira conditions (Pd(II), CuI, DIPEA) with trimethylsilylacetylene affords the symmetrically substituted exTTF 2 in good yield. Further removal of the trimethylsilyl group is easily achieved by treatment with potassium carbonate, yielding the free terminal alkyne 3 in quantitative yield. The introduction of the two anchor groups in 3 is carried out through a second Sonogashira reaction. The coupling of 2,6-diethynyl-exTTF 3 with 1-acetylthio-4-iodobenzene (4), in the presence of Pd(II) and copper iodide and DIPEA in THF, led to the target molecule 5 in moderate yield (35%). Compound 4 was obtained in one step and with high yield from 1,4-diiodobenzene by reaction with n-BuLi in ether, followed by reaction with S8, and eventually, acetyl chloride .
All new compounds were fully characterized by spectroscopic and electrochemical means. Interestingly, compounds 2, 3 and 5 showed similar 1H NMR spectra due to their symmetry resulting in relatively simple spectra, which confirms the proposed structures. In particular, compound 5 exhibits the methyl groups as a singlet at 2.45 ppm, and the protons corresponding to the 1,3-dithiole rings appear as a singlet at 6.37 ppm. In the 13C NMR of the target molecule 5, the carbonyl groups appear at 192.4 ppm and the carbons of the alkyne moieties at 88.0 and 90.4 ppm, with the terminal methyl groups at 28.7 ppm.
The redox properties of compound 5 (0.2 mM) were determined by cyclic voltammetry at room temperature in THF using TBAPF6 (0.1 M) as a supporting electrolyte under argon atmosphere and at a scan rate of 0.1 Vs−1. The electrochemical cell consisted of a glassy carbon working electrode, Ag/AgNO3 reference electrode and a Pt wire counter electrode. It is worth mentioning that ferrocene was not employed as the inner reference since its oxidation potential overlaps with that of the exTTF unit. Similar to pristine exTTF, the new exTTF derivative 5 exhibited only one quasi-reversible oxidation peak, involving a two-electron process to form the dication state [25,26]. This oxidation peak appears at Epa = 217 mV (ΔE = 285 mV, peak-to-peak separation), which is quite similar to the oxidation peak found for pristine exTTF (Epa = 244 mV, ΔE = 350 mV) (Figure 2).
We first tried to form neutral-molecule molecular junctions. Simply, compound 5 can be seen as an analogue of an oligo(phenylene ethynylene), specifically an OPE3-dithiol compound (where 3 indicates the total number of phenyl rings), in which the central phenyl ring has been substituted by an exTTF unit. We could, therefore, expect compound 5 to form molecular junctions in a similar way to the OPE3-dithiol, which readily does so either in solution, or under solvent-free conditions as recently reported . Hence, as a starting point, we followed the same sample preparation conditions as previously with the OPE3-dithiol . We prepared 0.1–1 mM solutions of compound 5 in both 1,2,4-trichlorobenzene (TCB) or mesitylene/tetrahydrofuran (Mes/THF 4:1), and exposed a clean gold substrate to the solution for a period of approximately 30 min. The sample was then dried under a flow of argon and mounted inside the scanning tunneling microscope (STM). All experiments were then performed under solvent-free, ambient conditions. In order to form molecular junctions of compound 5, we followed the break junction technique . During the experiment, the variation in conductance (G) is recorded while an STM tip is moved vertically (z) in and out of contact with a gold substrate, forming and breaking gold nanocontacts (G vs z trace). When the two gold electrodes (the STM tip and the substrate) are in contact, or in close enough proximity, one or more molecules of the compound adsorbed on the surface can bridge the two electrodes. This binding occurs through one of the terminal groups (thiols in this case) to each of the electrodes, hence forming a molecular junction. In this configuration, when separating the electrodes (larger z values), we observe conductance plateaus while the gold nanocontact or a molecular junction remains intact during the pulling. The plateaus then end abruptly when the junction is broken.
Figure 3c and 3d show the 2D histograms consisting of several thousand individual G vs z traces recorded in a break junction experiment in the presence of compound 5. These are compared with those recorded on an unmodified gold substrate and in the presence of OPE3-dithiol. These histograms include all the measured traces without filtering out the tunneling-only traces from those with plateaus. One can see that the colored region of the 2D histogram of OPE3-dithiol extends significantly to larger inter-electrode distances than the unmodified gold substrate. Also, it does so within a narrow band of conductance values (roughly between log(G/G0) = −3 and −4.5). This 2D profile exemplifies the presence of plateaus occurring in a narrow region of conductance. In the 2D histograms of compound 5, we can also see that the colored region extends to higher electrode interdistance values (z) than for the unmodified gold substrate. However, as opposed to the case of OPE-dithiol, there is no clear protuberance in the histogram that would indicate the concentration of plateaus at a given conductance value.
Therefore, in order to determine whether the signature of compound 5 is simply very weak, we performed a trace separation process to build histograms of only the traces displaying plateaus. Figure 4 shows the results of this separation for the traces considered in Figure 3. In particular, for compound 5, two separation steps were carried out, resulting in three categories of traces: (1) those with only a smooth exponential decay (Figure 4a and Figure 4d, labeled “No Plat.”), (2) those with poorly-defined plateaus (Figure 4b and Figure 4e, labeled “Plat. 1”), and (3) those with well-defined plateaus (Figure 4c and Figure 4f, labeled “Plat. 2”). We considered the conductance range between log(G/G0) = −5.3 and −0.5 and established that a trace has a plateau whenever a z interval (Δz) longer than 0.15 nm is needed to observe a change of conductance Δlog(G/G0) = 0.12 along the trace. Amongst these traces, we separate the well-defined plateaus as those traces for which a Δz larger than 0.2 nm is needed to observe a Δlog(G/G0) = 0.1 along the trace. These are empirical criteria which gave us the best separation results. We see that the histograms built from the trace of the third group present protuberances, suggesting that molecular junctions are successfully formed from compound 5.
From a close inspection of Figure 4c and Figure 4f, we see that the plateau conductance in these experiments varies by more than 2 orders of magnitude (between log(G/G0) = −3 and −5), significantly greater than OPE3-dithiol, and also that they extend to between 1 and 2 nm. However, it is important to note here that the percentage of traces with well-defined plateaus (given in parentheses in Figure 4) is very low. Even with a total number of 5000 measured traces, the final number of traces with well-defined plateaus is less than 50 traces. We also note here that these results correspond to the most successful measurement runs, and the percentage of traces with plateaus was even lower in other cases. We stress that under the same experimental conditions, we obtained percentages of around 35% of well-defined plateaus for OPE3-dithiol.
We have recently measured the conductance of a series of molecules based on a similar OPE3 backbone, also terminated with thiols, but with differing numbers of dithiafulvalene (DTF) substituents placed at various positions along the OPE backbone . One of the main structural differences with these molecules compared to compound 5 is that conjugation is maintained along the OPE3 backbone. For this series of compounds, we could observe a clear signature for each, which was similar to the unmodified OPE3. We showed that the presence of the DTF side groups does not influence the low-bias conductance and secondly, and that their presence does not significantly hinder molecular junction formation. Different to the measurements carried out thus far on compound 5, we carried out these measurements using dichloromethane as the deposition solvent, and since the results were clear, we decided to applied this also to compound 5. The sample was prepared by exposing the gold to a solution of compound 5 in dichloromethane (DCM) for 1 h, followed by drying the sample with a flow of N2. The results of these measurements are shown in Figure 5.
As with the measurements on compound 5 using the TCB drop-casting method, after dip-casting from DCM, we could observe the presence of G(z) plateaus, in this case down to just above 10−7 G0 facilitated by the use of larger gains. The percentage of traces displaying these plateaus was again low, however, and similar to that found before, with less than 2% fitting the well-defined criteria. From the 2D histograms in Figure 4, we see again that no favored region of conductance develops after separating the junctions displaying plateaus (labeled Plat. 1 and Plat. 2). This, we believe, indicates that whilst the formation of molecular junctions is possible, it is unclear exactly how the molecule binds in these junctions. The low percentage of molecular junction formation can be considered consistent with our observations on the OPE3 series with DTF side groups, in that the presence of sulfur atoms in the center of the molecule slightly reduces the probability of forming junctions. This is likely because the mobility of the molecule over the surface is reduced, which in turn makes it difficult for molecules to diffuse into the freshly created nanogaps.
The lack of a clear molecular signature may ultimately be due to several reasons. Aside from reducing the probability of junction formation, the sulfur atoms in the center of the molecule can also bind to the electrodes during the evolution of the junctions, preventing the wire from becoming fully stretched. This is a distinct possibility for this molecule due to the well-known interaction of the parent extended-TTF and gold . An alternative explanation for the lack of a clear signal for compound 5 is that its end-to-end conductance is too low for it to be observed in our setup. If this is the case, this would mean the conductance is lower than 10−7 G0. This is quite likely due to the cross-conjugated nature of the exTTF unit, and the known effect this has on conductance.
Despite the lack of a clear signal for the neutral molecule, we decided to proceed by trying to form the charge-transfer (CT) complex of compound 5 with 7,7,8,8-tetracyanoquinodimethane (TCNQ). When combining two equimolar solutions of the two compounds in DCM, we observed no color change, suggesting no complexation. The solvent was then changed to acetonitrile, a more polar solvent better able to stabilize complex formation, but this too did not give the anticipated strong color change. Only when the solution was heated at reflux for 3 h under ambient conditions and a large excess of TCNQ was added (5 equiv), a dark green color developed. UV–vis absorption spectroscopy confirmed the formation of a CT complex (Supporting Information File 1, Figure S1). To ensure the quantitative formation of the CT complex, a more straightforward method compatible with the conductance measurements was considered. Specifically, we switched from TCNQ to the stronger acceptor 2,3,5,6-tetrafluoro-7,7,8,8-tetracyanoquinodimethane (F4TCNQ). This produced a clear change of color as soon as the first drops of the acceptor were added to the solution of compound 5 in DCM, yielding a green solution. We monitored the formation of the CT complex by recording the UV–vis spectrum of the solution. As successive amounts of the acceptor are added to the solution of the donor, the peak at 448 nm of the neutral donor species decreases and new peaks appear between 600–900 nm that grow with the addition of more acceptor and also methanol (see Supporting Information File 1, Figure S2). Specifically, we observe three peaks at 689 nm, 764 nm and 867 nm, which can be assigned to the radical anion F4TCNQ species . For the break junction measurements, it was necessary to avoid having an excess of the acceptor because F4TCNQ forms molecular junctions itself due to its four terminal cyano groups (see Supporting Information File 1, Figure S7). Hence, we formed the charge-transfer (CT) complex by adding approximately 0.5 mL of a 10−4 M solution of F4TCNQ in DCM to a 1 mL 10−4 M solution of compound 5. In this way, we can avoid having significant amounts of free acceptor on the surface. This, however, means we do not exactly know the ratio of donor to acceptor in our case. It is known that the CT complex of the parent exTTF with TCNQ forms in a 1:4 ratio . However, as we have evidence of the formation of a radical anion of F4TCNQ, we believe a donor to acceptor ratio of 1:2 is more likely, for which there is a precedent in a substituted exTTF complex with TCNQ . We then allowed 24 h for the molecules to adsorb onto the gold in order to obtain as high a density of molecules as possible. This increases the possibility that a significant fraction of the solution is still the free donor species. Although there will still be some of the neutral donor species present on the surface, as we have shown that this molecule does not give a clear signal in break junction experiments, there will be no problem of signal overlap with the CT complex.
In contrast to the neutral form of compound 5, a significant percentage of conductance plateaus for the CT complex sample were observed. They were found fall into two main groups, labeled as high and low conductance. Figure 6a shows examples of individual G(z) traces displaying conductance plateaus. Firstly, we separate traces showing only tunneling (Figure 6b) from those containing plateaus (Figure 6c) using the following criteria: z interval (Δz) longer than 0.12 nm needed to observe a change of conductance Δlog(G/G0) = 0.16. We then further divided the traces into two more groups for those with plateaus above or below log(G/G0) = −3.8. Dividing the traces using a value slightly above or below this does not change the separation significantly as the difference between the two types of plateaus is clear. By fitting a Gaussian function to the histograms we find peak positions of log(G/G0) = −3.0 for the high group, and log(G/G0) = −4.7 for the low group (see Supporting Information File 1, Figure S6 for the 1D histograms). We also measured the junction break-off distance for the two groups (which we define as the separation required to move from log(G/G0) = −0.5 to either log(G/G0) = −4 for the high plateaus, or log(G/G0) = −6.1 for the low plateaus). We found mean values of 0.65 ± 0.25 nm and 1.14 ± 0.30 nm for the high and low groups, respectively (see Supporting Information File 1, Figure S3 for break-off histograms). We repeated the measurements using a freshly prepared CT complex and gold electrodes, and obtained a very similar result for the low conductance group (see Supporting Information File 1, Figures S4 and S5). In the repeated measurement, however, we did not observe a signal in the high conductance region, above 10−3 G0.
As can be seen from the separation of plateaus into high and low groups, the two types generally occur independently. The origin of the high conductance state is difficult to be totally sure about. The two groups may arise from independent chemical species, in which case it would be natural to label the high group as junctions for which transport takes place through a molecule of F4TCNQ. The length and conductance are similar to the control test we carried out on only this molecule (see Supporting Information File 1, Figure S7). It may also be possible that the high conductance state arises through contact to the central part of the CT complex and one of the terminal thiol groups. We observed a similar feature for the neutral molecule (see Figure 4b and Figure 4c) although for the neutral molecule, this signal was never very reproducible. The low group, on the other hand, fits well with conductance taking place across the whole molecule (i.e., thiol to thiol). We calculate an S–S distance of 2.4 nm for the CT complex (substituting the dihydroanthracene core of compound 5 for anthracene). We estimate the amount of gold retraction to be 0.5 nm, which then gives a real Au–Au separation of 1.64 ± 0.30 nm for the mean breaking distance of the low group, corresponding well with the length of the molecule. The fact that the molecule does not seem to fully stretch inside the junctions (unlike the OPE3-dithiol) may be due to the bulky nature of the complex, which contains many groups with the potential to interact with gold.
We cannot be totally sure of the species which gives rise to the low conductance state in so much as we do not know its exact charge and ratio of donor to acceptor molecules. It is conceivable that some charge is transferred to the electrodes upon molecular junction formation, as is common for molecular junctions . However, it is well known that the cationic (1+) state of exTTF is less thermodynamically stable than the dication (2+) , resulting in inverted oxidation potentials of the cation and dication. Additionally, it is known that the resonance between the neutral and dication is impossible due to the strong difference in geometries. Thus, it is more than plausible that the molecule retains its cationic state in the junction. We further point out that in several recent studies, a very similar conductance was found for a molecule that is structurally similar and contains only the anthracene central group. In a study by Hong et al., a conductance peak was found at log(G/G0) = −4.5, which is close to the low conductance group we observe . As it is known that the oxidation of exTTF compounds results in the planarization and aromatization of the central three rings to the anthracene unit, the similarity in conductance is perhaps indicative of these changes. We have also recently shown that adding substituents to the central phenyl rings of OPE wires has no noticeable influence on the molecular conductance (groups were tested ranging from electron-withdrawing fluorine to electron-donating methoxy groups) . It is also rational to suppose that the, albeit charged, dithiole side groups of our CT complex would not strongly modify the conductance of the otherwise neutral anthracene-containing wire. Finally, we point out that recent conductive AFM measurements carried out on self-assembled monolayers of anthracene and anthraquinone OPE derivatives showed a difference of approximately two orders of magnitude between their respective conductance values . This would fit with our hypothesis that the conductance of the neutral exTTF molecule, which has a similar bonding pattern to anthraquinone, is below our experimental resolution, based on the value we measure for the CT complex.
To gain insight into the conduction mechanism of compound 5, we performed theoretical calculations based on a combination of density functional theory (DFT) and Green’s functions techniques within the framework of the Landauer theory for coherent transport. The complete technical details are reported in Supporting Information File 1.
We first optimized the geometry of the molecule in the gas phase. The HOMO was found to be localized on the exTTF unit, which presented the expected butterfly shape, while the LUMO appeared delocalized over the whole molecule (Figure 7).
Subsequently, we constructed two kinds of metal–molecule–metal junctions, where the molecule is bound to a gold cluster in a top and hollow position, respectively (Figure 8). We then computed the zero-bias electron transmission following the procedure explained in Supporting Information File 1. The corresponding transmission curves are shown in Figure 9. Notice that the molecular HOMO–LUMO gap was corrected, following the procedure previously reported . The HOMO and all other occupied orbital energies were shifted by Σocc=−IP − εH +Δocc, while the LUMO and all other unoccupied orbital levels were shifted by Σvirt=−EA − εL +Δvirt. Here, Δocc (Δvirt) is the image charge correction for the occupied (unoccupied) states, εH (εL) is the Kohn–Sham energy of the gas phase HOMO (LUMO), and IP(EA) is the gas phase ionization potential (electron affinity). All quantities are reported in Table 1 for both binding geometries.
Table 1: Kohn–Sham HOMO and LUMO, ionization potential (IP), electron affinity (EA) and image charge correction for occupied Δocc and Δvirt unoccupied orbitals. All quantities are in eV.
The alignment of the Breit–Wigner resonances related to both the HOMO and LUMO (at approximately −1 and 2.7 eV from the Fermi level, respectively) do not show a strong dependence on the binding geometry (Figure 9). The electron transport is dominated by the HOMO, although interference features (resonance–antiresonance pairs) appear in the energy region close to the Fermi level. This is not surprising, given the spatial extension of the frontier orbitals, in particular of the HOMO, which is localized on the ex-TTF unit only. In fact, Fano resonances are known to arise when a “pendant” orbital, which is weakly coupled to the electrodes, is coupled to an orbital delocalized over the main axis of the molecular wire [38,39].
The computed conductance values, evaluated as the transmission at the Fermi level, are in the range 10−6 to 10−7 G0. These low conductance values are at the limit of what we can observe experimentally. This, therefore, would be consistent with the idea that the real conductance of the molecule is too low to be recorded in the experiments. Nevertheless, we cannot absolutely rule out the possibility that the molecule does not form molecular junctions in the experiments.
We have synthesized a molecular wire containing a π-extended tetrathiafulvalene (exTTF) group and studied its single-molecule electrical transport properties along with those of its charge-transfer complex with F4TCNQ. Within the accessible conductance range (10 to 10−7 G0) we did not observe a clear conductance signature of the neutral parent molecule. This alone could suggest either that its conductance is too low or that it does not form stable junctions. We did, however, find a clear conductance signature in the experiments carried out on the charge-transfer complex. As complexation with the acceptor oxidizes the molecule by removing two electrons from the exTTF group, thus converting it from a buckled and cross-conjugated group into a planar aromatic group, we predict the CT species to have a higher conductance than the neutral molecule. This, we believe, supports the idea that the conductance of the neutral molecule is very low, below our measurement sensitivity. This would make the conductance difference between the neutral and CT species at least two orders of magnitude. This can be considered as favorable for the use of single molecules as chemical sensors, in which analyte molecules may bind to a backbone to alter its conductance. Further combinations of donors and acceptors should be explored in order to evaluate this potential.
|Supporting Information File 1: Detailed experimental procedures for the synthesis and characterization of 5, break junction experiments and theoretical methods.|
|Format: PDF||Size: 625.2 KB||Download|
We thank Dr. Jose Manuel Santos Barahona for help with preparing the CT complex with TCNQ, and recording and interpreting the UV–vis spectra. Financial support by the European Commission (EC) FP7 ITN “MOLESCO” Project No. 606728, the European Research Council (ERC-436 2012 ADG_20120216-Chirallcarbon), the CAM (PHOTOCARBON project S2013/MIT-2841, NANOFRONTMAG-CM project S2013/MIT-2850 and MAD2D project S2013/MIT-3007), FP7-ENERGY-2012-1-2STAGE-number 309223 (PHOCS) and the Spanish MICINN/MINECO through the programs MAT2011-25046, MAT2014-57915-R and PRI-PIBUS-2011-1067 is acknowledged. MB was partly supported by a FY2012 (P12501) Postdoctoral Fellowship for Foreign Researchers from the Japan Society for Promotion of Science (JSPS) and by a JSPS KAKENHI, “Grant-in-Aid for JSPS Fellows”, grant no. 24·02501. YA is also thankful to another KAKENHI, “Grant-in-Aid for Scientific Research on Innovation Areas, Molecular Architectonics: Orchestration of Single Molecules for Novel Functions” (#25110009). LAZ was supported by the Spanish MICINN under Grant MAT2011-23627. FP acknowledges support by the Carl-Zeiss foundation and the Collaborative Research Center 767 “Controlled Nanosystems: Interaction and Interfacing to the Macroscale”.