High conductance values in [pi]-folded molecular junctions

Abstract

Folding processes play a crucial role in the development of function in biomacromolecules. Recreating this feature on synthetic systems would not only allow understanding and reproducing biological functions but also developing new functions. This has inspired the development of conformationally ordered synthetic oligomers known as foldamers. Herein, a new family of foldamers, consisting of an increasing number of anthracene units that adopt a folded sigmoidal conformation by a combination of intramolecular hydrogen bonds and aromatic interactions, is reported. Such folding process opens up an efficient through-space charge transport channel across the interacting anthracene moieties. In fact, single-molecule conductance measurements carried out on this series of foldamers, using the scanning tunnelling microscopy-based break-junction technique, reveal exceptionally high conductance values in the order of 10^-1 G₀ and a low length decay constant of 0.02 Å^-1 that exceed the values observed in molecular junctions that make use of through-space charge transport pathways.

ARTICLE

Received 17 Aug 2016 | Accepted 9 Mar 2017 | Published 18 May 2017

Marco Carini1, Marta P. Ruiz2,3, Imanol Usabiaga4, Jos A. Fernndez4, Emilio J. Cocinero4, Manuel Melle-Franco5,6, Ismael Diez-Perez2,3 & Aurelio Mateo-Alonso1,7

Folding processes play a crucial role in the development of function in biomacromolecules. Recreating this feature on synthetic systems would not only allow understanding and reproducing biological functions but also developing new functions. This has inspired the development of conformationally ordered synthetic oligomers known as foldamers. Herein, a new family of foldamers, consisting of an increasing number of anthracene units that adopt a folded sigmoidal conformation by a combination of intramolecular hydrogen bonds and aromatic interactions, is reported. Such folding process opens up an efcient through-space charge transport channel across the interacting anthracene moieties. In fact, single-molecule conductance measurements carried out on this series of foldamers, using the scanning tunnelling microscopy-based break-junction technique, reveal exceptionally high conductance values in the order of 10 1 G0 and a low length decay constant of 0.02 1 that exceed the values observed in molecular junctions that make use of through-space charge transport pathways.

1 POLYMAT, University of the Basque Country UPV/EHU, Avenida de Tolosa 72, E-20018 Donostia-San Sebastian, Spain. 2 Department of Materials Science and Physical Chemistry, Institute of Theoretical and Computational Chemistry, University of Barcelona, Mart I Franqus 1, 08028 Barcelona, Spain. 3 Institute for Bioengineering of Catalonia (IBEC), Mart I Franqus 1, 08028 Barcelona, Spain. 4 Departamento de Qumica Fsica, Facultad de Ciencia y Tecnologa, Universidad del Pas Vasco (UPV/EHU), Apartado 644, E-48080 Bilbao, Spain. 5 CICECOAveiro Institute of Materials, Department of Chemistry, University of Aveiro, 3810-193 Aveiro, Portugal. 6 Centro ALGORITMI, 4710-057 Braga, Portugal. 7 Ikerbasque, Basque Foundation for Science, Bilbao, Spain. Correspondence and requests for materials should be addressed to M.M.-F. (email: mailto:[email protected]

Web End [email protected] ) or to I.D.-P.(email: mailto:[email protected]

Web End [email protected] ) or to A.M.-A. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:15195 | DOI: 10.1038/ncomms15195 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

DOI: 10.1038/ncomms15195 OPEN

High conductance values in p-folded molecular junctions

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15195

The high level of performance of proteins and nucleic acids in terms of function is not only the result of their chemical (primary) structure but also of the implicit

non-covalent intramolecular interactions between specic units in their sequence that promote folding into a well-dened (secondary and tertiary) three-dimensional structure. Recreating this feature on synthetic systems would not only allow reproducing biological functions but also developing new functions that align with our technological needs. This has inspired the development of conformationally ordered synthetic oligomers known as foldamers18, which are able to adopt folded conformations through non-covalent intramolecular interactions encoded in their primary structure. In the last years, there have been remarkable advances in the design of foldamers that self-assemble into well-dened structural motifs18, such as helices, pillars and sheets, among others. Many foldamers have been explored in diverse elds including protein mimicking915, molecular recognition1618, self-assembly1921, and photoinduced electron and energy transfer2225, with very promising results.

Charge transport in single molecules is of utmost importance in the development of molecular-scale electronic components. Among these, molecular wires2631 show a prominent position since they allow interconnecting different elements giving rise to complex molecular circuitry and to the transduction of biochemical events. A lot of effort2631 has been devoted to understanding through-bond charge transport across molecular wires with different lengths and conjugation to quantify their charge transport capabilities and to establish the charge transport mechanism. However, in higher-order systems, such as biological systems and organic solids, long-range electron transport does not only depend on the inherent properties of the individual molecules imposed by their chemical structure and p-conjugation but also on the spatial organization of such molecules32,33. In fact, a strong electronic coupling between the frontier orbitals of adjacent molecules, which is often associated with an excellent spatial overlap between their p-systems, is necessary to ensure efcient charge transport. This has motivated the study of synthetic molecular wires that impose through-space charge transport schemes either by a combined though-bond/ through-space channel3438 or by a through-space channel3942, and also of DNA4349 in which base pairs are stacked at virtually the graphitic interlayer distance.

Herein, we report a new family of foldamers consisting of an increasing number of anthracene units, which fold by a combination of intramolecular hydrogen bonds and aromatic interactions. Such folding process opens up an efcient through-space charge transport channel across the interacting anthracene moieties (Fig. 1). As a matter of fact, single-molecule conductance measurements carried out on this series of foldamers, using the scanning tunnelling microscopy-based break junction (STM-BJ) technique50,51, reveal exceptionally high conductance values (G) in the order of 10 1 G0 and a low length decay constant () of0.02 1. These values exceed the G and values observed in molecular junctions that make use of through-space charge transport pathways34,37,42 (G 10 4 to 10 2 G0, 0.21 to

1.16 1), including DNA43,44,48 (G 10 4 to 10 2 G0,

0.18 to 0.43 1).

ResultsDesign and synthesis. As our foldamers, we have designed and synthesized a series of oligomeric structures that combine 9,10-dimethyleneanthracene units stitched together by 2,6-pyridinedicarboxamides (Fig. 2). The series shows an increasing number of anthracene moieties: dimer-NH2; trimer-NH2; tetramer-

NH2; and pentamer-NH2. Anthracenes have been selected among

other acenes because of their well-balanced properties in terms of charge transport, stability and solubility52. The presence of the 2,6-pyridinedicarboxamides53 promotes a sigmoidal conformation because of the favourable intramolecular hydrogen-bonding interactions. This folded conformation is reinforced by aromatic interactions between the anthracene units, which are also favoured by the optimal parallel arrangement and the exibility provided by the methylene groups that connect the amide and the anthracene moieties. Triisopropylsilyl (TIPS) groups have been introduced in the pyridine units to ensure the solubility of the oligomers. Terminal amine groups have been selected as they have proven to be effective anchoring groups for attaining single-molecule junctions54.

The oligomers were rst synthesized as tert-butoxycarbonyl-protected amines (-NHBoc), which were deprotected afterwards in acidic media to make available the terminal amines (-NH2) as anchoring groups. The structure of both protected and deprotected foldamer series has been conrmed by nuclear magnetic resonance (NMR) and matrix-assisted laser desorption ionization (coupled to a time-of-ight (TOF) analyser) mass spectrometry. All details of the synthesis and characterization of the oligomers are given in the supporting information (Supplementary Note 1).

Determination of the folded structure. A combination of X-ray crystallography, NMR, ultravioletvisible absorption and photoluminescence (PL) spectroscopy, infrared ion dip spectroscopy and theoretical calculations provide strong evidence that the oligomers adopt a folded structure.

Crystals suitable for X-ray crystallography could be obtained for trimer-NHBoc from dimethylsulfoxide (DMSO) solutions. The X-ray diffraction analysis conrms a sigmoidal and folded

Folding

Anchoring

Charge transport channel

Figure 1 | p-Folded molecular junctions. Schematic representation of the folding and anchoring processes needed to obtain p-folded molecular junctions from a representative member of the foldamer family studied in this work.

2 NATURE COMMUNICATIONS | 8:15195 | DOI: 10.1038/ncomms15195 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15195 ARTICLE

Si

O O

N

ZHN

NHZ

H H

N

N

N

H

H

H

N N

N

H H

N N

O O O

NHZ

N

O

Si

Si

Si

Si Si

Z = Boc;

Z = H; Dimer-NH2

Dimer-NHBoc TFA Z = Boc; Trimer-NH2

Z = H;

Trimer-NHBoc TFA

O

O O

N H H

O O

N

N

N

ZHN ZHN

ZHN

H

N

NHZ

N H H

N

N

N

H H

N

H H

N

N N N

NHZ

O O

O

O O

H H

N N N

O O

H H

N N N

O O

Si Si Si Si

Pentamer-NHBoc TFA

Figure 2 | Chemical structures of the foldamers. Structures of the complete -NHBoc and -NH2 foldamer series studied in this work: dimer-NHBoc, dimer-NH2, trimer-NHBoc, trimer-NH2, tetramer-NHBoc, tetramer-NH2, pentamer-NHBoc, and pentamer-NH2. Arrows indicate the transformation of the -NHBoc series into the -NH2 series.

Z = Boc;

Z = H; Tetramer-NH2

Tetramer-NHBoc TFA Z = Boc; Pentamer-NH2

Z = H;

structure of trimer-NHBoc and also provides an insight into the interactions between the anthracene moieties in the solid state (Fig. 3a). Remarkably, even considering the interference of DMSO with the folding processthe X-ray structure reveals the presence of two DMSO molecules bound to the 2,6-pyridinedicarboxamides by hydrogen bondsthe high predisposition of the anthracene moieties to interact between each other is evidenced. In fact, the central anthracene adopts a CHp conformation by ipping over the DMSO molecules to stick to the two external anthracene moieties.

To study the folding process in solution, we carry out all the experiments with the NHBoc foldamer series as it represents a closer model to the break-junction scenario, in which the free amines are anchored to the gold tips. 1,1,2,2-Tetrachloroethane (TCE) is selected as a solvent to carry out the studies because the oligomers are highly soluble in it; it has a high boiling point (146 C) optimal for the break-junction studies; and it has a low hydrogen bond basicity (SbH 0.08)55, which favours

the formation of intramolecular hydrogen bonds that hold together the 2,6-pyridinedicarboxamide backbone in a sigmoidal conguration. The NMR anthracene aromatic signals of dimer-NHBoc and trimer-NHBoc in TCE-d2 at 25 C appear shifted upeld and broadened in comparison to the Boc-protected 9,10-bis(aminomethyl)anthracene (A-NHBoc) that is used as a reference (Fig. 3b,c). This is consistent with the close proximity between anthracene moieties, which results in the observed shielded signals as an effect of anisotropy. Moreover, a more substantial broadening is clearly observed in the NMR spectra of tetramer-NHBoc and pentamer-NHBoc that leads to an almost complete disappearance of the aromatic

signals. We rst ascribed this effect to the aggregation of the tetramer-NHBoc and pentamer-NHBoc in solution at NMR concentrations (10 3 to 10 4 M). However, this contradicts the high solubility of the whole NHBoc foldamer series in TCE.

Two-dimensional diffusion-ordered spectroscopy (DOSY) NMR studies at 25 C reveal that the broadening of the aromatic signals is indeed not related to aggregation in solution, as the signals of the oligomers diffuse at the same rate, which means that the oligomers are monodisperse species in solution (Fig. 3d). Therefore, the complex NMR observed for tetramer-NHBoc and pentamer-NHBoc indicates the coexistence of multiple conformers at room temperature in a complex dynamic equilibrium. This is supported by variable temperature NMR (VT-NMR) experiments carried out from 40 to 140 C that show how the complex ne structure of the pentamer-NHBoc coalesces at increasing temperatures into a simpler NMR spectrum that can be easily assigned to the structure (Fig. 3e), as the high temperatures disrupt the intramolecular interactions, unfolding the oligomers.

The ultravioletvisible electronic spectra in TCE showed the typical features of anthracene derivatives with four vibronic bands on the main electronic absorption (Fig. 3f). An electronic interaction between anthracene moieties is clearly observed, as the absorption bands appear increasingly bathochromically shifted with an increasing number of anthracene units in the oligomers. When dimer-NHBoc was compared with reference A-NHBoc, a bathochromic shift (Dl) of 1.0 nm was observed.

The electronic interactions between anthracenes become more evident when increasing the number of anthracenes in the oligomers as the bathochromic shifts become increasingly larger

NATURE COMMUNICATIONS | 8:15195 | DOI: 10.1038/ncomms15195 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15195

a

b c

A-NHBoc

Dimer-NHBoc

Trimer-NHBoc

Tetramer-NHBoc

Pentamer-NHBoc

Pentamer-NHBoc

NHBoc

BocHN

NHBoc

A-NHBoc

9

140 C 130 C 120 C 110 C 100 C

90 C 80 C 70 C 60 C 50 C 40 C 30 C 20 C 10 C

0 C 10 C 20 C 30 C 40 C

8 7 6 5

d e

1E-4

1E-5

1E-6

1E-7

Pentamer-NHBoc

A-NHBoc Dimer-NHBoc Trimer-NHBoc Tetramer-NHBoc Pentamer-NHBoc

Diffusion units

9

8

7 6 5

9 8 7 6 5

f g

1.2

A-NHBoc Dimer-NHBoc Trimer-NHBoc Tetramer-NHBoc Pentamer-NHBoc

Normalized absorbance

1

0.8

0.6

0.4

0.2

0

Normalized photoluminescence

1

0-0

0-1 0-2 0-3

0.8

0.6

0.4

0.2

0

300 320 340 360Wavelength (nm)

(Experimental)

(Simulated)

380 400 420

400 450 500 550 600 650

Wavelength (nm)

h i

(Experimental)

(Simulated)

Dimer-NHBoc

Trimer-NHBoc

2,800 3,000 3,200 3,400Wavenumber (cm1)

3,600 3,800

2,800 3,000 3,200 3,400 3,600 3,800

Wavenumber (cm1)

Figure 3 | Structural characterization of the foldamers. (a) Solid-state structure of trimer-NHBoc with two DMSO molecules (highlighted in yellow). (b) Chemical structure of A-NHBoc. (c) NMR spectra of the foldamer NHBoc series at 25 C. (d) Two-dimensional DOSY NMR of pentamer-NHBoc at 25 C. (e) VT-NMR of pentamer-NHBoc. (f) Absorption and (g) PL spectra of the NHBoc foldamer series at 25 C. Infrared ion dip spectrumof (h) dimer-NHBoc and (i) trimer-NHBoc. The triangles in be indicate the signals that correspond to anthracene protons.

in trimer-NHBoc (Dl 1.5 nm), tetramer-NHBoc (Dl 2.0 nm)

and pentamer-NHBoc (Dl 5.5 nm), when compared with

A-NHBoc, respectively. Almost invariable changes were observed in the absorption spectra on dilution (from 10 4 to 10 6 M), which is consistent with the intramolecular nature of the interactions between anthracene units.

The intramolecular aromatic interactions between anthracenes are also evidenced in the ultravioletvisible PL studies. The PL spectra of the -NHBoc foldamer series show the typical emission features of anthracene derivatives with four vibronic bands (Fig. 3g). Overall, the observed PL bands are bathochromically

shifted as the number of anthracene moieties in the oligomers increases, in line with the absorption measurements. When dimer-NHBoc was compared with A-NHBoc, using the 0-1 transition as a reference, a bathochromic shift (Dl 1.0 nm)

is clearly observed. The bathochromic shifts become increasingly larger in trimer-NHBoc (Dl 1.5 nm), tetramer-NHBoc

(Dl 2 nm) and pentamer-NHBoc (Dl 6.0 nm). These energy

changes are coupled to a change on the relative intensity of the 0-2 and 0-3 transitions that increase together with the number of anthracene moieties in the oligomer. Such increase of the relative intensity is due to an underlying transition that corresponds to

4 NATURE COMMUNICATIONS | 8:15195 | DOI: 10.1038/ncomms15195 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15195 ARTICLE

a

the anthracene excimer (B500 nm)56. The presence of the anthracene excimer in the PL spectra evidences the face-to-face overlap between anthracenes in solution. Most importantly, almost invariable changes were observed in the PL spectra on dilution (from 10 5 to 10 7 M), which also supports the intramolecular nature of the interactions between anthracenes.

To obtain a much clearer view of the folding process, conformation-specic vibrational spectroscopic studies in the gas phase were carried out as they allow elucidating the intrinsic structures of moderately large molecules (typically o500600 Da)

without the interferences of packing or solvation. Remarkably, dimer-NHBoc (984.33 Da) and trimer-NHBoc (1,532.10 Da) could be studied even if they roughly duplicate and triplicate, respectively, the typical molecular weight for this kind of measurements57,58. The full experimental description is given in Supplementary Information. In brief, an infrared laser pulse (1,064 nm) was used to transfer the intact foldamer molecules from solid state to gas phase. Then, the target vapourized sample was diluted and expanded in an argon jet and skimmed to generate a cold collimated beam. In the extraction and acceleration regions of the linear TOF, an ultraviolet dye laser ionized the foldamers through a resonant two-photon ultraviolet excitation (B34,000 cm 1). Conformation-selective spectra could be obtained from a third tunable optical parametric oscillator/optical parametric amplier infrared laser through infrared ion dip spectroscopy. The conformationally specic vibrational spectra could be interrogated in the CH and NH stretch regions, which are extremely sensitive to the local hydrogen-bonded conformational environment. The infrared experimental spectra of dimer-NHBoc and trimer NH-Boc are very similar (Fig. 3h,i) with a set of bands in the 2,8503,100 cm 1 region, which correspond to the CH stretch vibrational modes. With regard to NH modes, the experimental band about 3,430 cm 1 corresponds to the symmetric and antisymmetric coupled modes of the amide nitrogens. The presence of a single red-shifted band points out that the amide nitrogens are equivalent and that they are weakly hydrogen bonded to the nitrogen atoms of pyridine rings. A second set of NH bands should be observed, which would correspond to the carbamate nitrogens. However, these modes are not observed in the experimental spectra. This is not strange since these atoms are not forming hydrogen bonds and hence, these bands should be weaker. The spectra of several minimized conformers of dimer-NHBoc and trimer-NHBoc were calculated and compared with the experimental spectra to assign the intrinsic folded conformations. To ensure that the assignment is correct, the spectra were calculated from the minimized structure dimer-NHBoc and trimer-NHBoc exhibiting TIPS groups. The good agreement between experimental and calculated spectra (Fig. 3h,i) illustrates that the foldamers adopt a pp conformation between the anthracene moieties (see below). The simulated spectrum for the alternative CHp conformation for dimer-

NHBoc does not t with the experimental data conrming univocally that the observed conformation is pp (Supplementary

Fig. 1). In the case of trimer-NHBoc the simulations of the pp and the CHp conformers are similar and both could account for the experimental spectrum (Supplementary Fig. 2). However, the higher correlation of the simulated pp conformer with the experimental spectrum, the clear observation of the pp conformation for dimer-NHBoc and the larger stability found for the pp conformers from modelling (see below) indicate a preferential pp conformation for trimer-NHBoc, as well.

Overall, from the combination of all different characterization techniques, we can safely conclude that the foldamers adopt a folded conformation in the solid state, in solution and in the

gas phase. X-ray crystallography of trimer-NHBoc shows a folded structure in the solid state even in the presence of two crystallization DMSO molecules that impose a CHp conformation (Fig. 3a). NMR, ultravioletvisible absorption and PL measurements in solution (Fig. 3cg) show the close proximity of the anthracene moieties (shielding of the NMR anthracene signals) and an intramolecular electronic interaction among them (batochromically shifted absorption and emission features) that gives rise to a complex intramolecular dynamic equilibrium between different pp and CHp conformers (a single diffusing DOSY signal for the foldamers and also denaturation of the foldamer at high temperature observed in VT-NMR, with a preferential pp conformation as the number of anthracenes increases (anthracene excimer emission evidencing a face-to-face overlap between anthracenes). Conformationspecic vibrational spectroscopic studies conrm a more favourable pp conformation for dimer-NHBoc and trimer-NHBoc in the gas phase (Fig. 3h,i).

Theoretical calculations were carried out to provide a detailed picture of the fundamental interactions determining the energetics and structure of the folded molecules (Fig. 4a). The TIPS groups have a marginal effect on the geometrical and electronic structure and were omitted. The effect of the (implicit) solvent was modelled and was found to have a negligible effect on the geometries of dimer-NH2 and therefore was not used in

CHTT

TTTT

b

LUMO

HOMO

Figure 4 | Computed structure of the foldamers. (a) Representative CHp and pp structures for the different foldamers. (b) HOMO and LUMO orbitals for CHp and pp conformers of pentamer-NH2.

NATURE COMMUNICATIONS | 8:15195 | DOI: 10.1038/ncomms15195 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15195

the geometry optimizations. On the basis of the spectroscopic and X-ray data described above two different types of folded conformations were explored: CHp and pp. On one hand, the CHp conformations were based on the crystal structure of trimer-NHBoc using partial optimizations at the B3LYP/6-31 g(d,p) level in which only the amines and hydrogens were allowed to relax, as full optimizations lead to pp structures.

Trimer-NHBoc crystallizes in a CHp conformation yet, when the co-crystallizing solvent molecules were removed and the structure was optimized in vacuum, it changed into a pp structure, which indicates that interconversion between

CHp and pp folded structures is energetically feasible. On the other hand, classical molecular dynamics simulations were run followed by optimizations at the BLYP/def2-TZP-D3BJ level to obtain the pp conformers studied (see Supplementary

Information for full details). Overall, density functional theory molecular mechanics and semi-empirical models are found to favour compact pp conformations, in line with our experimental observations. The rst CHp stabilized structure of dimer-NH2,

which also shows the highest CHp character of the whole foldamer series, is predicted to be 2.1 kcal mol 1 above the lowest lying pp conformer. The distance between stacked anthracenes for the pp dimer-NH2 is 3.26 . The structures of trimer-NH2, tetramer-NH2 and pentamer-NH2 are substantially more complex. As a matter of fact, the computed structures conrm that the foldamers are held together by a continuous interplay of aromatic interactions between anthracene units, consistent with the VT-NMR and DOSY experiments that conrm a complex conformational equilibrium.

We computed the electronic structure of CHp and pp conformations at the B3LYP/6-311 g(2d,p) level surrounded by

implicit solvent (Fig. 4b). The electronic structure of the frontier orbitals for the four foldamers is qualitatively similar for the pp and the CHp series. The HOMO orbitals for all foldamers have consistently large densities on anthracene moieties, while the LUMO orbitals have often, but not exclusively, larger densities on the pyridine-2,6-dicarboxamide residues (Fig. 4b and Supplementary Figs 36). In addition, in all studied cases, the terminal amines are electronically coupled to the terminal anthracenes in the HOMO orbitals. In the case of the CHp conformers, the anthracene units are electronically decoupled at the HOMO level. In contrast, the pp conformers, due to the stacking conformations, show HOMO orbitals with densities that spread over more than one anthracene moiety. Furthermore, in all foldamers, the CHp and pp conformations show nearly degenerate occupied orbitals for each anthracene unit with energies from 5.3 to 5.6 eV (Supplementary Table 1). In particular, the HOMO degeneracy of the pp conformers extends to the HOMO1 (dimer-NH2, trimer-NH2, tetramer-NH2 and pentamer-NH2), HOMO2 (trimer-NH2, tetramer-NH2 and pentamer-NH2), HOMO3 (tetramer-NH2 and pentamer-NH2) and HOMO4 (pentamer-NH2 only), and in most cases with densities that spread over more than one anthracene moiety (Supplementary Figs 57). All foldamers also share very similar ionization potentials ranging from 5.3 to 5.4 eV. A comparison of the HOMO and HOMO1 for the CHp and pp conformers in the simplest foldamer, dimer-NH2, illustrates the differences between conformers (Supplementary Fig. 7). For the pp conformer, there is an electronic density path linking the terminal amines and both anthracene moieties in the HOMO and HOMO1. Furthermore, for this conformer both molecular orbitals present very similar electron densities and energies (110 meV). In contrast, for the CHp conformer, the HOMO and the HOMO1 orbitals show localized electron densities on different anthracenes and a slightly larger energy difference (190 meV). Overall, the pp and CHp conformers for the whole foldamer series exhibit degenerate

occupied frontier orbitals with large electron density on the anthracene moieties and with ionization potentials which are closely aligned to the Fermi energy of gold ( 5.5 eV). This

provides a viable pathway for charge transport across the molecule through the anthracene moieties.

Single-molecule charge transport studies. Single-molecule charge transport characterization of each oligomer was carried out in single-molecule junctions using the STM-BJ technique50,51. The experiments were conducted in TCE at very low concentrations (o10 9 M) to ensure an efcient anchoring of individual molecules between the electrodes and to avoid any aggregation in solution. Briey, the STM-BJ is based on driving a Au STM tip in and out of contact to/from a Au(111) substrate functionalized with the target foldamer. During the contact process, individual oligomers can spontaneously bridge between both biased (few mV) electrodes via the two -NH2 terminal groups (Fig. 1, bottom). The current is then recorded for each pulling stage in the form of current versus time/displacement (Fig. 5a), and the traces displaying molecular plateau features below the quantum conductance (G0 2e2h 1) were used to

build the conductance histograms (Fig. 5b). The observed maxima in the histograms represent the most probable conductance values for the formed single-molecule contact (see Supplementary Information for details). The most intense (higher counting) conductance peak, which corresponds to the most probable molecular junction conguration, shows very large conductance values in the range of 10 1 G0 for all foldamers:0.215 G0 for dimer-NH2; 0.195 G0 for trimer-NH2; 0.178 G0 for tetramer-NH2; and 0.168 G0 for pentamer-NH2. Such high conductance values decrease exponentially with an increasing number of anthracenes in the molecular backbone. By tting the conductance data to a simple tunnelling transport regime, G Ae L, and using the distance between p-stacked

anthracenes estimated from the simulations (3.26 ), a very low length decay constant 0.02 1 is obtained (Fig. 4c).

To experimentally conrm that charge transport takes place preferentially through the interacting anthracene moieties, a truncated trimer (#-trimer-NH2) was synthesized, in which the central anthracene ring has been substituted by a non-aromatic cyclohexane ring (Fig. 5d,e). Conductance histograms for #-trimer-NH2 showed that the highest conductance peak(0.017 G0) is one order of magnitude smaller than the corresponding conductance observed for the homologous trimer-NH2 (0.195 G0). This substantial conductance drop can be ascribed to the interrupted conduction channel because of the lack of an anthracene moiety in #-trimer-NH2 that interferes with charge transport. Consequently, the high conductance values observed are undoubtedly owed to the presence of a charge transport channel across the interacting anthracene units.

To validate the above proposed tunnelling charge transport mechanism in the studied single-molecule junctions, we have measured temperature-dependent single-molecule conductance for the longer pentamer-NH2 moiety (Fig. 5f,g). To avoid any possible unfolding, we covered a temperature range below room temperature from 6 to 25 C. The invariance of the single-molecule conductance versus temperature around the working room temperature for the longer backbone (Fig. 5g) supports a tunnelling transport regime throughout all the conductance series in Fig. 5b,d. It is worth noticing in Fig. 5g that while the conductance dispersion (represented as error bars) is fairly similar within the studied temperature range, the single-molecule experimental yield (dened as the percentage of current traces displaying single-molecule events or plateaus over the total collected traces) drops down by a factor of four from the

6 NATURE COMMUNICATIONS | 8:15195 | DOI: 10.1038/ncomms15195 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15195 ARTICLE

a

G(G 0)

0.2

0.1

Dimer-NH2

Trimer-NH2

Tetramer-NH2

Pentamer-NH2

0.2 nm 0.4 nm

0.2 nm 0.2 nm

Electrode distance (nm)

b c

Dimer-NH2 Trimer-NH2

#-Trimer-NH2

0.215 G0

0.195 G00.178 G00.168 G0

Normalized counts (a.u.) Normalized counts (a.u.) Normalized counts (a.u.)

[afii9826] = 0.02 1

0.4

0.1

Tetramer-NH2 Pentamer-NH2

G(G0 )

0.1

G(G0)

6.0 12.0

d e

Si

O O

O

N

0.017 G0

H N

N

H H

N

N

H

H

N N

O

NH

#-Trimer-NH2

0.01

0.1 0.0034 0.0035

1/T (K1)

Si

G(G0)

G(G0)

f g

6 C 10 C 15 C 20 C RT

0.1

G(G0 )

0.0036

Figure 5 | Single-molecule conductance measurements. (a) Individual pulling traces and (b) conductance histograms of the different foldamers. All conductance values were extracted from Gaussian ts. (c) Calculated value for the pp foldamers. Error bars capture the variability in peak position.

(d) Conductance histograms for #-trimer-NH2. (e) Structure of #-trimer-NH2. (f) Conductance histograms for the pentamer-NH2 at variable temperatures. (g) Arrhenius plot of the single-molecule conductance for the pentamer-NH2. Error bars represent the s.d. of the experimental conductance values extracted from the full width at half maximum of the Gaussian ts. Applied voltage bias were 5 mV for the trimer-NH2pentamer-NH2 series and 50 mV for dimer-NH2.

lowest to the highest working temperature, which results in the observed counts decrease in Fig. 5f histogram peaks from top to bottom. This behaviour has been observed before in single-molecule transport studies59 and is ascribed to thermal instabilities at the molecule/electrode contact as the temperature is raised.

Overall, the conductance measurements carried out point towards a coherent tunnelling charge transport process across the anthracene units as supported by the presence of a charge transport channel across the interacting anthracene unitsthis is consistent with the calculated electronic structure and evidenced by the conductance drop observed when

comparing trimer-NH2 and #-trimer-NH2 (Fig. 5b,d)and by the invariance of the single-molecule conductance versus temperature (Fig. 5g), which supports a tunnelling transport regime.

The moleculeelectrodes electronic coupling strength (G) has been estimated from a single-channel resonant tunnelling model60,61 yielding values between 50 and 60 meV from pentamer-NH2 to dimer-NH2 (Supplementary Note 2). Even if such values support the observed coherent tunnelling behaviour60, they cannot be used to conrm whether charge transport is resonant, but nevertheless, they point out to a high coupling scenario62.

NATURE COMMUNICATIONS | 8:15195 | DOI: 10.1038/ncomms15195 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15195

DiscussionWe have reported a series of synthetic oligomers consisting of an increasing number of anthracene moieties that fold by a combination of hydrogen bonds and aromatic interactions as demonstrated by a series of studies in the solid, liquid and gas phases, and by theoretical calculations. Remarkably, such folding process opens up an efcient through-space charge transport channel across the interacting anthracene moieties that has been conrmed both theoretically and experimentally. First, the computed electronic structure of the frontier orbitals of the foldamer series illustrates the degeneracy of anthracene occupied frontier orbitals, which provides a viable pathway for charge transport across the anthracene moieties. Second, charge transport studies at the single-molecule level show exceptionally high conductance values in the range of 10 1 G0 for all foldamers and a low length decay constant of 0.02 1. Third, comparative conductance studies carried out with #-trimer-NH2, in which the central anthracene ring has been substituted by a non-aromatic cyclohexane ring, reveal conductance values one order of magnitude smaller than the homologous trimer-NH2, conrming that charge transport takes place preferentially through the channel constituted by the interacting anthracene moieties. While a complete picture of the transport mechanism at play is still to be determined, the invariance of the single-molecule conductance versus temperature supports a tunnelling transport regime. Overall, this work illustrates that carefully designed foldamers provide a new and efcient means for charge transport at the single-molecule level and also new application perspectives for foldamers.

Methods

Synthesis and structural characterization. All commercially available reagents and solvents were used without further purication. Anhydrous tetrahydrofuran (THF) and toluene were dried through a SPS solvent purication system. All reactions of the compounds were carried out under a nitrogen atmosphere and oven-dried glassware. Column chromatography was carried out using Silica gel 60 (4060 mm; 230400 mesh) from VWR. Thin layer chromatography was performed to follow the reaction process by using sheets (20 20) of aluminium pre-coated

with Silica gel 60 F254 from Merck. Ultraviolet-active compounds were detected with an ultraviolet lamp from CAMAG at wavelength l 254 or 366 nm. 1H-NMR

and 13C-NMR spectra were recorded on Bruker Avance 400 or 500 spectrometer at 298 K using partially deuterated solvents as internal references. Chemical shifts are reported in (d) parts per million and referred to the residual solvent peak. DOSY and VT-NMR measurements were performed at the NMR General Services of the University of the Basque Country by Dr Jos Ignacio Miranda. Matrix-assisted laser desorption ionization (coupled to a TOF analyser) experiments were recorded on Bruker REFLEX spectrometer at POLYMAT by Dr Antonio Veloso. X-ray diffraction measurements were performed at the X-ray diffraction unit of the General Services of the University of the Basque Country by Dr Leire San Felices. Intensity data were collected on an Agilent Technologies Super-Nova diffractometer, which was equipped with monochromated Cu ka radiation(l 1.54184 ) and Atlas charge-coupled device detector. Measurement was

carried out at 150.00(10) K with the help of an Oxford Cryostream 700 PLUS temperature device. Data frames were processed (united cell determination, analytical absorption correction with face indexing, intensity data integration and correction for Lorentz and polarization effects) using the Crysalis software package. The structure was solved using Olex2 and rened by full-matrix least-squares with SHELXL-97. Final geometrical calculations were carried out with Mercury and PLATON as integrated in WinGX.

Foldamer handling. The -NHBoc and -NH2 foldamer series are sensitive to oxygen and light and decompose in solution to the corresponding anthracene endoperoxides in about 24 h. All anthracene-containing compounds were characterized from fresh samples in a light-free environment. In the case the experiments required more than 2 h they were carried out in a nitrogen- and light-free environment (either under an nitrogen stream, using airtight glassware or handled in a glovebox).

Optical characterization. Absorption and PL spectra were recorded ona Perkin-Elmer Lambda 950 spectrometer and a LS55 Perkin-Elmer Fluorescence spectrometer, respectively, at 298 K.

Infrared ion dip spectroscopy. Approximately 40 mg of foldamer (dimer-NHBoc or trimer-NHBoc) were used in the experiments. The sample was mixed with nanotubes and deposited on a graphite cylindrical rod (15 mm long and 4 mm diameter). A stepper motor rotated and translated the sample. So, the infrared pulses (1,064 nm, B8 ns) from a Q-switched Nd:YAG laser found fresh sample.

A full description of technique is described here63. The foldamers exhibit broad absorption electronic spectra due to several chromophore anthracene groups. The infrared spectra were recorded at 34,343 and 34,013 cm 1 for dimer-NHBoc and trimer-NHBoc, respectively. Finally, the experimental infrared spectra were compared with theoretical spectra. The frequencies of the NH and CH stretch modes, expressed in wavenumbers, were corrected for anharmonicity using the multipliers, 0.995 and 0.98, respectively.

Computer modelling. We studied different conformers for the dimer-NH2 and trimer-NH2 foldamers. The TIPS groups have a marginal effect on the geometric and electronic structure of these systems and were omitted for computational efciency. To generate the starting structures, we rst performed molecular dynamics simulations at 298 K with the molecular mechanics MMFF94 force eld as implemented in the software Tinker. Collisions with solvent molecules were considered implicitly through the use of stochastic dynamics with a friction coefcient of 2 ps to improve conformational sampling. Simulations were run for 10 ns and snapshots of the trajectory were selected at 50 ps intervals yielding200 different conformations, which were optimized at MMFF94 level and at the PM6-D3H4 level64. The so-obtained geometries were ltered by an algorithm explicitly considering topological symmetry with a root mean squared deviation cut-off of 0.5 (ref. 65). The unique conformations were subsequently minimized at the BLYP/def2-SVP-D3BJ (also used for infrared spectra simulations) and BLYP/def2-TZVP-D3BJ (geometries and total energies) levels66 within the RI approximation with the programme Orca 3.03 (ref. 67). The effect of the (implicit) solvent was modelled with the programme Gaussian 09 (ref. 68) and found to have a negligible effect on the geometries of the dimer-NH2 and was not used in the geometry optimizations. The aforementioned procedure produces prevalently conformers with pp stacked anthracenes. The lowest energy minimum for the dimer-NH2 and the minima presenting the largest contact between anthracene moieties in the trimer were selected as representative for the pp conformations. The pentamer-NH2 pp conformation was produced by fusing two trimer-NH2 molecules followed by optimization, while the pp tetramer-NH2 was obtained by cutting out one monomer from the pentamer-NH2 followed by optimization. CHp conformers were derived from the trimer-NHBoc crystalline structure by, rst, substituting the NHBoc groups by amine groups and, second, using partial optimizations at the B3LYP/6-31g(d,p) levels where only the amines and hydrogen atoms were allowed to relax. This was necessary as full optimizations favour pp stacked structures. The CHp pentamer-NH2 was built by fusing two CHp trimer-NH2 molecules, but in this case only the central CH2antraceneCH2 unit was optimized. The CHp tetramer-NH2 was obtained cutting out one monomer from the pentamer. The electronic structure for neutral and charged foldamers was computed at the B3LYP/6-311 g(2d,p) level, as with this Hamiltonian HOMO

and LUMOs eigenvalues are similar to the adiabatic ionization potential and electron afnity of conjugated molecules in low-polarity solvents69. The adiabatic ionization potential and electron afnity were explicitly computed for all foldamers and match the HOMO and LUMOs eigenvalues for these systems.

STM break junction. Foldamer-NH2 TFA salt was dissolved in MeOH (saturated) and the free amine was precipitated by addition of triethylamine (TEA), the solid was collected by centrifugation and was washed with MeOH/TEA 20:1 and twice with MeOH and nally vacuum-dried until constant weight. All the conductance measurements were carried out in the dark in anaerobic conditions with a mechanically and electronically isolated PicoSPM I microscope head controlled by a Picoscan 2500 electronics (all from Agilent) and using a homemade PTFE STM cell. Data captures were acquired using a NI-DAQmx/BNC-2110 National Instruments (LabVIEW data acquisition System) and analysed with LabVIEW code. The procedure of a typical break-junction experiment is based on bringing the STM tip to tunnelling distance over a at clean Au (111) surface area as a rst step. The STM feedback is then turned off and the tip is driven into and out of contact with the substrate at a speed of 12 V s 1. These two-point feedback loop is used to collect thousands of current decays (5,0006,000). Single-molecule conductance (G) was determined using the expression G Istep/UBIAS, where I is the current and

U is the voltage difference between the two junction electrodes. The current decays are accumulated to semi-logarithmic conductance histograms. The observed plateaus in the individual current decays result in the observed peaks in the conductance histograms and provide an averaged value of the single-molecule conductance. Transient curves that are either noisy or that showed smooth exponential decay because of the absence of molecular bridge formation were rejected when building the histograms using an automatic selection procedure driven by a code written in LabVIEW. The histograms were compiled by applying the same automated selection criteria to each set of the recorded decay curves. The selection procedure allows current traces showing counts exceeding a dened threshold to be added to the conductance histogram. The percentage decay curves that showed clear molecular steps were typically 815% and were all selected for building the histograms. This selection process made peaks in the conductance

8 NATURE COMMUNICATIONS | 8:15195 | DOI: 10.1038/ncomms15195 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15195 ARTICLE

histograms more prominent above the tunnelling background and also allowed a quantitative measure of the yield of molecular junction formation.

Data availability. Full details on the synthesis and characterization of the foldamers and additional details regarding infrared ion dip spectroscopy, computer modelling and STM break junction are given in the Supplementary Information. Crystallographic data for trimer-NHBoc are deposited with the Cambridge Crystallographic Data Centre under reference number CCDC-1542042.

References

1. Gellman, S. H. Foldamers: a manifesto. Acc. Chem. Res. 31, 173180 (1998).2. Hill, D. J., Mio, M. J., Prince, R. B., Hughes, T. S. & Moore, J. S. A eld guide to foldamers. Chem. Rev. 101, 38934012 (2001).

3. (eds Hecht, S. & Huc, I. Foldamers: Structure, Properties, and Applications (Wiley, 2007).

4. Guichard, G. & Huc, I. Synthetic foldamers. Chem. Commun. 47, 59335941 (2011).

5. Zhang, D. W., Zhao, X., Hou, J. L. & Li, Z. T. Aromatic amide foldamers: structures, properties, and functions. Chem. Rev. 112, 52715316 (2012).

6. Das, A. & Ghosh, S. Supramolecular assemblies by charge-transfer interactions between donor and acceptor chromophores. Angew. Chem. Int. Ed. 53, 20382054 (2014).

7. Hartley, C. S. Folding of ortho-phenylenes. Acc. Chem. Res. 49, 646654 (2016).8. Yu, Z. & Hecht, S. Remote control over folding by light. Chem. Commun. 52, 66396653 (2016).

9. Porter, E. A., Wang, X., Lee, H.-S., Weisblum, B. & Gellman, S. H. Antibiotics: non-haemolytic [beta]-amino-acid oligomers. Nature 404, 565565 (2000).

10. Horne, W. S. et al. Structural and biological mimicry of protein surface recognition by a/b-peptide foldamers. Proc. Natl Acad. Sci. USA 106, 1475114756 (2009).

11. Lee, E. F. et al. High-resolution structural characterization of a helical a/b-peptide foldamer bound to the anti-apoptotic protein Bcl-xL. Angew.

Chem. Int. Ed. 48, 43184322 (2009).12. Shirude, P. S. et al. Macrocyclic and helical oligoamides as a new class of G-quadruplex ligands. J. Am. Chem. Soc. 129, 1189011891 (2007).

13. De Poli, M. et al. Conformational photoswitching of a synthetic peptide foldamer bound within a phospholipid bilayer. Science 352, 575580 (2016).

14. Brown, R. A., Diemer, V., Webb, S. J. & Clayden, J. End-to-end conformational communication through a synthetic purinergic receptor by ligand-induced helicity switching. Nat. Chem. 5, 853860 (2013).

15. Laursen, J. S., Harris, P., Fristrup, P. & Olsen, C. A. Triangular prism-shaped beta-peptoid helices as unique biomimetic scaffolds. Nat. Commun. 6, 7013 (2015).

16. Juwarker, H., Suk J.-M. & Jeong, K.-S. Foldamers with helical cavities for binding complementary guests. Chem. Soc. Rev. 38, 33163325 (2009).

17. Tanatani, A., Hughes, T. S. & Moore, J. S. Foldamers as dynamic receptors: probing the mechanism of molecular association between helical oligomers and rodlike ligands. Angew. Chem. Int. Ed. 41, 325328 (2002).

18. Gan, Q. et al. Helix-rod host-guest complexes with shuttling rates much faster than disassembly. Science 331, 11721175 (2011).

19. Berl, V., Huc, I., Khoury, R. G., Krische, M. J. & Lehn, J. M. Interconversion of single and double helices formed from synthetic molecular strands. Nature 407, 720723 (2000).

20. Kwon, S. et al. Magnetotactic molecular architectures from self-assembly of beta-peptide foldamers. Nat. Commun. 6, 8747 (2015).

21. Collie, G. W. et al. Shaping quaternary assemblies of water-soluble non-peptide helical foldamers by sequence manipulation. Nat. Chem. 7, 871878 (2015).

22. Marcos Ramos, A. et al. Supramolecular control over donor-acceptor photoinduced charge separation. J. Am. Chem. Soc. 126, 96309644 (2004).

23. Zeidan, T. A., Wang, Q., Fiebig, T. & Lewis, F. D. Molecular wire behavior in pi-stacked donor-bridge-acceptor tertiary arylureas. J. Am. Chem. Soc. 129, 98489849 (2007).

24. Wolffs, M. et al. Helical aromatic oligoamide foldamers as organizational scaffolds for photoinduced charge transfer. J. Am. Chem. Soc. 131, 48194829 (2009).

25. Li, X. et al. Photoinduced electron transfer and hole migration in nanosized helical aromatic oligoamide foldamers. J. Am. Chem. Soc. 138, 1356813578 (2016).

26. Nitzan, A. & Ratner, M. A. Electron transport in molecular wire junctions. Science 300, 13841389 (2003).

27. Tao, N. J. Electron transport in molecular junctions. Nat. Nanotechnol. 1, 173181 (2006).

28. Genereux, J. C. & Barton, J. K. Mechanisms for DNA charge transport. Chem. Rev. 110, 16421662 (2010).

29. Grozema, F. C. & Siebbeles, L. D. A. (eds) Charge and Exciton Transport through Molecular Wires (Wiley, 2011).

30. Aradhya, S. V. & Venkataraman, L. Single-molecule junctions beyond electronic transport. Nat. Nanotechnol. 8, 399410 (2013).

31. Guldi, D. M., Nishihara, H. & Venkataraman, L. Molecular wires. Chem. Soc. Rev. 44, 842844 (2015).

32. Bredas, J. L., Beljonne, D., Coropceanu, V. & Cornil, J. Charge-transfer and energy-transfer processes in pi-conjugated oligomers and polymers: a molecular picture. Chem. Rev. 104, 49715004 (2004).

33. Coropceanu, V. et al. Charge transport in organic semiconductors. Chem. Rev. 107, 926952 (2007).

34. Seferos, D. S., Blum, A. S., Kushmerick, J. G. & Bazan, G. C. Single-molecule charge-transport measurements that reveal technique-dependent perturbations.J. Am. Chem. Soc. 128, 1126011267 (2006).35. Molina-Ontoria, A. et al. [2,2]paracyclophane-based pi-conjugated molecular wires reveal molecular-junction behavior. J. Am. Chem. Soc. 133, 23702373 (2011).

36. Wielopolski, M. et al. Blending through-space and through-bond pi-pi-coupling in [2,2]-paracyclophane-oligophenylenevinylene molecular wires.J. Am. Chem. Soc. 135, 1037210381 (2013).37. Chen, L. et al. Multichannel conductance of folded single-molecule wires aided by through-space conjugation. Angew. Chem. Int. Ed. 54, 42314235 (2015).

38. Ullmann, K. et al. Single-molecule junctions with epitaxial graphene nanoelectrodes. Nano Lett. 15, 35123518 (2015).

39. Kang, Y. K., Rubtsov, Iovine, P. M., Chen, J. & Therien, M. J. Distance dependence of electron transfer in rigid, cofacially compressed, pi-stacked porphyrin-bridge-quinone systems. J. Am. Chem. Soc. 124, 82758279 (2002).

40. Vura-Weis, J. et al. Crossover from single-step tunneling to multistep hopping for molecular triplet energy transfer. Science 328, 15471550 (2010).

41. Batra, A. et al. Quantifying through-space charge transfer dynamics in pi-coupled molecular systems. Nat. Commun. 3, 1086 (2012).

42. Schneebeli, S. T. et al. Single-molecule conductance through multiple pi-pi-stacked benzene rings determined with direct electrode-to-benzene ring connections. J. Am. Chem. Soc. 133, 21362139 (2011).

43. Guo, X., Gorodetsky, A. A., Hone, J., Barton, J. K. & Nuckolls, C. Conductivity of a single DNA duplex bridging a carbon nanotube gap. Nat. Nanotechnol. 3, 163167 (2008).

44. Dulic, D. et al. Direct conductance measurements of short single DNA molecules in dry conditions. Nanotechnology 20, 115502 (2009).

45. van Zalinge, H. et al. Variable-temperature measurements of the single-molecule conductance of double-stranded DNA. Angew. Chem. Int. Ed. 45, 54995502 (2006).

46. Slinker, J. D., Muren, N. B., Renfrew, S. E. & Barton, J. K. DNA charge transport over 34 nm. Nat. Chem. 3, 228233 (2011).

47. Livshits, G. I. et al. Long-range charge transport in single G-quadruplex DNA molecules. Nat. Nanotechnol. 9, 10401046 (2014).

48. Artes, J. M., Li, Y., Qi, J., Anantram, M. P. & Hihath, J. Conformational gating of DNA conductance. Nat. Commun. 6, 8870 (2015).

49. Di Ventra, M. & Taniguchi, M. Decoding DNARNA and peptides with quantum tunnelling. Nat. Nanotechnol. 11, 117126 (2016).

50. Xu, B. & Tao, N. J. Measurement of single-molecule resistance by repeated formation of molecular junctions. Science 301, 12211223 (2003).

51. Chen, F. & Tao, N. J. Electron transport in single molecules: from benzene to graphene. Acc. Chem. Res. 42, 429438 (2009).

52. Anthony, J. E. Functionalized acenes and heteroacenes for organic electronics. Chem. Rev. 106, 50285048 (2006).

53. Odriozola, I., Kyritsakas, N. & Lehn, J. M. Structural codons: linearity/helicity interconversion by pyridine/pyrimidine exchange in molecular strands. Chem. Commun. 6263 (2004).

54. Venkataraman, L., Klare, J. E., Nuckolls, C., Hybertsen, M. S. & Steigerwald, M.L. Dependence of single-molecule junction conductance on molecular conformation. Nature 442, 904907 (2006).55. Abraham, M. H., Andonian-Haftvan, J., Whiting, G. S., Leo, A. & Taft, R. S. Hydrogen bonding. Part 34. The factors that inuence the solubility of gases and vapours in water at 298 K, and a new method for its determination.J. Chem. Soc., Perkin Trans. 2, 17771791 (1994).56. Sugino, M. et al. Elucidation of anthracene arrangement for excimer emission at ambient conditions. Cryst. Growth Des. 13, 49864992 (2013).

57. Cocinero, E. J., arabal, P., Vaden, T. D., Simons, J. P. & Davis, B. G. Sensing the anomeric effect in a solvent-free environment. Nature 469, 7679 (2011).

58. Len, I., Milln, J., Cocinero, E. J., Lesarri, A. & Fernndez, J. A. Shaping micelles: the interplay between hydrogen bonds and dispersive interactions. Angew. Chem. Int. Ed. 52, 77727775 (2013).

59. Hines, T. et al. Transition from tunneling to hopping in single molecular junctions by measuring length and temperature dependence. J. Am. Chem. Soc. 132, 1165811664 (2010).

60. Frisenda, R. & van der Zant, H. S. J. Transition from strong to weak electronic coupling in a single-molecule junction. Phys. Rev. Lett. 117, 126804 (2016).

NATURE COMMUNICATIONS | 8:15195 | DOI: 10.1038/ncomms15195 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15195

61. Komoto, Y., Fujii, S., Nishino, T. & Kiguchi, M. High electronic couplings of single mesitylene molecular junctions. Beilstein J. Nanotechnol. 6, 24312437 (2015).

62. Matsuhita, R., Horikawa, M., Naitoh, Y., Nakamura, H. & Kiguchi, M. Conductance and SERS measurement of benzenedithiol molecules bridging between Au electrodes. J. Phys. Chem. C 117, 17911795 (2013).

63. Rijs, A. & Oomens, J. in Gas-Phase IR Spectroscopy and Structure of Biological Molecules (eds Rijs, A. & Oomens, J.) (Springer, 2015).

64. MOPAC2016, version: 16.035L (Stewart Computational Chemistry, Colorado Springs, CO, USA, 2016) Available at http://OpenMOPAC.net

Web End =http://OpenMOPAC.net .

65. OBoyle, N. M. et al. Open Babel: an open chemical toolbox. J. Cheminformatics 3, 114 (2011).

66. Strutyski, K., Melle-Franco, M. & Gomes, JANF New parameterization scheme of DFT-D for graphitic materials. J. Phys. Chem. A. 117, 28442853 (2013).67. Neese, F. The ORCA program system. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2, 7378 (2012).

68. Frisch, M. J. et al. Gaussian 09 (Gaussian, Inc., Wallingford, Connecticut, USA, 2009).

69. Marco, A. B. et al. Bis(triisopropylsilylethynyl)-substituted pyrene-fused tetraazaheptacene: synthesis and properties. Phys. Chem. Chem. Phys. 18, 1161611619 (2016).

Acknowledgements

We thank Dr J.I. Miranda (University of the Basque Country) for assistance on NMR characterization. We are grateful to the Basque Science Foundation for Science (Ikerbasque), POLYMAT, the University of the Basque Country (SGIker and UFI11/23), the Deutsche Forschungsgemeinschaft (AU 373/3-1 and MA 5215/4-1), Gobiernode Espana (Ministerio de Economa y Competitividad CTQ2016-77970-R, CTQ2015-71936-REDT, CTQ2015-71406-ERC, CTQ2015-64579-C3-3P, CTQ-2014-54464-R and CTQ-2015-68148), Gobierno Vasco (BERC programme, Consolidated Groups IT520-10 and PC2015-1-01(06-37)), Diputacin Foral de Guipzcoa (OF215/2016(ES)), CICECOAveiro Institute of Materials, POCI-01-0145-FEDER-007679 (FCT ref. UID/CTM/50011/ 2013), ON2 (NORTE-07-0162-FEDER-000086) and the European Union (ERA Chemistry, Marie Curie Career Integration Grant No. 618247, and FET-Open (2D-INK) Grant No. 664878).

Author contributions

M.C. carried out the synthesis and characterization of the foldamers under the supervision of A.M.-A.; M.C. and M.P.R. carried out the STM break-junction studies under the supervision of I.D.-P.; I.U., J.A.F. and E.J.C. carried out the ion dip infrared spectroscopy measurements; M.M.-F. carried out the calculations to establish the conformation, the electronic structure and the infrared frequencies of the whole foldamer series; A.M.-A. conceived and supervised the study and wrote the paper with the input from all authors.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Carini, M. et al. High conductance values in p-folded molecular junctions. Nat. Commun. 8, 15195 doi: 10.1038/ncomms15195 (2017).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2017

10 NATURE COMMUNICATIONS | 8:15195 | DOI: 10.1038/ncomms15195 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications