Single-molecule detection of dihydroazulene photo-thermal reaction using break junction technique

Abstract

Charge transport by tunnelling is one of the most ubiquitous elementary processes in nature. Small structural changes in a molecular junction can lead to significant difference in the single-molecule electronic properties, offering a tremendous opportunity to examine a reaction on the single-molecule scale by monitoring the conductance changes. Here, we explore the potential of the single-molecule break junction technique in the detection of photo-thermal reaction processes of a photochromic dihydroazulene/vinylheptafulvene system. Statistical analysis of the break junction experiments provides a quantitative approach for probing the reaction kinetics and reversibility, including the occurrence of isomerization during the reaction. The product ratios observed when switching the system in the junction does not follow those observed in solution studies (both experiment and theory), suggesting that the junction environment was perturbing the process significantly. This study opens the possibility of using nano-structured environments like molecular junctions to tailor product ratios in chemical reactions.

ARTICLE

Received 19 Sep 2016 | Accepted 30 Mar 2017 | Published 22 May 2017

Cancan Huang1,2, Martyn Jevric3, Anders Borges3, Stine T. Olsen3, Joseph M. Hamill2, Jue-Ting Zheng1, Yang Yang1, Alexander Rudnev2, Masoud Baghernejad2, Peter Broekmann2, Anne Ugleholdt Petersen3, Thomas Wandlowski2,z, Kurt V. Mikkelsen3, Gemma C. Solomon3, Mogens Brndsted Nielsen3& Wenjing Hong1,2

Charge transport by tunnelling is one of the most ubiquitous elementary processes in nature. Small structural changes in a molecular junction can lead to signicant difference in the single-molecule electronic properties, offering a tremendous opportunity to examine a reaction on the single-molecule scale by monitoring the conductance changes. Here, we explore the potential of the single-molecule break junction technique in the detection of photo-thermal reaction processes of a photochromic dihydroazulene/vinylheptafulvene system. Statistical analysis of the break junction experiments provides a quantitative approach for probing the reaction kinetics and reversibility, including the occurrence of isomerization during the reaction. The product ratios observed when switching the system in the junction does not follow those observed in solution studies (both experiment and theory), suggesting that the junction environment was perturbing the process signicantly. This study opens the possibility of using nano-structured environments like molecular junctions to tailor product ratios in chemical reactions.

1 State Key Laboratory of Physical Chemistry of Solid States, College of Chemistry and Chemical Engineering, Collaborative Innovation Center of Chemistry for Energy Materials, Xiamen University, 361005 Xiamen, China. 2 Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland. 3 Department of Chemistry and Nano-Science Center, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen , Denmark. Correspondence and requests for materials should be addressed to W.H. (email: mailto:[email protected]

Web End [email protected] ) or to M.B.N. (email: mailto:[email protected]

Web End [email protected] ) or to G.C.S. (email: mailto:[email protected]

Web End [email protected] ) or to K.V.M. (email: mailto:[email protected]

Web End [email protected] ).

zDeceased.

NATURE COMMUNICATIONS | 8:15436 | DOI: 10.1038/ncomms15436 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

DOI: 10.1038/ncomms15436 OPEN

Single-molecule detection of dihydroazulene photo-thermal reaction using break junction technique

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15436

The switching of photochromic molecules is accompanied by changes in electronic, structural and/or chemical properties, making photoswitches versatile building blocks

for potential applications in materials science, electronics and biotechnology13. Progress in this area requires detailed understanding of these molecules and various analytical techniques have been used to analyse the photoreaction process, such as ultravioletvisible (UVVis) spectroscopy, uorescence spectroscopy4,5 and Raman spectroscopy6,7. As a potential analytical tool, the single-molecule break junction technique has been employed to study the conductancestructure correlation of a large range of molecular components in the past years814. More importantly, the break-junction techniques could also provide a means to statistically quantify the photo-reaction kinetics of a single-molecule device connected between two electrodes, which may offer some new understanding beyond measuring an ensemble of molecules as in other spectroscopies.

We became interested in exploring the dihydroazulene (dha)/vinylheptafulvene (vhf) system15,16 as previous studies have shown light-controlled conductance switching of dha derivatives in a silver nanogap fabricated using sublimable molecules and also bias-induced switching for a junction operating in the Coulomb blockade regime17,18. Fabrication of junctions by sublimation puts some constraints on the molecules that can be studied in regard to molecular weight and stability. The break-junction technique employs a solution of the molecules and it has allowed us to achieve charge transport properties of larger dha molecules by incorporating anchoring groups ( SAc) at each end of the molecule to control the

anchoring conguration.

We employ the mechanically controllable break junction (MCBJ) technique to detect the photo-thermal reaction processes of a tailor-made photochromic dha/vhf system. Three distinguishable conductance states are observed experimentally,

while reversible and distinct changing of the single molecular conductance is observed between two of the states. Beneting from the well-distinguished conductance states, the reaction kinetics and reversibility could be studied via statistical analysis of the conductancedistance traces, and it is found that the reaction in the junction does not follow those observed in solution, which agrees well with the energy calculation. Density functional transport calculations reproduce the relative conductance of these molecular states and identify a changing destructive interference effect as being responsible for the large span in conductance.

ResultsSingle-molecule conductance measurement. We measured the single-molecule conductance of the target molecules using the MCBJ technique under ambient conditions19,20. The MCBJ measurements proceeded as follows: the contacted point of the two electrodes was repeatedly broken and reformed in the solution containing 0.2 mM target molecule under the control of piezo-stack movement. During the breaking process, the anchoring groups (thioacetates) can contact the surface of gold electrodes via AuS bond and bridge the molecule between the separated electrodes to form the Au-molecule-Au junction. Under a xed bias (100 mV), the current through the junction was recorded and is used for the further analysis. Figure 1b shows the typical conductance traces plotted on a logarithmic scale. The plateaus with a conductance value below G0 (conductance quantum) are assigned to molecular junctions. Hundreds to thousands of conductance traces are used to construct the one-dimensional (1D) histograms and two-dimensional (2D) histograms.

To probe the photo-conversion process of dha-6, we introduced the in situ UV irradiation system to the MCBJ setup. In the initial state (dha-6 solution without UV irradiation), all the

a

In situ UV irradiation

c

0

30

15

dha-6

vhf

2

0

log(G/G 0)log(G/G 0)

4

1

6

8

0 1 2 3

z (nm)

2

d

0

b

2

1

5

0

5

dha-6 junction vhf junction

2

0

0

2

1

2

4

2

log(G/G 0)

6

4

4

log(G/G 0)

6

6

8

8

0

8 8 6 4 2 0

log(G/G0)

Figure 1 | Single-molecule conductance measurements. (a) Schematic of the mechanically controllable break junction (MCBJ) measurements during the photo-conversion of dihydroazulene under the in situ ultraviolet (UV) irradiation. The solution contains dihydroazulene (dha-6) and vinylheptafulvene (vhf). (b) Typical individual conductancedistance traces recorded in break junction measurements; blue for the dha-6 junction and red for the vhf junction. The applied bias potential is 100 mV. (c) 2D conductancedisplacement histogram constructed from 1,000 conductancedistance traces, the conductance region 10 3.0B10 4.5 G0 is dened as region 1, and the conductance region 10 4.5B10 6.0 G0 is dened as region 2. Linear scaling is applied for the colour bar. (d) 2D covariance histogram constructed from 1,000 conductancedistance traces.

2 NATURE COMMUNICATIONS | 8:15436 | DOI: 10.1038/ncomms15436 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15436 ARTICLE

a

10

M b

1.0

0

heating dha-6

0.8

10

Counts

M

P

0.6

0

0.4

M

10

0.2

0

50

100

150

0

10

20

30

Time (min)

0 0 1 2 3

z (nm)

z (nm)

d

c

4,000

20

0

Counts

Counts

2,000

20

0 0 2 0 2 0 2

0 2

0 2

z (nm)

0

8

6

4

2

0

log(G/G0)

Figure 2 | Photo-thermal reaction kinetics. (a) The typical relative displacement distributions of MCBJ conductance measurements with the in situ UV irradiation (30 min). The conductance regions of the top, middle and below section are 10 4.5B10 6 G0 (vhf), 10 3.0B10 4.5 G0 (dha-6) and whole region of 10 3.0B10 6 G0, respectively. The Gaussian tting is used to determine the area ratio of M and T, in which the peak marked as T is ascribed to the tunneling traces and M is for the molecular junction traces. (b) The dha-6 percent versus the reaction time. The upright triangle presents the data points determined from the in situ UV irradiation process, and tted with the exponential function, P P0 A exp ( R0t) (black solid curve), to calculate

the time decay constant. The inverted triangle shows the data points from the ex situ heating process, and is tted with P 1 (1 P0) exp( R0t) (black

dashed curve). (c) One-dimensional (1D) conductance histogram from different states. Blue is for the initial dha-6 state marked with blue solid box in b; red is for the state after in situ UV irradiation for 120 min, marked as red solid box in b; black is from the state after the heating for 30 min connected with the state marked as black solid box in b. The overlapping of blue and black curves suggests a highly reversible dha-62vhf conversion. Due to the limited time for recording during the in situ UV-MCBJ measurements, the numbers of conductance curves for the statistical analysis are limited to around 500.(d) Plateau distribution of the break junctions from different reaction times connected to b. The conductance region of vhf is 10 4.5B10 6 G0, and that of dha-6 is 10 3.0B10 4.5G0.

molecular junctions are formed by trapping the dha-6 into the nanogaps (dha-6 junction). The conductance traces show the plateau in the range 10 3.0B10 4.5 G0, and the maximum conductance peak around 10 3.70.1 G0 is obtained by a

Gaussian tting of the constructed 1D histogram. With UV irradiation (365 nm), the dha-6 undergoes a ring-opening process of the ve-membered ring to form vhf. Thus, both dha-6 and vhf exist in the solution simultaneously. During the conductance measurements, two kinds of conductance traces are observed, one with a plateau around 10 3.0B10 4.5 G0 (marked as region 1)

assigned to dha-6 junction (blue traces) and the other with a plateau around 10 4.5B10 6 G0 (marked as region 2), which is ascribed to vhf junctions. The 1D conductance histogram is plotted in Fig. 1b (right part), showing two distinct conductance peaks. The constructed 2D conductancedisplacement histogram (Fig. 1c) also shows two well-distinguished conductance features with more than one order of magnitude (1.4) conductance difference between them and almost the same plateau length. Moreover, we analysed the correlation between these two kinds of conductance traces by compiling the 2D covariance histogram in Fig. 1d (ref. 21). The signicant anti-correlation in the intersection of region 1 and region 2 indicates that the molecular junction is constructed either by dha-6 or vhf for each stretching cycle. Therefore, we can extract each type of conductance trace to quantitatively determine the percentages of dha-6 and vhf in solution.

Photo-thermal reaction kinetics. To further investigate the photoreaction process of the dha-6, we studied the reaction kinetics with in situ UV irradiation by statistical analysis of the percentage of dha-6 as a function of time. Based on the anti-correlation between dha-6 junction and vhf junction mentioned above, the ratio of dha-6/vhf or the content of each component (dha-6 or vhf) could be determined from counting the traces presenting each conductance character, although there are around 510% traces exhibiting direct tunnelling without forming a molecular junction.

Figure 2a shows the typical relative displacement distributions of conductance from MCBJ measurements with in situ UV irradiation (30 min). To determine the time-dependent junction evolution, 500 conductance curves are used for the statistical analysis for each time period. Using the conductance range of 10 3.0B10 4.5 G0 (region 1) and 10 4.5B10 6 G0 (region 2), we are able to construct the relative displacement distribution for dha and vhf, respectively. In the relative displacement distributions20 shown in Fig. 2a, the Gaussian tting of the two peak distributions suggests that the traces with molecular junction (marked with M) are typically longer than 1 nm, while the tunnelling traces without molecular plateau are typically at around 0.5 nm (marked with T). In this way, we could determine the peak area ratio of dha-6 or vhf junctions from the Gaussian tting of displacement distribution of all conductance traces including the dha-6 (or vhf) junction traces and tunnelling traces

NATURE COMMUNICATIONS | 8:15436 | DOI: 10.1038/ncomms15436 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15436

a b

dha-7 dha-6

3

0

30

0.8

Percent of dha-6

UV heating

4

log(G/G 0)

2

0

5

0.4

log(G/G 0)

log(G/G 0)

4

6

6

0.0

0

8 0 1 2 3 4 5 6

c

d

3

0.8

Percent of dha-6

4

Counts

4,000

5

2,000

0.4

6

0

0.0

8

6

4

2

0

7 0 1 2 3 4 5 6

log(G/G0)

Cycles

Figure 3 | Photo/thermal conversion of the isomer. (a) Plot of reversible conversion of dha-62vhf in six cycles. The solid circle means the conductance, related to left scale, while the hollow circle is the percentage of dha-6, related to right scale. (b) Schematic and 2D conductancedisplacement histograms for the conversion of dha-7 with the treatments by UV irradiation and heating. (c) One-dimensional (1D) conductance histogram from different states. The black is from the initial dha-7 state, the red is for the state after UV irradiation to vhf and the blue is from the conductance measurements for solution after heating process. (d) Plot of conductance switch cycles starting from dha-7. The error bars of the conductance are determined from the conductance analysis based on Gaussian Function (s.d.). The solid circle means the conductance, related to left scale, while the hollow circle is the percentage of dha-6, related to right scale.

(see Supplementary Fig. 22 for more details). Therefore, the percentage of dha-6 junctions in all molecular junctions could be determined for the further kinetics analysis.

We use the combination of the in situ UV irradiation and conductance measurements to investigate the photo-conversion kinetics. The typical stretching distance distributions are presented in Fig. 2d, which are constructed from different states during the conversion process. The percentage of dha-6, along with the irradiation time, is plotted in Fig. 2b, showing that the vhf content reached saturation after 120 min of UV irradiation. The data (upright triangles) were tted with the exponential function, P P0 A exp( R0t), where P0 is the nal value of

the percentage, R0 is the decay constant of the best-t curve with units of time 1 used for the calculation of time constant by taking the inverse. Therefore, the time constant for dha-6-vhf photoreaction is calculated as 455 min, which is higher than the previous results in the solution (less than 5 min)22,23. Two factors could be responsible for the longer time to equilibrate. Firstly, during the UV irradiation, heating caused by the light is inevitable, which leads to the back reaction (vhf-dha-6). More importantly, in the molecular junction, a localized surface plasmon resonance will form in nano-gaps (1.53 nm) between two Au electrodes under bias potential, leading to a large buildup of photon intensity and high concentration of energetic electrons11,24. During the in situ UV-MCBJ measurements, the localized surface plasmon resonance excitation effect could compete with the absorption of the photon from the irradiation to accelerate the back reaction. Figure 2b also illustrates that there is no complete conversion from dha-6 to vhf, that is, we see the decay curve plateau at just under 40% dha-6. We therefore speculate that the heating effect from irradiation and plasmon non-radiative decay must increase the back reaction rate under the in situ UV-MCBJ measurements. A consequence of this is

that longer irradiation cannot further decrease the percentage of dha-6. It also illustrates that we could not obtain the complete conversion from dha-6 to vhf. Therefore, the heating effect from irradiation and plasmon non-radiative decay could increase the back reaction rate under the in situ UV-MCBJ measurements. It is also conrmed by the experiments with longer time of UV irradiation (Supplementary Fig. 23).

To understand the reversible ring-closing process, the ex situ heating of B60 C was applied. As plotted in Fig. 2b (inverted triangle), an increasing percentage of dha-6 with heating time was observed and the vhf fully switched back to dha-6 (noted as dha-60 below for clarity) in 30 min. Furthermore, the function of

P 1 (1 P0) exp( R0t) was applied for the tting of the

thermal process, and the disagreement between the tting and the experiments suggests that the constraints of the junction may alter the kinetics of the back reaction. The 1D histograms from three states, the initial dha-6, vhf after in situ UV irradiation for 120 min, and the dha-60 after heating 30 min from vhf state, are constructed and displayed in Fig. 2c. It is shown that the conductance peak of dha-60 overlaps perfectly with dha-6, suggesting a highly reversible dha-62vhf conversion.

Photo/thermal conversion of the isomer. For the conversion of dha-62vhf, the reversible switching cycles are performed sequentially. We demonstrate that there is no signicant conversion attenuation in conductance ratio; thus, Gdha-6/Gvhf

remained approximately 1.4, as shown in Fig. 3a. Therefore, dha-62vhf can be considered as a highly promising single-molecule switch. However, the conversion ratio of dha-6 shows slight attenuation from the analysis of plateau distribution,i.e., we see slightly less dha-6 after every switching cycle. This may be due to the heat accumulation in the nano-gap regime, which

4 NATURE COMMUNICATIONS | 8:15436 | DOI: 10.1038/ncomms15436 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15436 ARTICLE

a

b

NC

CN

a c

dha-6

vhf

dha-7

AcS

SAc

TS7 (28.88)

S

TS6 (28.58)

No interference

Shifted interference

N

N

heating

UV

UV

dha-6

dha-7

E-vhf

Relative gibbs free

energy (kcal mol1 )

10

R

R

S

20

R

CN

CN

S

NC CN

SAc

AcS

AcS

CN SAc

NC

CN

CN

Interference

1

8

7

R

R

N

N

CN

CN

CN

R

R R

Rot.

CN

2

S

(5.15)

(5.19) (5.18)

3 4 5

6

R

R

0

(0.39) (3.12) (3.26) (0.0)

CN CN

NC

S

N

Azulene numbering

NC

R R

N

(s-cis)

S

AcS SAc

b

100

101

102

NC

CN

Figure 4 | The structure of the molecules studied in this work.(a) The numbering used around azulene. (b) The molecular structures of the systems studied and the conversions observed in junction. The system can switch among dha-6, vhf and dha-7, where the numbers 6 and 7 refer to the position of anchoring group substitution in the seven-membered ring of dha.

30

Transmission

103

104

1053 2 1 0

E EF (eV)

1 2 3

signicantly shifts the reaction equilibrium of the dha-62vhf conversion. Together, these two results suggest that while two clear, reversible conductance states are observable at the single-molecule level, incomplete conversion might complicate the operation of a device utilizing an ensemble of molecules, and further work is required to optimize the percentage conversion in a junction.

As introduced in Fig. 4, switching between vhf and dha-6 is not the only possible reaction for this system as the reverse reaction can, in principle, lead to dha-7. As we did not see evidence of this conversion when starting from dha-6, dha-7 was prepared by chemical synthesis25 and taken as the precursor. During the conductance measurements, the operations of UV irradiation and thermal heating were also applied to dha-7. The 2D conductance displacement and 1D conductance histograms from one photo-thermal cycle are shown in Figs 3b and 3c. The conductance measurements demonstrated that the conductance of dha-7 (10 5.80.1 G0) is lower than that of dha-6, and it is also even lower than that of the ring-opened vhf state. More interestingly, in the case of the dha-7, after one photo/heating cycle, we only obtained one conductance peak, at G 10 3.70.1 G0, implying

an almost complete conversion back to dha-6. Thus we conclude that the isomerization dha-7- vhf-dha-6 occurred during the UV irradiation and heating reaction cycle. As shown in Fig. 3b, the UV irradiation triggers the ring-opening reaction and vhf is formed, followed by equilibration between vhf isomers. During the ring-closing reaction under the thermal heating, however, only dha-6 is formed. This can be seen in the 2D histogram, which exhibits an intensive conductance cloud around G 10 3.70.1 G0, ascribed to dha-6. After the rst complete

switching cycle, almost no conductance features of dha-7 are observed, shown in Fig. 3b (right). Afterwards, six photo-thermal cycles were performed. It was clear that the conductive properties are quite similar with dha-6, including the reversibility and switching, as depicted in Fig. 3d. The conductance switched between high conductance (dha-6, 10 3.70.1G0) and low conductance (vhf, 10 5.10.1G0) with ratio of around 25.

While complete conversion of dha-7 to dha-6 was observed in the junction, which is probably promoted by a more favorable anchoring conguration of dha-6 or by a difference in the vhf to dha transition state dipole moments (vide infra), no such preference for dha-6 has been observed in solution. Irradiating

a solution of dha-7 in CD3CN at 365 nm for ca. 2 h gave a mixture of E/Z isomeric vhfs according to 1H-NMR spectroscopy (see Supplementary Information). After 2 days in the dark, the vhfs underwent conversion to a mixture of dha-7 and dha-6 in a ratio of 4:3. Isomerization was also conrmed by UVVis absorption spectra in solution. Thus, subjecting a pristine sample of either dha-6 (lmax 393 nm) or dha-7 (lmax 372 nm) to one

opening-closure cycle in CH3CN gave a mixture of isomers indicated by an intermediate absorption maximum at ca. 380 nm.

Density functional theory calculations. In an effort to gain insight into the reaction kinetics of the switch and the discrepancy between the behaviours in solution and junction in regard to isomerization, we investigated the potential energy surface connecting the vhf and dha isomers using density functional theory calculations. Information on the details is given in the Supplementary Information. The energy prole of the system (employing acetonitrile as solvent) is shown in Fig. 5a. We nd that the free energies of dha-6 and dha-7 are very similar which is in good agreement with the outcome from experimental solution studies. The vhf isomers are connected to the dha species via transition states with similar energy barriers exceeding 20 kcal mol 1. The barriers for conversion between vhf isomers are much lower and we can therefore assume chemical equilibrium between the vhf isomers at all times in solution. Since the conversion from vhf to either dha is thermally activated, it is highly sensitive to the barrier heights. Any asymmetry in the barrier height would therefore lead to the promotion of one dha isomer over the other before chemical equilibrium is reached. The agreement between the experimental and theoretical solution

Figure 5 | Calculated energy and transport properties. (a) Potentialfree energy prole of the dha and vhf species in acetonitrile solution.(b) Calculated transmissions for structures are shown on the right. Zero-bias conductance indicates G(dha-6)4G(vhf)4G(dha-7), in agreement with experiment. (c) Example diagrams. No interference is expected at the

Fermi energy for dha-6 because all diagrams without onsite loops have the same sign. For dha-7, all diagrams contain the same onsite loop so destructive interference is expected at the Fermi energy. For vhf, diagrams with onsite loops on N shift the interference away from the Fermi energy.

NATURE COMMUNICATIONS | 8:15436 | DOI: 10.1038/ncomms15436 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15436

Figure 6 | X-Ray crystal structure. Molecular structure of dha-7 with displacement ellipsoids of 50% for non H-atoms (CCDC 1501743); red: oxygen; yellow: sulfur; blue: nitrogen. Crystals were grown from dichloromethane and heptane.

studies suggests that it is something peculiar to the junction environment that leads to the switching asymmetry. One notable difference between metal electrodes and solution is the ability of electrodes to polarize and screen charges, thereby stabilizing polar/charged systems. We nd that the dipole moment of the transition state of Z-vhf-dha-7 is slightly smaller than for E-vhf-dha-6. We therefore speculate that screening of the dipole moment from electrons in the leads or enhanced electric elds associated with the Au nano-gap might promote the local formation of dha-6 relative to dha-7. Since spontaneous ring-opening should occur at time-scales much longer than the experiment, this would explain why dha-6 is preferentially observed in the break-junction experiment.

For comparing the energies of the two molecules dha-6 and dha-7 when placed in a junction (ensuring interactions between molecule and metal leads), we performed combined QM/MM calculations26 using the program Dalton QCP27. The gold electrodes are built as two hemispheres from a fcc unit cell and described using molecular mechanics. Each gold cluster consists of 262 Au atoms, where each gold atom is assigned a polarizability of 31.04 a.u. The molecule faces the [111] surface with the axis through the two sulfur atoms oriented perpendicular to the gold surface. The density functional theory geometry-optimized molecules were inserted between the two gold clusters and then single point energy calculations were performed (CAM-B3LYP//cc-pVDZ). Interestingly, we nd that dha-6 is more stable than dha-7 by 3.04 kcal mol 1 in the junction, while the difference was only 0.39 kcal mol 1 for the molecules in solution (acetonitrile). A smaller difference between the reactive s-cis Z-vhf and s-cis E-vhf isomers was obtained, with the Z-vhf being 0.72 kcal mol 1 more stable than the E-vhf. Although the slightly more energetic E-vhf is the precursor for dha-6, it pays off overall from a thermodynamic point of view to form dha-6 in the junction. We were not able to perform reliable TS calculations in the junction in order to ascertain whether the kinetics of ring closure also agrees with experiment. Nevertheless, we have previously observed that the larger the dha vhf energy

difference, the faster is the vhf-to-dha conversion2729.

Despite the chemical similarities of dha-6, vhf and dha-7, their calculated transport properties are vastly different and followed the trend dha-64vhf4dha-7 in the experiments. A theoretical analysis using a diagrammatic approach rst used by Markussen et al.30,31 combined with density functional theory transport calculations shows that dha-6 has the highest conductance because it is almost planar and does not show destructive interference (Fig. 5b). Conversely, dha-7 and vhf are bent and show destructive interference. The difference between these two systems results from the different energies at which the

interference feature occurs. The diagrammatic description predicts no destructive interference for dha-6 because it is possible to pair up all the non-traversed p-orbitals along the path connecting the electrodes for these molecules. However, this is not possible for dha-7 and vhf as illustrated in Fig. 5c. In dha-7, the anchoring group is placed at C7, which causes signicant deviation of the p-system from planarity as also conrmed by

X-ray crystallographic analysis (Fig. 6). Moreover, the cyano groups are electronically decoupled from the rest of the p-system by saturated bonds and destructive interference therefore occurs in the middle of the band gap. This decoupling is not present in vhf and the anti-resonance is therefore shifted away from the middle of the band gap due to the electron-accepting character of CN. The calculated features of dha-6 and dha-7 are in line with their different optical properties determined experimentally (vide supra): dha-6 exhibits a red-shifted longest-wavelength absorption maximum relative to dha-7.

DiscussionTo conclude, we have developed an approach to detect the photo-thermal reaction of dha/vhf system using the single-molecule break junction technique. We have demonstrated that this technique is able to track the photo-thermal reaction, and to evaluate the reversibility and multiple switching possibilities during the reaction using the well-distinguished conductance states of single-molecule dha/vhf junctions. The combined transport calculations suggest that these distinct conductance states were ascribed to the introduction and shifting of interference features in the molecular transmission.

More importantly, beyond simply the application as a molecular switch, this study raises questions about tuning reaction kinetics, barriers and ultimately product ratios in the novel environment of molecular junctions. While we saw clear and reversible switching, the product ratio obtained does not correspond with that observed in solution studies. It remains a question for future work whether the product ratios can be further tuned by junction engineering or if this principle can be employed to tune product ratios for other reactions.

Methods

Single-molecule conductance measurement. Conductance measurementswere performed using the MCBJ technique with a home-built MCBJ setup19,20. To perform the conductance measurement, the solution contains 0.2 mM target molecules in mixture solvent of THF:Mesitylene (TMB) 1:4. A blank experiment

of solvent without target molecule is presented in Supplementary Fig. 21.

UV irradiation and heating experiment. To have the in situ UV irradiation experiment, one UV LED at 365 nm is built with the power of 300 mW cm 2.

The home-built UV LED was put on the top of the solution with a distance of B1 cm. As for the photo-conversion kinetic studies, the UV LED were turned on during the conductance measurements. In the ex situ heating process, the solution was put in a water bath with constant temperature (65 C). After setting waiting time (0, 5, 15, 25 and 30 min), we took 200 ml solution into the liquid cell of MCBJ for conductance measurements.

Synthesis and characterization. Synthesis and NMR characterization are provided in Supplementary Figs 15 and Supplementary Note 1. UVVis and NMR spectroscopic studies in solution are provided in Supplementary Figs 610 and Supplementary Note 2. Details on the theoretical work are provided in Supplementary Figs 1120, Supplementary Table 1 and Supplementary Notes 36.

Data availability. The X-ray crystallographic coordinates for the structure reported in this article have been deposited at the Cambridge Crystallographic Data Centre (CCDC), under deposition number CCDC 1501743. These data can be obtained free of charge from The Cambridge Crystallographic Data Centre via http://www.ccdc.cam.ac.uk/data_request/cif

Web End =www.ccdc.cam.ac.uk/data_request/cif . The data that support the ndings of this study are available from the corresponding author upon reasonable request.

6 NATURE COMMUNICATIONS | 8:15436 | DOI: 10.1038/ncomms15436 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15436 ARTICLE

References

1. Irie, M., Fukaminato, T., Sasaki, T., Tamai, N. & Kawai, T. Organic chemistry: a digital uorescent molecular photoswitch. Nature 420, 759760 (2002).

2. Jia, C. C. et al. Covalently bonded single-molecule junctions with stable and reversible photoswitched conductivity. Science 352, 14431445 (2016).

3. Tian, H. & Yang, S. J. Recent progresses on diarylethene based photochromic switches. Chem. Soc. Rev. 33, 8597 (2004).

4. Vaya, I., Lhiaubet-Vallet, V., Jimenez, M. C. & Miranda, M. A. Photoactive assemblies of organic compounds and biomolecules: drug-protein supramolecular systems. Chem. Soc. Rev. 43, 41024122 (2014).

5. van Thor, J. J. Photoreactions and dynamics of the green uorescent protein. Chem. Soc. Rev. 38, 29352950 (2009).

6. Li, J. F. et al. Shell-isolated nanoparticle-enhanced Raman spectroscopy. Nature 464, 392395 (2010).

7. Nie, S. M. & Emery, S. R. Probing single molecules and single nanoparticles by surface-enhanced Raman scattering. Science 275, 11021106 (1997).

8. Aragones, A. C. et al. Electrostatic catalysis of a Diels-Alder reaction. Nature 531, 8891 (2016).

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

10. Jia, C. & Guo, X. Molecule-electrode interfaces in molecular electronic devices. Chem. Soc. Rev. 42, 56425660 (2013).

11. Cheng, Z. L. et al. In situ formation of highly conducting covalentAu-C contacts for single-molecule junctions. Nat. Nanotechnol. 6, 353357 (2011).

12. Sun, L. et al. Single-molecule electronics: from chemical design to functional devices. Chem. Soc. Rev. 43, 73787411 (2014).

13. Huang, C., Rudnev, A. V., Hong, W. & Wandlowski, T. Break junction under electrochemical gating: testbed for single-molecule electronics. Chem. Soc. Rev. 44, 889901 (2015).

14. Frisenda, R., Janssen, V. A. E. C., Grozema, F. C., van der Zant, H. S. J. & Renaud, N. Mechanically controlled quantum interference in individual p-stacked dimers. Nat. Chem 8, 10991104 (2016).

15. Boggio-Pasqua, M., Bearpark, M. J., Hunt, P. A. & Robb, M. A. Dihydroazulene/vinylheptafulvene photochromism: a model for one-way photochemistry via a conical intersection. J. Am. Chem. Soc. 124, 14561470 (2002).

16. De Waele, V., Schmidhammer, U., Mrozek, T., Daub, J. & Riedle, E. Ultrafast bidirectional dihydroazulene/vinylheptafulvene (DHA/VHF) molecular switches: photochemical ring closure of vinylheptafulvene proven by a two-pulse experiment. J. Am. Chem. Soc. 124, 24382439 (2002).

17. Lara-Avila, S. et al. Light-triggered conductance switching in single-molecule dihydroazulene/vinylheptafulvene junctions. J. Phys. Chem. C 115, 1837218377 (2011).

18. Broman, S. L. et al. Dihydroazulene photoswitch operating in sequential tunneling regime: synthesis and single-molecule junction studies. Adv. Funct. Mater. 22, 42494258 (2012).

19. Hong, W. et al. Single molecular conductance of tolanes: experimental and theoretical study on the junction evolution dependent on the anchoring group.J. Am. Chem. Soc. 134, 22922304.20. Hong, W. et al. An MCBJ case study: the inuence of pi-conjugation on the single-molecule conductance at a solid/liquid interface. Beil. J. Nanotechnol 2, 699713.

21. Makk, P. et al. Correlation analysis of atomic and single-molecule junction conductance. ACS Nano 6, 34113423 (2012).

22. Li, T. et al. Ultrathin reduced graphene oxide lms as transparent top-contacts for light switchable solid-state molecular junctions. Adv. Mater. 25, 41644170 (2013).

23. Pathem, B. K. et al. Photoreaction of matrix-isolated dihydroazulenefunctionalized molecules on Au{111}. Nano Lett. 13, 337343 (2013).

24. Chen, X.-J., Cabello, G., Wu, D.-Y. & Tian, Z.-Q. Surface-enhanced Raman spectroscopy toward application in plasmonic photocatalysis on metal nanostructures. J. Photochem. Photobiol. C 21, 5480 (2014).

25. Jevric, M., Broman, S. L. & Nielsen, M. B. Palladium-mediated strategies for functionalizing the dihydroazulene photoswitch: paving the way for its exploitation in molecular electronics. J. Org. Chem. 78, 43484356 (2013).

26. Olsen, S. T., Arcisauskaite, V., Hansen, T., Kongsted, J. & Mikkelsen, K. V. Computational assignment of redox states to Coulomb blockade diamonds. Phys. Chem. Chem. Phys. 16, 1747317478 (2014).

27. Aidas, K. et al. The DALTON quantum chemistry program system. WIREs Comput. Mol. Sci 4, 269284 (2014).

28. Cacciarini, M. et al. Towards solar energy storage in the photochromic dihydroazulene-vinylheptafulvene system. Chem. Eur. J. 21, 74547461 (2015).

29. Skov, A. B. et al. Towards storage of solar energy in photochromic molecules: benzannulation of the dihydroazulene/vinylheptafulvene couple. ChemPhotoChem. doi: 10.1002/cptc.201600046.

30. Markussen, T., Stadler, R. & Thygesen, K. S. The relation between structure and quantum interference in single molecule junctions. Nano Lett. 10, 42604265 (2010).

31. Markussen, T., Stadler, R. & Thygesen, K. S. Graphical prediction of quantum interference-induced transmission nodes in functionalized organic molecules. Phys. Chem. Chem. Phys. 13, 1431114317 (2011).

Acknowledgements

This work was generously supported by the University of Copenhagen, the Danish e-Infrastructure Cooperation, the European Union Seventh Framework Programme (FP7/2007-2013) under the ERC grant agreement no. 258806, the Danish Council for Independent ResearchNatural Sciences, the Carlsberg foundation, NSFC (21673195, 21503179), EC FP7 ITNs MOLESCO project numbers 606728, and the Young Thousand Talent Project of China.

Author contributions

W.H. and M.B.N. originally conceived the concept and designed the experiments. C.H., K.V.M., G.C.S., M.B.N., J.-T.Z., Y.Y. and W.H. prepared the manuscript using feedback from other authors. Break-junction measurements were carried out in the laboratory of T.W. and W.H. by C.H., J.H., A.R. and M.B.; synthetic work and spectroscopic studies in solution of the molecules were carried out in the laboratoryof M.B.N. by M.J. and A.U.P.; calculations were carried out in the group of K.V.M. and G.C.S. by A.B. and S.T.O. X-ray crystallographic analysis was performed by A.U.P. All authors have given approval to the nal version of the manuscript.

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: Huang, C. et al. Single-molecule detection of dihydroazulene photo-thermal reaction using break junction technique. Nat. Commun. 8, 15436doi: 10.1038/ncomms15436 (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

NATURE COMMUNICATIONS | 8:15436 | DOI: 10.1038/ncomms15436 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7