Markus
Gantenbein‡
a,
Xiaohui
Li‡
b,
Sara
Sangtarash‡
*c,
Jie
Bai
b,
Gunnar
Olsen
a,
Afaf
Alqorashi
c,
Wenjing
Hong
*b,
Colin J.
Lambert
*c and
Martin R.
Bryce§
*a
aDepartment of Chemistry, Durham University, DH1 3LE, Durham, UK. E-mail: m.r.bryce@durham.ac.uk
bState Key Laboratory of Physical Chemistry of Solid Surfaces, NEL, iChEM, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China. E-mail: whong@xmu.edu.cn
cDepartment of Physics, Lancaster University, LA1 4YB, Lancaster, UK. E-mail: c.lambert@lancaster.ac.uk; s.sangtarash@lancaster.ac.uk
First published on 23rd October 2019
We report the synthesis of a series of oligophenylene-ethynylene (OPE) derivatives with biphenylene core units, designed to assess the effects of biphenylene antiaromaticity on charge transport in molecular junctions. Analogues with naphthalene, anthracene, fluorene and biphenyl cores are studied for comparison. The molecules are terminated with pyridyl or methylthio units. Single-molecule conductance data were obtained using the mechanically controllable break junction (MCBJ) technique. It is found that when electrons pass from one electrode to the other via a phenylene ring, the electrical conductance is almost independent of the nature of the pendant π-systems attached to the phenylene ring and is rather insensitive to antiaromaticity. When electrons pass through the cyclobutadiene core of the biphenylene unit, transport is sensitive to the presence of the relatively weak single bonds connecting the two phenylene rings of biphenylene, which arise from partial antiaromaticity within the cyclobutadiene core. This leads to a negligible difference in the molecular conductance compared to the fluorene or biphenyl analogues which have standard single bonds. This ability to tune the conductance of molecular cores has no analogue in junctions formed from artificial quantum dots and reflects the quantum nature of electron transport in molecular junctions, even at room temperature.
Important molecular parameters are the length, the conformation, the alignment of the molecular orbitals relative to the Fermi level of the metal leads, and the binding geometry at the molecule–metal contacts. Molecules with an oligo(phenyleneethynylene) (OPE) backbone and various anchor groups have been widely studied as single-molecule bridges between two metal electrodes. Para-Linked OPE-3 systems (3 refers to the number of phenylene rings in the backbone) are benchmark molecules in this context.15–24 OPEs are synthetically versatile and their conjugative and functional properties have been systematically tuned across many parameters.25–27 The alkyne bonds in para-linked OPEs serve two main purposes: (i) they ensure a length-persistent rigid-rod structure with no possibility of geometrical isomerization (unlike oligophenylenevinylenes) (OPVs) and (ii) they space the phenyl rings apart, which enables the rings to rotate freely and achieve coplanarity, thereby maximizing the frontier orbital overlap along the molecule. The OPE-type framework is therefore an ideal test-bed for probing the effects on single-molecule conductance of incorporating Hückel aromatic, non-aromatic or antiaromatic core units into a molecular backbone.
The role of aromatic and heteroaromatic rings in molecules wired into metal–single-molecule–metal assemblies has been well explored and in some cases the low-bias conductance is shown to be sensitive to the extent of aromaticity. Independent theoretical studies by Solomon and coworkers,28,29 and experimental STM-BJ studies by Venkataraman et al.30 concluded that increased aromaticity at the core of a molecular wire decreases the conductance of the molecular junction. For example, the experimental conductance of a series of molecules with amine anchors clearly followed the sequence 2,5-disubstituted-thiophene < -furan < -cyclopentadiene.30 In contrast, a study on a comparable series of 2,5-disubstituted-furan, -pyrrole, -thiophene and -cyclopentadiene derivatives with pyridyl anchors found no statistically significant dependence of the conductance on the aromaticity of the core.31 For molecules with tricyclic cores, it has been shown that increasing aromaticity at the core decreases the conductance for the para-linked molecules (dibenzothiophene < carbazole < dibenzofuran < fluorene) with pyridyl anchors,32 in agreement with Venkatarman et al.30 However, the sequence is different for the isomers where the tricyclic core is meta-linked into the backbone (dibenzothiophene ≈ dibenzofuran < carbazole ≈ fluorene).32
It is therefore of particular interest to explore the effect of incorporating 4n π-electron antiaromatic rings into the backbone of molecules that are wired into junctions. There are very few reported measurements on molecules of this type as they are usually demanding to synthesize and they are often unstable under ambient conditions. However, there are recent experimental precedents that antiaromaticity can be harnessed to enhance the conductance of single-molecule junctions. Yin et al. reported a single-molecule switch with a 9,9′–biindeno[2,1–b]thiophenylidene core that becomes antiaromatic with 6–4–6π-electrons upon electrochemical oxidation: a concomitant increase in conductance (on–off ratio of ca. 70) was observed.33 Fujii et al. reported that the conductance of an antiaromatic 16π-electron norcorrole-based nickel complex is more than one order of magnitude higher than that of its aromatic 18π-electron nickel-porphyrin based analogue.34
These two recent studies33,34 concern structurally rather complex molecules for which there are few appropriate model systems. Further work on new families of molecules is therefore timely in order to achieve a better understanding of the relationship between aromaticity, antiaromaticity and conductance in molecular junctions. For this study we focus on OPE derivatives with a biphenylene core unit. Biphenylene is a classic example of a stable molecule containing a 4-π-electron ring.35 The extent of resonance stabilization, electron delocalization and aromaticity of the peripheral benzene rings, and the related antiaromaticity of the central 4-π cyclobutadiene ring of biphenylene has been widely debated.36–40 The experimental and theoretical evidence shows that there is antiaromaticity in the central ring, but this antiaromaticity is partly alleviated by a degree of bond fixation in the benzene rings analogous to Kekulé-type structures. We are aware of only one report of biphenylene derivatives in a molecular junction. Biphenylene was 2,7-disubstituted with amine and cyclic thioether anchoring groups, and no significant increase in conductance was found in STM experiments, compared with the fluorene analogues.40 It is noteworthy that neutral biphenylene is isoelectronic with the highly-conductive 6–4–6 π-electron cation referred to above.33
We now report the synthesis, single-molecule conductance measurements and theoretical studies of eight OPE-based compounds shown in Fig. 1, with particular focus on the new biphenylene derivatives 1–4. Key molecular design features are as follows: (i) all the molecules have terminal pyridyl22,41 or methylthio22,42 anchor groups at both ends. These anchors are known to bind efficiently to gold22,41,42 and they have good chemical stability during the synthetic steps. They were chosen in preference to thiol anchors, because the thiol group would require additional protection/deprotection steps during the synthesis/assembly onto gold.43 (ii) The biphenylene cores are incorporated into the backbone with two different connectivities (1,4- or 2,7-difunctionalized).
![]() | ||
Fig. 1 (a) The structures of molecules 1–8 studied in this work. (b) Schematics of the MCBJ technique and the single-molecule junction of molecule 1. |
The 1,4–disubstitution (molecules 1 and 3) allows for a conduction pathway between the leads through a typical OPE-3 framework, with the cyclobutadiene ring as a pendant feature. On the other hand, the 2,7-disubstitution (molecules 2 and 4) dictates that the conduction pathway is through the entire biphenylene core. (iii) Molecules 5–8 are studied as model OPE analogues. At the outset the main question we sought to address was: Can biphenylene antiaromaticity lead to a measurable effect on charge transport in this series of molecules?
The route to the 2,7-difunctionalized biphenylene derivatives 2 and 4 starting from the readily-available commercial reagent 9 is shown in Scheme 2. Compound 16 was synthesized by a literature route in three steps from commercial 2,5-dibromonitrobenzene.45 Model compounds 4–8 were synthesized by analogous methods from the dihalogenated core units (see ESI†).
The mechanically controllable break junction (MCBJ) technique was used to characterize the single-molecule conductance of compounds 1–8, as shown in Fig. 1b. The experiments were carried out under ambient conditions by employing a homebuilt I–V converter with a sampling rate of 10 kHz.46 During the measurements, the breaking/connecting process of a notched gold wire was performed under the control of a combination of a stepping motor and a piezo stack. In this way, single-molecule junctions were formed by the repeatedly breaking/connecting of gold point contacts in solution (tetrahydrofuran:
1,3,5-trimethylbenzene = 1
:
4 v/v) with molecules at 0.1–0.4 mM concentration. Meanwhile the evolution of conductance characteristics was recorded at a fixed bias voltage of 100 mV, and then more than 1000 curves were used for further statistical analysis. More details are reported in our previous paper.47
Moreover for molecules 2 and 4, (with 2,7-connectivity) the single bond length of the cyclobutadiene core is longer than the single bonds within the phenylene rings (see Fig. S4 in the ESI†). Consequently, the electronic coupling in the former is weaker than in the latter. Therefore, transport is sensitive to the presence of the relatively weak single bonds connecting the two phenylene rings. This decreases the transmission within the gap and hence decreases the electrical conductance (see Fig. S6 in the ESI†).
To determine the effect on charge transport of a pendant cyclobutadiene unit on the OPE-3 system, molecules 5 and 7 comprising a similar core structure to 3 were investigated (Fig. 2d for 5 and 7). No significant difference in conductance values was observed among molecules 3, 5 and 7, which is evidence that the pendant cyclobutadiene of 3 has essentially no effect on charge transport. The low conductance state with broader peak width of 7 is attributed to either: (i) π-stacked dimer junctions through intermolecular interactions20,49–51 or (ii) the same single-molecule junction with different contact geometries.52 Molecules 6 and 8 which are analogues of the 2,7-disubstituted biphenylene 4, without the antiaromatic core, were also synthesized and studied, as shown in Fig. 2d. Biphenyl derivative 6 shows slightly lower conductance than that of 8, because of the larger torsion angle between the two phenyl rings of 6, compared to the planar fluorene unit of 8.31,53–55 Furthermore, the conductance of antiaromatic molecule 4 is comparable with that of 6 and 8, and no enhancement in charge transport due to the biphenylene unit in 4 is observed. This agrees with a previous comparison of conductance through a biphenylene and a fluorene core reported by Venkataraman et al.40
To reveal the evolution of the stretching process, 2D conductance-distance clouds were constructed by normalizing more than 1000 typical conductance traces to a relative zero point at 10−0.3G0 and plotted as intensity graphs54 as shown in Fig. 2e, f and S3.† The features at 1 G0 correspond to the construction of atomic gold–gold contacts just before the breaking process. The distinct high-density clouds between 10−4.0 and 10−6.0G0 are ascribed to the formation of single-molecule junctions. The relative stretching distance Δz histograms were constructed from 10−0.3G0 to the end of the conductance peaks in the 1D conductance histograms, which is one order of magnitude lower than the most probable conductance value. The most probable absolute stretching distance z* is obtained by adding the snap-back distance zcorr = 0.5 nm to the most probable relative stretching distance Δz*: z* = Δz* + zcorr.41 These values are in good agreement with the molecular length and the results are summarized in Table 1.
Compounds | Calculated conductance/log(G/G0)a | Measured conductance/log(G/G0)a | Measured length z*/nm |
---|---|---|---|
a Most probable conductance values and the error bars are based on the standard deviation in the Gaussian fitting of the 1D conductance histograms. | |||
1 | −4.05 | −4.6 ± 0.41 | 1.35 ± 0.17 |
2 | −4.90 | −6.1 ± 0.90 | 1.91 ± 0.23 |
3 | −3.47 | −4.4 ± 0.52 | 1.87 ± 0.22 |
4 | −3.65 | −4.9 ± 0.82 | 1.96 ± 0.32 |
5 | −3.35 | −4.2 ± 0.46 | 1.50 ± 0.23 |
6 | −3.60 | −4.8 ± 0.50 | 1.90 ± 0.27 |
7 | −3.40 | −4.1 ± 0.49 | 1.49 ± 0.17 |
−5.9 ± 0.77 | 1.90 ± 0.29 | ||
8 | −3.75 | −4.7 ± 0.44 | 2.34 ± 0.29 |
In order to model charge transport across these molecular junctions, we used scattering theory combined with density functional theory (DFT). The optimal geometry and ground state Hamiltonian were obtained using the SIESTA56 implementation of DFT and the transmission coefficients T(E) of electrons with energy E passing from one electrode to another through the molecules shown in Fig. S9† were calculated using the Gollum quantum transport code.57 Details of the computational methods are reported in the ESI.†
Fig. 3 shows the calculated T(E) of compounds 1–8, from which their electrical conductances are given by G/G0 ≈ T(EF).Since the Fermi energy of the electrodes EF relative to the frontier orbital energies is not usually predicted accurately by DFT, electrical conductances are plotted as functions of EF relative to the DFT-predicted value EDFTF. The highlighted area shows the Fermi energy at which the calculated conductances are in qualitative agreement with the experimental findings. In this region of EF, Fig. 3b and c show that molecules 3, 5 and 7, and molecules 4, 6 and 8, respectively, have similar conductances. Furthermore, Fig. 3a shows that the conductance of 3 is slightly higher than that of 4, in agreement with experiment, while the ratio of conductances of the pyridyl-terminated molecules 1 and 2 is significantly higher than that of the SMe-terminated molecules 3 and 4. Interestingly, this difference in the conductance ratios of 1:
2 compared with 3
:
4 would not be predicted by a simple tight binding model, in which all bonds are assigned the same hopping integral, because the cyclobutadiene unit enforces different bond lengths within the cores of these molecules. This aspect is discussed in detail in the ESI.†
![]() | ||
Fig. 3 DFT results of the transmission coefficients for (a) compounds 1–4, (b) 3, 5 and 7 and (c) 4, 6 and 8. The highlighted area shows the Fermi energy at which the calculated conductances are in qualitative agreement with the experimental findings. The relaxed structures are shown in Fig. S9 in the ESI.† |
For the future, it would be of interest to examine variants of these molecules with alternative connectivities, since both of the biphenylene cores considered here have odd–even connectivities (1,4- for 1 and 3, or 2,7- for 2 and 4). Fig. S6† shows a numbering system for the pz orbitals and transmission curves of the biphenylene core with a variety of connectivities. The calculations reveal that molecules with even–even (such as 2,8) or odd–odd (such as 1,7) connectivities exhibit a strong destructive interference feature within the HOMO–LUMO gap, independent of the degree of antiaromaticity of the cyclobutadiene core. However, these alternative connectivities of substituents onto the biphenylene core pose significant synthetic chemistry challenges. This ability to tune the conductance of molecular cores has no analogue in junctions formed from artificial quantum dots and reflects the quantum nature of electron transport in molecular junctions, even at room temperature.
Footnotes |
† Electronic supplementary information (ESI) available: Synthesis and characterization of compounds 1–8; methods and results for the single-molecule conductance measurements of compounds 1–8; computational methods for calculating the transmission coefficients. See DOI: 10.1039/c9nr05375a |
‡ M. G., X. L. and S. S. contributed equally to this work. |
§ M. R. B. coordinated the writing of the manuscript with contributions from all authors. All authors have given approval to the final version of the manuscript. |
This journal is © The Royal Society of Chemistry 2019 |