Reverse Bond-length Alternation in Cumulenes: Candidates for

Publication Date (Web): August 1, 2018 ... progress in cumulene synthesis, we discuss substituent strategies that may increase reverse bond-length alt...
0 downloads 0 Views 2MB Size
Subscriber access provided by University of Winnipeg Library

Feature Article

Reverse Bond-length Alternation in Cumulenes: Candidates for Increasing Electronic Transmission with Length Marc H. Garner, William Bro-Jørgensen, Pernille D. Pedersen, and Gemma C. Solomon J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b05661 • Publication Date (Web): 01 Aug 2018 Downloaded from http://pubs.acs.org on August 3, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Reverse Bond-length Alternation in Cumulenes: Candidates for Increasing Electronic Transmission with Length Marc H. Garner, William Bro-Jørgensen, Pernille D. Pedersen, Gemma C. Solomon* Department of Chemistry and Nano-Science Center, University of Copenhagen, Universitetsparken 5, DK-2100, Copenhagen Ø, Denmark.

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 47

Abstract. Single-molecule conductance generally decays exponentially with the length of the molecule when the transport mechanism is a coherent tunneling process. However, it was recently found that this length dependence can be reversed in linear conjugated molecules if the bond-length alternation is reversed. In this work we show that even-carbon cumulenes show this behavior as the bond-lengths are reversed for the dominant π system compared to the equivalent polyenes and polyynes. We explore the electronic origins of the reversed bond-length alternation in cumulenes and its relation to the length dependence of the electronic transmission. Through density functional theory and non-equilibrium Greens function calculations we predict that cumulenic wires have reverse decay of transmission with length, i.e., the decay constant β is found to be negative. As a direct consequence of the reversed bond-length alternation, the electronic transmission increases with length as the HOMO-LUMO gap rapidly narrows. Based on recent progress in cumulene synthesis, we discuss substituent strategies that may increase reverse bond-length alternation. Cumulenes stand out as promising candidates for a series of molecules that may show reverse decay of single-molecule conductance with increasing length.

2 ACS Paragon Plus Environment

Page 3 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Introduction The tunneling probability of a quantum mechanical particle decays exponentially with the length of the potential energy barrier, a trend that is seen in electron transfer and transport experiments for a range of molecular systems1-7. In single-molecule junction experiments the low-bias conductance between two metal electrodes is measured through a bridging molecule8-10, essentially constituting the tunneling barrier for the system. The attenuation of the single-molecule conductance, G, with length is generally described by  =  · exp (− · )

(1)

where g0 is a pre-factor that depends, in part, on the interface between the molecule and the electrodes, L is length in arbitrary length units or number of atoms, and β is the exponential attenuation factor11-16. When the tunneling transmission between two electrodes through a linear πconjugated molecule is modelled at the Hückel level of theory, β can be related to the bond-length alternation  = −2 ln   

(2)

where β is the exponential attenuation factor, and ts and td are the resonance integrals of the formal single and double bonds, respectively.17 However, Tsuji et al. noted this relationship has a peculiar implication: If the bond-length alternation of the molecule is reversed, β will be negative, and the tunneling transmission will increase with increasing molecular length.17 In linear conjugated molecules the ground-state geometry has formal single bonds that are longer than the double bonds, and consequently ts < td. In a geometry without bond-length alternation (ts = td), the molecule is metallic and there is no decay with length. Beyond the metallic limit, a regime of reverse tunneling decay may exist for molecules where ts > td, and thus the tunneling

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 47

transmission increases with the length of the molecule. Since the first suggestion of a reverse transmission decay system by Tada and Yoshizawa18, different reverse and decay-free systems have been suggested, including armchair graphene nanoribbons19, ethyndiyl-linked and fused oligomers of porphyrin20-24, and quinoidal molecules25-27. Still, reverse decay of single-molecule conductance with length has yet to be verified experimentally and the search for stable molecules with reverse bond-length alternation continues.28-29 Here we suggest cumulenes as a potential candidate for a series of molecules with reverse bondlength alternation in their stable ground-state geometries. While it has been duly noted in a number of studies that cumulenes show weak bond-length alternation30-41, the fact that the alternation is reversed and may correspond well to the situation described by Tsuji et al.17 has not yet been explored. In this work, we examine the possibility of reverse decay of tunneling transmission in evencarbon cumulenes due to reverse bond-length alternation. We study the molecules through density functional theory (DFT) and Hückel model calculations. All DFT calculations presented in the manuscript are done in the Atomic Simulation Environment42 (ASE) with the PBE exchangecorrelation functional43 and double-ζ plus polarization basis set for all atoms as implemented in GPAW.44-45 A benchmark study for this method is included in Supporting Information part A. The paper proceeds as follows. First, follows an introduction to the electronic structure of cumulenes. Second, we explore reverse bond-length alternation and conductance decay with simple Hückel models. Next, we present DFT calculations of the single-molecule transmission for alkene, cumulene, and alkyne series. This is followed by a brief exploration and discussion of chemical strategies for increasing the reverse bond-length alternation effect in cumulenes, and finally we summarize the findings.

4 ACS Paragon Plus Environment

Page 5 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Electronic structure of cumulenes Cumulenes are linearly π-conjugated molecules with all-double bond Kekulé structures. Odd n [n]cumulenes (n = 3, 5 …), where n is the number of cumulated double bonds so that the number of carbon atoms is n+1, have been synthesized and characterized up to [9]cumulene.46-52 We will denote these simply as odd [n]cumulenes. We only treat the odd [n]cumulenes here, i.e., the evencarbon series; the odd-carbon cumulenes will be treated in separate studies due to their helical orbital structure.53-55 Scheme 1: Two Orthogonal π systems of [3]cumulene

The basis of p orbitals in cumulenes has two important features. First, as shown in Scheme 1, the two π systems are orthogonal. Different from the even [n]cumulenes, for which these can mix and form helical π orbitals as we have described in detail elsewhere55, the π systems of odd [n]cumulenes remain separated by symmetry also when the molecule is substituted. The second important feature is that the two π systems are displaced by one carbon atom relative to each other.32 The long π-system (blue in Scheme 1) spans n+1 atoms and lies out of the hydrogen plane of the molecule; the short π-system (red in Scheme 1) spans n–1 atoms and lies in the plane of the

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 47

hydrogens. In this sense, the central double bond in [3]cumulene lies in the plane of the hydrogens, while the two double-bonds at the edges are equivalent to those we know from butadiene. In electron transport experiments the electrodes will typically primarily couple into one of the two π-systems as they are separated by symmetry. Though dependent on the specific molecular design and experimental method, it is most likely that the transport properties will be dominated by the out-of-plane π-system as this extends fully onto the terminal carbon atoms. As we shall see, it is in this π-system the bond-lengths are “reversed” in the sense that unlike the equivalent alkene and alkyne the shortest bonds are no longer at the edges of the molecule. Bond-length Alternation in Cumulenes. In Figure 1 we show the calculated structures for alltrans alkene, odd [n]cumulene, and alkyne molecules. Note that we use non-standard nomenclature for the alkenes and alkynes. For clarity, we use the same n-index for all three series, corresponding to the number of cumulated double bonds in the cumulene; the total number of carbon atoms is n+1 in all cases. In Figure 1b, the optimized structures and bond lengths are plotted for [3]alkene (butadiene), [3]cumulene (butatriene), and [3]alkyne (butadiyne). The difference between these three linear carbon molecules is just two hydrogen atoms between each. While the alkene and alkyne have the shortest bond at the edges as expected from their Kekulé structures, the cumulene has the longest bond at the edges and the shortest bond in the middle.

6 ACS Paragon Plus Environment

Page 7 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1. DFT-calculated bond lengths. (a) Chemical structures of [n]alkene, [n]cumulene, and [n]alkyne, where x is a positive integer and n = 2x+1. (b) Optimized structures and bond lengths in Ångström for [3]alkene, [3]cumulene, and [3]alkyne. (c) Bond lengths of [11]alkene, [11]cumulene, and [11]alkyne. In Figure 1c, we show the calculated bond lengths for n = 11 for the three series of molecules. Again, the first and last bonds are the shortest for [11]alkene and [11]alkyne in agreement with the Kekulé structures. For [11]cumulene the bond length alternation is very small in the center of the molecule, and the molecule approaches the true cumulenic structure with equal bond lengths. The weak alternation continues for the longer members of the cumulene series and is consistently reverse compared with alkenes and alkynes31, 33, 50. Worth noting, though all bonds in the Kekulé structures of cumulenes are denoted as double-bonds, the carbon atoms are sp-hybridized and consequently the bond lengths are very similar to those of the alkynes.

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 47

The calculated (reverse) bond-length alternation in the cumulene series is somewhat underestimated compared with calculations at higher levels of theory as presented in Supporting Information part A and crystallographic data.50 The results presented herein should therefore be considered a lower-boundary of the reverse conductance decay effect. The reverse bond-length alternation trend is consistent throughout the cumulene series studied here (up to n=19). Alkynes and cumulenes are both molecular forms of the “infinite” 1D material carbyne56-59. These have received much attention recently and their mechanical and electronic properties have been studied extensively60-65. In relation to carbyne, we note that reverse bond-length alternation is an edge effect; in a periodic structure there is either bond-length alternation or no bond-length alternation.66-68 In this sense, reverse bond-length alternation is purely a chemical concept, because chemists expect alternation between single and double bonds as drawn in Kekulé structures. Consequently, chemists expect linear conjugated molecules to start and end with double or triple bonds; if not, the molecule is a (di)radical. Frontier Molecular Orbitals. What is the origin of the reverse bond-length alternation in cumulenes? In Figure 2a we examine the frontier molecular orbitals of [5]cumulene, which consists of two orthogonal π-systems. Again, the long π-system is n+1 atoms long; the short π-system in n– 1 atoms long, now with some hyperconjugation onto the hydrogen atoms. These two π-systems are separated by symmetry and their nodal structure is displaced in space.55

8 ACS Paragon Plus Environment

Page 9 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2. Electronic structure calculated with DFT. (a) Molecular structure and frontier molecular orbitals of [5]cumulene animated with iso-value = 0.02. (b) HOMO-LUMO gap as function of length for [n]alkene, [n]cumulene, and [n]alkyne. In alkenes and alkynes, bond-length alternation stabilizes the highest occupied molecular orbital (HOMO) and destabilizes the lowest unoccupied molecular orbital (LUMO), in equivalence to Peierls distortion in infinite atomic chains.69 Normal bond length alternation opens the HOMOLUMO gap and consequently the gaps are larger for alkynes than for alkenes as plotted in Figure 2b. Reverse bond length alternation has the opposite effect and therefore the HOMO-LUMO gap for cumulenes is smaller than that of alkenes. Though alkynes and cumulenes have almost the same stucture, the calculated HOMO-LUMO gap is over half an eV smaller for the cumulenes.

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 47

In cumulenes a distortion of the bond lengths that stabilizes the HOMO will destabilize the neardegenerate HOMO–1 because its nodal structure is displaced by one atom.32 Furthermore, the HOMO–1 has significant hyperconjugation to the hydrogen atoms, thus making the normal bond length distortion more expensive. The antibonding hyperconjugative interaction means that the inplane and out-of-plane MOs switch place for the longer members of the series. The extra cost of the hyperconjugative interaction appears to be the reason the in-plane π system wins, and the bond length alternation is reversed as seen in Figure 1. Consequently, the long π system which we expect to dominate the electron transport has its bond lengths reversed, and correspond to the situation described by Tsuji et al.17

Model Calculations We now turn to model calculations to understand the relation between electronic structure and transmission decay with the Hückel method (tight binding model with only nearest neighbor interactions). We set the coulomb integral (on-site energy) to 0 eV, and for a double-bond we use resonance integral t = –3 eV. We use different variations of t based on the DFT-calculated bond lengths for the different molecules. We make three models. The first two are shown diagrammatically in Scheme 2, and are based on the bond lengths of the alkene series as described in Supporting Information part B. Here, for the formal double bonds we use td = t, and for the formal single-bonds we use ts = 0.935·t; we refer to this model as Normal-BLA because t d > t s. Scheme 2: Normal-BLA and Reverse-BLA models

10 ACS Paragon Plus Environment

Page 11 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

To understand the opposite limit, we reverse the model so td = 0.935·t and t s = t. This model we refer to as Reverse-BLA because ts > td, that is, the formal single-bonds are “shorter” than the formal double-bonds. However, from the data presented in Figure 1 we note that the cumulene bond lengths are not completely reversed compared with the alkene and alkyne. It is primarily the first bond that is long, after that there is only very weak bond-length alternation. Shown diagrammatically in Scheme 3, we make a model which corresponds to the bond lengths of the cumulene series, which we refer to as Cumulene-BLA. Here, there is only weak bond-length alternation in the central bonds. We use t1 = t and t2 = 0.990· t, while for the longer edge bonds we use ted = 0.968·t. This model corresponds only to the cumulene π-system that is out of hydrogen plane, i.e., the long π-system which we expect to dominate the transport properties. Scheme 3: Cumulene-BLA model.

We calculate the Landauer transmission with the non-equilibrium Green’s functions (NEGF) approach with wide-band electrodes attached to the terminal sites in our models with non-zero matrix elements of ΓL/R set to γ = 0.5 eV. In Figure 3 we examine the transmission decay. The transmission decays exponentially with length for the Normal-BLA model, and increases exponentially with length for the Reverse-BLA model as previously predicted17; the cumulene model is in-between as one would expect from equation 2. Comparing Figure 3a and 3b, there are two differences worth noting. First, the smaller HOMO-LUMO gap in Reverse-BLA moves the resonances closer to the Fermi energy, which is arbitrarily set to zero in the model. And more

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 47

subtly, the resonance shape is wider (less decaying) for Reverse-BLA. The resonance position and broadness are the two effects that control the off-resonant transmission.

Figure 3. Transmission calculated with Hückel method for different model lengths (values of n). (a) Transmission as function of energy for the normal-BLA model. (b) Transmission as function of energy for the reverse-BLA model. (c) Transmission at the Fermi energy for reverse-BLA (orange), cumulene-BLA(red), and normal-BLA models (purple). Zeroth-order Greens Function. For Hückel models the transmission at the Fermi energy, T(EF), can be approximated using the zeroth-order Green’s function as described in the work of Yoshizawa and coworkers 18, 70-71. This analysis yields the following proportionality.

12 ACS Paragon Plus Environment

Page 13 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

( ) ∝  

, ∙ ," $ # − 

%

(3)

Here ci,L/R is the orbital coefficient of the i'th molecular orbital on the site of contact of the left (L) or right (R) electrode, and εi is the eigenvalue of the i'th molecular orbital. The numerator and denominator are not independent of each other (the eigenvectors and eigenvalues of the Hamiltonian). The decay of transmission is essentially a competition between the orbital coefficients at the contact points becoming smaller, and the HOMO-LUMO gap becoming smaller (making the inverse of the molecular orbital eigenvalues larger). There is an important link to the diradical character of the molecules here, which in the MO picture loosely corresponds to frontier orbitals with large orbital coefficients on the terminal atoms and small HOMO-LUMO gap.72-75 The connection between diradical character and transmission has been discussed in detail elsewhere,29, 72-73, 76-79

as well as its relation to the coupling between the associated magnetic centers73, 80-84.

Following this type of analysis, Tada and Yoshizawa found that reverse transmission decay in modified graphene appears due to the frontier orbital localizing on the edges.18-19, 85-86 To explore this suggestion for a cumulenic system, we decompose the zeroth order Green’s function (eq. 3) for the HOMO of the three Hückel models. In Figure 4a the square of the term for the HOMO in equation 3 is plotted. It is, to a very good approximation, proportional to the transmission at the Fermi energy (Figure 3c), as expected by considering the Coulson-Rushbrooke pairing theorem.87-88 We now look at the two factors: The square of the inverse HOMO eigenvalue (panel b, mind that EF is mid-gap) and the square of the product of the orbital coefficients (panel c). While the inverse of the HOMO eigenvalue results in a large difference between the three models, there is only a small difference between the orbital coefficients of the three. In all three cases the terminal orbital coefficients still decay significantly with length, and reverse bond-alternation only dampens this effect.

13 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 47

Figure 4. Decomposition of the contribution from the HOMO to the transmission by analysis of the zeroth-order Green’s function at the Fermi energy as function of length. (a) The term of the HOMO in the zeroth-order Green’s function has the same trend as the transmission in Figure 3. (b) Square of the inverse of the HOMO eigenvalue, i.e., the denominator in the zeroth-order Green’s function. (c) Square of the product of the HOMO orbital coefficients at the left and right contact sites in, i.e., the numerator in the zeroth-order Green’s function. By direct comparison of Figure 4b and 4c, it is clear that the main contributor to the reverse decay of transmission is the smaller HOMO-LUMO gap.27 Algethami et al. also attributed the increasing transmission with length in fused porphyrin oligomers to a rapid decrease in the HOMO-LUMO gap as the oligomer length is increased.24 In this context a vanishing HOMO-LUMO gap generally

14 ACS Paragon Plus Environment

Page 15 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

corresponds to increasing diradical character, which may provide some chemical intuition for reverse conductance decay. It is clear that as the HOMO-LUMO gap vanishes the conductance will saturate at a value close to G0, the quantum of conductance, and the transport will become resonant; thus forming the upper limit for reverse conductance decay. In cumulenes the increasing diradical character with length can be observed in the reduced energetic cost of axial torsion as the groundstate and diradical transition-state geometries approach each other energetically.30,

89-92

The

importance of the terminal orbital coefficients and the HOMO-LUMO gap was also discussed in detail by Stuyver et al.29 on the topic of conductance decay, and Proppe and Herrmann73 who looked at electronic communication through molecular bridges in a broader context. It is clear that the increasing diradical character and vanishing HOMO-LUMO gap are intimately linked to the edge-localization of frontier orbitals described for functionalized graphene.18-19,

85-86

Further

exploration of this connection between the one- and two-dimensional cases may lead to additional insight into the physical mechanism of reverse transmission decay and to more candidate molecular systems with this intriguing property.

DFT Transmission Let us now examine the electron transport with DFT combined with the NEGF approach93. Single-molecule junctions are built by placing the molecule between two four-atom Au pyramids (tetrahedrons) on 4x4 Au(111) surfaces with periodic boundary conditions in the plane of the surfaces. The molecule in the junction is relaxed to a force threshold of 0.05 eV/Å with all Au atoms kept fixed during the optimization, and the transmission is calculated with a k-point sampling over a 4x4x1 Monkhorst-Pack mesh in the first Brilloun zone. Many studies have focused on the transport properties of atomic chains of carbon, but none of them correspond to the molecular geometries typically probed in STM break-junction experiments.41, 84, 94-98 Prasongkit et al.94 found the transmission of even-carbon cumulenic wires to

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 47

be almost independent of length, and similar trends were seen for the magnetic coupling through such systems by Sarbadhikary et al.84 However, the type of electrode-molecule contact considerably influences the transport properties of π-conjugated molecules.99 Here we construct Au-molecule-Au junctions where the molecule is weakly bound through a donor-acceptor bond between amine groups and Au electrodes, where the contact geometries are generally well-defined and the transport is known to be off-resonant.100 In Figure 5a, optimized junction structures are shown for the eightcarbon species.

Figure 5. DFT transmission of amine-linked single-molecule junctions. (a) Optimized junction structures of amine-linked [7]alkyne, [7]cumulene, and [7]alkene. (b) Transmission plotted semilogarithmically against energy for odd [n]cumulenes, n = 3-9. c) Transmission plotted semilogarithmically against energy for [n]alkynes, n = 3-9. (d) Transmission at the Fermi energy as function of length (n) for [n]alkene (purple), [n]cumulene (red), and [n]alkyne (blue), for n = 3-19. All transmission plots are included in SI part C. The transmission of the four shortest members of the alkyne and cumulene series are plotted in Figure 5b and 5c. As expected from the model calculations the transmission decays exponentially in the middle of the HOMO-LUMO gap (around EF) for the alkyne series, while the transmission does not show clear decay for the cumulene series as the HOMO-LUMO rapidly becomes narrow with

16 ACS Paragon Plus Environment

Page 17 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

increasing length. Plotting the transmission at the Fermi energy against length (n) for the three series in Figure 5d, we see the trend we expect from the bond-length alternation (Figure 1); the alkyne a has larger decay constant, β, than the alkene series, while the cumulene series has a slightly negative β, i.e., reverse decay of transmission with length. It is important to note that β is subject to the energy level alignment relative to the Fermi energy which will depend on the anchor group used, and for calculations it will be method dependent. We note that there is reasonable agreement between the alkyne calculations and reported experimental values of β.101-107 Given that we underestimate the reverse bond-length alternation with the method used, in particular for the longer members of the cumulene series (see Fig. S1 and S2), we expect the result to persist at higher-levels of theory.

Increasing Reverse Bond-length Alternation in Cumulenes The bond-length alternation and conductance decay are reversed for simple functionalized cumulenes; however, these effects are fairly weak. Therefore, it is imperative to explore ways to increase these effects. The electronic structure can be manipulated chemically by using electron withdrawing and donating substituents. If donors and acceptors are placed at opposite ends of the molecule a push-pull system is formed and the charge distribution is manipulated. Consequently, resonance structures with zwitterionic character will become increasingly significant. This strategy was recently used to manipulate the cross-conjugated character of oligo(phenyleneethynylene) wires.108 The push-pull strategy has recently also successfully been applied to [3]cumulenes resulting in increasingly reversed bond-length alternation as the zwitterionic resonance structure becomes more pronounced than the cumulenic one; the two contributors shown in the top panel of Scheme 4 represent a simple (and perhaps less realistic) example of such a molecule.109-113 These

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 47

resonance structures provide an instructive picture as to how the redistribution of charge in the molecule directly affects the bond-length alternation, as a reverse single-triple bond alternation becomes predominant in the push-pull substituted species. Scheme 4: Cumulenic and Zwitterionic Resonance Contributors of Substituted Cumulenes.

There are several substituent strategies that achieve this effect. Replacing the entire end-groups by different odd-carbon rings, forming a pentaheptafulvalene-expanded cumulene will do exactly this114-115; such species were previously found to have a substantial dipole moment as the electron structure adapts towards 6 electrons in each ring, as shown in Scheme 4 (middle).116-117 A number of substituted cumulenes show increased reverse bond-length alternation as the substituents may help stabilize different resonance contributors, such as the pentafulvalene-expanded cumulene shown at the bottom of Scheme 4.37, 50, 118-121 For this purpose a captodative strategy was explored theoretically for stabilizing quinoidal resonance structures which give rise to reverse conductance decay in recent work by Stuyver et al.122

18 ACS Paragon Plus Environment

Page 19 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 6. DFT-calculated bond lengths of substituted [n]cumulenes. (a) Chemical structures of [n]cumulene,

(CN)2[n]cumulene(NH2)2,

[n]cumulene-extended

pentaheptafulvalene,

and

[n]cumulene-extended pentafulvalene, where is x is a positive integer and n = 2x+1. (b) Optimized structures and bond lengths in Ångström for the four [3]cumulene species. Note the fulvalene structures are scaled smaller due to their size. (c) Bond lengths of the four [11]cumulene species. The calculated structures for the substituted cumulenes shown in scheme 4 and Figure 6 all show increased reverse bond-length alternation. In all three substituted [3]cumulenes the edge bonds are longer and the central bond is shorter compared to unsubstituted [3]cumulene. A trend that persists for the longer members of each series. The pentafulvalene expanded cumulene shows slightly weaker bond-length alternation than the other systems, but still considerably stronger than the unsubstituted cumulene. These values are in good agreement with crystallographic data of similar systems.109-110, 113, 120, 123

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 47

Figure 7. DFT transmission of amine-linked single-molecule junctions. (a) Optimized junction structures of amine-linked [5]pentaheptafulvalene, [5]pentafulvalene, and [7]cumulene. (b) Transmission plotted semilogarithmically against energy for odd [n]pentaheptafulvalene, n = 3-9. c) Transmission plotted semilogarithmically against energy for odd [n]pentafulvalene, n = 3-9. (d) Transmission at the Fermi energy as function of length (n) for [n]cumulene (red), [n]pentaheptafulvalene (blue), and [n]pentafulvalene (purple), for n = 3-15. By increasing the reverse bond length alternation, the reverse conductance decay trend is increased. In Figure 7, the transmission of the penteheptafulvalene- and pentafulvalene-expanded cumulenes are calculated with DFT as in the previous section. The (CN)2–(NH2)2 substituted system is not included as we cannot make junctions with comparable anchor groups. The transmission trend of the expanded cumulenes is dependent on the energy level alignment relative to EF. Still, the trend we expect from the model calculations is clear. Stronger reverse bondlength alternation results in stronger reverse decay of conductance with length. We see this in the more negative β for the pentaheptafulvalene and pentafulvalene series in Figure 7d. The conductance decay in push-pull systems is likely to be bias and solvent dependent as the molecules are likely to show rectification.124-126 The electric field that the bias induces will, pending the orientation of the molecule, further enhance the push-pull effect on the bond-length alternation.

20 ACS Paragon Plus Environment

Page 21 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

However, it must be noted that the effect of the electric field on the bond-length alternation is predicted to be modest at low bias.127-129

Prospective Systems and Outlook It is clear that stability of the longer members of a series will be one of the main challenges for realizing reverse decay on conductance. For cumulenes this in an inherent challenge as the longer members are less stable and cannot be protected by normal strategies such as protection through sterically crowding the molecule. Ingenious strategies for protecting long linear carbons molecules are being undertaken; polyynes and cumulenes protected with rotaxane rings have been reported in recent years.49, 130-133 Multipodal platform approaches holds promise for systematic control of the contact geometry of linear one-dimensional wires, which may help avoiding dimerization by spatial separation of the wires.134-137 The general stability challenge for any molecular series showing reverse conductance decay is the increased diradical character that seem to be its prerequisite.29 Push-pull and captodative substituent strategies are promising ways to accommodate (or reduce the) diradical character.122 It remains the central challenge to accommodate molecules with increasingly small HOMO-LUMO gaps which is a common denominator for highly conducting molecules, particularly those with reverse bondlength alternation. This is in analogy to the suggested negative correlation between conductance and aromaticity in small cyclic systems.138-142 In the aromatic limit the molecules are stable with a big HOMO-LUMO gap, in the anti-aromatic limit the molecules are unstable diradicals with vanishing HOMO-LUMO gaps. Scheme 5: Quinodal and Aromatic Resonance Contributors of the Quinodimethane Unit.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 47

Despite the afore-mentioned stability challenges, recent synthetic efforts have successfully seen quinodimethane units fused into cumulenic chains. Oligo-quinodimethanes themselves show strong reverse decay of transmission as the aromatic sextet of each phenyl ring is an important resonance contributor in addition to the quinoidal structure as shown in Scheme 5 (top).29,

143

As length

increases the quinodal structure loses out and diradical character increases, therefore only few oligo-quinodimethane derivates have been successfully isolated.144-145 The aromatic sextets promote reverse bond length alternation when incorporated into cumulenes as shown in Scheme 5. Recent synthetic reports of quinoid-expanded [3]cumulenes all display strong reverse bond alternation in the crystal structures; in all cases terminal substituent groups that can accommodate diradical character were used.146-149 It is clear that it will be a challenge to make long oligomers of such systems, still these are prospective systems for reverse conductance decay. Transmission calculations on such systems, included in Supporting Information part D, show clear reverse decay of transmission with length. Finally, we should mention a class of cumulenes we have not yet touched upon in this paper. Research in metallocumulenes (cumulenylidene complexes) has seen great progress and may hold promise for making longer cumulenic oligomers.150-152 From the perspective of the cumulene these show reverse bond-length alternation;153-158 some approach fully reversed alkyne bond lengths.159 As seen in polyynediyl complexes, metal centers may allow for stabilization of longer chains.160 Stability can further be enhanced by using large extended ligands, which can span the full length of the polyyne wire.161-162 The extra protection has already allowed the electron transport properties of long metallopolyynes to be probed experimentally.163-165 The role of one or more metal centers in

22 ACS Paragon Plus Environment

Page 23 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

the wire adds an extra level of complexity and the transport properties of metallocumulenes will be examined in future studies.166-167

Conclusions We have reexamined the peculiar result of Tsuji et al.17 that linear conjugated molecules with reverse bond length alternation show increasing tunneling transmission as the length of the molecule is increased. We have examined how cumulenes are a class of molecule that show this behavior. Despite their all-double bond structure, cumulenic bonds are not of equal length and the bond lengths are reversed for the dominant π system that span the full length of the carbon backbone. Consequently, we predict a negative decay constant, β, for several series of cumulenes. The transmission increases as the HOMO-LUMO gap rapidly becomes narrower with length; a direct consequence of reverse bond-length alternation. Reverse bond-length alternation in cumulenes further increases if suitable substituent strategies are applied. Based on recent synthetic progress in extended cumulenes, we find there is a great number of potential cumulene series that may show reverse decay of conductance with length. Moving forward, the central challenge in realizing reverse conductance decay with length will be stabilizing the longer members of the series. Corresponding Author *E-mail: [email protected] Supporting Information The Supporting Information is available free of charge on the ACS Publications website. Benchmark study of bond-length alternation; Hückel models, Transmission Plots; Cumuleneextended quinodimethane.

ACKNOWLEDGMENTS

23 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 47

We are grateful for funding from the Danish Council for Independent Research|Natural Sciences and the Carlsberg Foundation. References 1.

Wold, D. J.; Haag, R.; Rampi, M. A.; Frisbie, C. D. Distance Dependence of Electron

Tunneling Through Self-Assembled Monolayers Measured by Conducting Probe Atomic Force Microscopy:  Unsaturated Versus Saturated Molecular Junctions. J. Phys. Chem. B 2002, 106, 2813-2816. 2.

Chen, F.; Tao, N. J. Electron Transport in Single Molecules: From Benzene to Graphene.

Acc. Chem. Res. 2009, 42, 429-438. 3.

Wenger, O. S. How Donor−Bridge−Acceptor Energetics Influence Electron Tunneling

Dynamics and Their Distance Dependences. Acc. Chem. Res. 2011, 44, 25-35. 4.

Kirk, M. L.; Shultz, D. A.; Stasiw, D. E.; Lewis, G. F.; Wang, G.; Brannen, C. L.; Sommer,

R. D.; Boyle, P. D. Superexchange Contributions to Distance Dependence of Electron Transfer/Transport: Exchange and Electronic Coupling in Oligo(Para-Phenylene)- and Oligo(2,5Thiophene)-Bridged Donor–Bridge–Acceptor Biradical Complexes. J. Am. Chem. Soc. 2013, 135, 17144-17154. 5.

Winkler, J. R.; Gray, H. B. Long-Range Electron Tunneling. J. Am. Chem. Soc. 2014, 136,

2930-2939. 6.

Khoo, K. H.; Chen, Y.; Li, S.; Quek, S. Y. Length Dependence of Electron Transport

Through Molecular Wires - a First Principles Perspective. Phys. Chem. Chem. Phys. 2015, 17, 7796.

24 ACS Paragon Plus Environment

Page 25 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

7.

Schubert, C.; Margraf, J. T.; Clark, T.; Guldi, D. M. Molecular Wires - Impact of Pi-

Conjugation and Implementation of Molecular Bottlenecks. Chem. Soc. Rev. 2015, 44, 988-998. 8.

Xu, B.; Tao, N. J. Measurement of Single-Molecule Resistance by Repeated Formation of

Molecular Junctions. Science 2003, 301, 1221. 9.

Venkataraman, L.; Klare, J. E.; Tam, I. W.; Nuckolls, C.; Hybertsen, M. S.; Steigerwald, M.

L. Single-Molecule Circuits with Well-Defined Molecular Conductance. Nano Lett. 2006, 6, 458462. 10. Kaliginedi, V.; Moreno-García, P.; Valkenier, H.; Hong, W.; García-Suárez, V. M.; Buiter, P.; Otten, J. L. H.; Hummelen, J. C.; Lambert, C. J.; Wandlowski, T. Correlations between Molecular Structure and Single-Junction Conductance: A Case Study with Oligo(PhenyleneEthynylene)-Type Wires. J. Am. Chem. Soc. 2012, 134, 5262-5275. 11. McConnell, H. M. Intramolecular Charge Transfer in Aromatic Free Radicals. J. Chem. Phys. 1961, 35, 508-515. 12. Reimers, J. R.; Hush, N. S. Analytic Solutions to Resonant and Non-Resonant ThroughBridge Electronic Coupling. Nanotechnology 1996, 7, 417. 13. Magoga, M.; Joachim, C. Conductance and Transparence of Long Molecular Wires. Phys. Rev. B 1997, 56, 4722-4729. 14. Kushmerick, J. G.; Holt, D. B.; Pollack, S. K.; Ratner, M. A.; Yang, J. C.; Schull, T. L.; Naciri, J.; Moore, M. H.; Shashidhar, R. Effect of Bond-Length Alternation in Molecular Wires. J. Am. Chem. Soc. 2002, 124, 10654-10655.

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 47

15. Nitzan, A. Electron Transmission Through Molecules Aan Molecular Interfaces. Annu. Rev. Phys. Chem. 2001, 52, 681-750. 16. Hsu, L.-Y.; Rabitz, H. Theory of Molecular Conductance Using a Modular Approach. J. Chem. Phys. 2016, 145, 234702. 17. Tsuji, Y.; Movassagh, R.; Datta, S.; Hoffmann, R. Exponential Attenuation of ThroughBond Transmission in a Polyene: Theory and Potential Realizations. ACS Nano 2015, 9, 1110911120. 18. Tada, T.; Yoshizawa, K. Reverse Exponential Decay of Electrical Transmission in Nanosized Graphite Sheets. J. Phys. Chem. B 2004, 108, 7565-7572. 19. Li, S.; Gan, C. K.; Son, Y.-W.; Feng, Y. P.; Quek, S. Y. Anomalous Length-Independent Frontier Resonant Transmission Peaks in Armchair Graphene Nanoribbon Molecular Wires. Carbon 2014, 76, 285-291. 20. Sedghi, G.; Sawada, K.; Esdaile, L. J.; Hoffmann, M.; Anderson, H. L.; Bethell, D.; Haiss, W.; Higgins, S. J.; Nichols, R. J. Single Molecule Conductance of Porphyrin Wires with Ultralow Attenuation. J. Am. Chem. Soc. 2008, 130, 8582-8583. 21. Sedghi, G.; García-Suárez, V. M.; Esdaile, L. J.; Anderson, H. L.; Lambert, C. J.; Martín, S.; Bethell, D.; Higgins, S. J.; Elliott, M.; Bennett, N. et al., Long-Range Electron Tunnelling in OligoPorphyrin Molecular Wires. Nat. Nanotech. 2011, 6, 517. 22. Li, Z.; Park, T.-H.; Rawson, J.; Therien, M. J.; Borguet, E. Quasi-Ohmic Single Molecule Charge Transport Through Highly Conjugated Meso-to-Meso Ethyne-Bridged Porphyrin Wires. Nano Lett. 2012, 12, 2722-2727.

26 ACS Paragon Plus Environment

Page 27 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

23. Sedghi, G.; Esdaile, L. J.; Anderson, H. L.; Martin, S.; Bethell, D.; Higgins, S. J.; Nichols, R. J. Comparison of the Conductance of Three Types of Porphyrin-Based Molecular Wires: Β,Meso,Β-Fused Tapes, Meso-Butadiyne-Linked and Twisted Meso-Meso Linked Oligomers. Adv. Mater. 2012, 24, 653-657. 24. Algethami, N.; Sadeghi, H.; Sangtarash, S.; Lambert, C. J. The Conductance of PorphyrinBased Molecular Nanowires Increases with Length. Nano Lett. 2018. 25. Ramos-Berdullas, N.; Mandado, M. Electronic Properties of P-Xylylene and P-Phenylene Chains Subjected to Finite Bias Voltages: A New Highly Conducting Oligophenyl Structure. Chem. Eur. J. 2013, 19, 3646-3654. 26. Stuyver, T.; Fias, S.; De Proft, F.; Geerlings, P. The Relation Between Delocalization, Long Bond Order Structure Count and Transmission: An Application to Molecular Wires. Chem. Phys. Lett. 2015, 630, 51-56. 27. Ramos‐Berdullas, N.; Gil‐Guerrero, S.; Mandado, M. Transmission Channels in the Time‐Energy Uncertainty Relation Approach to Molecular Conductance: Symmetry Rules for the Electron Transport in Molecules. Int. J. Quantum Chem. 2018, 0, e25651. 28. Tada, T.; Yoshizawa, K. Molecular Design of Electron Transport with Orbital Rule: Toward Conductance-Decay Free Molecular Junctions. Phys. Chem. Chem. Phys. 2015, 17, 32099-32110. 29. Stuyver, T.; Fias, S.; Proft, F. D.; Geerlings, P.; Tsuji, Y.; Hoffmann, R. Enhancing the Conductivity of Molecular Electronic Devices. J. Chem. Phys. 2017, 146, 092310. 30. Hoffmann, R. Extended Hückel Theory—V: Cumulenes, Polyenes, Polyacetylenes and Cn. Tetrahedron 1966, 22, 521-538.

27 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 47

31. Mölder, U.; Burk, P.; Koppel, I. A. Quantum Chemical Calculations of Linear Cumulene Chains. J. of Mol. Struct.: THEOCHEM 2004, 712, 81-89. 32. Imamura, A.; Aoki, Y. Molecular Design of a Π-Conjugated Single-Chain Electronically Conductive Polymer. Int. J. Quantum Chem. 2006, 106, 1924-1933. 33. Innocenti, F.; Milani, A.; Castiglioni, C. Can Raman Spectroscopy Detect Cumulenic Structures of Linear Carbon Chains? J. Raman Spectrosc. 2010, 41, 226-236. 34. Cahangirov, S.; Topsakal, M.; Ciraci, S. Long-Range Interactions in Carbon Atomic Chains. Phys. Rev. B 2010, 82, 195444. 35. Imamura, A.; Aoki, Y. Electronic Structures and Molecular Structures of Polyynes. Int. J. Quantum Chem. 2013, 113, 423-427. 36. Neiss, C.; Trushin, E.; Görling, A. The Nature of One-Dimensional Carbon: Polyynic Versus Cumulenic. ChemPhysChem 2014, 15, 2497-2502. 37. Tommasini, M.; Milani, A.; Fazzi, D.; Lucotti, A.; Castiglioni, C.; Januszewski, J. A.; Wendinger, D.; Tykwinski, R. R. Π-Conjugation and End Group Effects in Long Cumulenes: Raman Spectroscopy and Dft Calculations. J. Phys. Chem. C 2014, 118, 26415-26425. 38. Li, Y.; Mondal, K. C.; Samuel, P. P.; Zhu, H.; Orben, C. M.; Panneerselvam, S.; Dittrich, B.; Schwederski, B.; Kaim, W.; Mondal, T., et al. C4 Cumulene and the Corresponding Air‐Stable Radical Cation and Dication. Angew. Chem. Int. Ed. 2014, 53, 4168-4172. 39. Casari, C. S.; Tommasini, M.; Tykwinski, R. R.; Milani, A. Carbon-Atom Wires: 1-D Systems with Tunable Properties. Nanoscale 2016, 8, 4414-4435.

28 ACS Paragon Plus Environment

Page 29 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

40. Milani, A.; Tommasini, M.; Barbieri, V.; Lucotti, A.; Russo, V.; Cataldo, F.; Casari, C. S. Semiconductor-to-Metal Transition in Carbon-Atom Wires Driven by Sp2 Conjugated End Groups. J. Phys. Chem. C 2017, 121, 10562-10570. 41. Lambropoulos, K.; Simserides, C. Electronic Structure and Charge Transport Properties of Atomic Carbon Wires. Phys. Chem. Chem. Phys. 2017, 19, 26890-26897. 42. Larsen, A. H.; Mortensen, J. J.; Blomqvist, J.; Castelli, I. E.; Christensen, R.; Dułak, M.; Friis, J.; Groves, M. N.; Hammer, B.; Hargus, C. et al., The Atomic Simulation Environment—a Python Library for Working with Atoms. J. Phys.: Condens. Matter 2017, 29, 273002. 43. Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. 44. Mortensen, J. J.; Hansen, L. B.; Jacobsen, K. W. Real-Space Grid Implementation of the Projector Augmented Wave Method. Phys. Rev. B 2005, 71, 035109. 45. Larsen, A. H.; Vanin, M.; Mortensen, J. J.; Thygesen, K. S.; Jacobsen, K. W. Localized Atomic Basis Set in the Projector Augmented Wave Method. Phys. Rev. B 2009, 80, 195112. 46. Leroyer, L.; Maraval, V.; Chauvin, R. Synthesis of the Butatriene C4 Function: Methodology and Applications. Chem. Rev. 2012, 112, 1310-1343. 47. Januszewski, J. A.; Wendinger, D.; Methfessel, C. D.; Hampel, F.; Tykwinski, R. R. Synthesis and Structure of Tetraarylcumulenes: Characterization of Bond-Length Alternation Versus Molecule Length. Angew. Chem. Int. Ed. 2013, 52, 1817-1821. 48. Januszewski, J. A.; Tykwinski, R. R. Synthesis and Properties of Long [N]Cumulenes (N ≥ 5). Chem. Soc. Rev. 2014, 43, 3184-3203.

29 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 47

49. Franz, M.; Januszewski, J. A.; Wendinger, D.; Neiss, C.; Movsisyan, L. D.; Hampel, F.; Anderson, H. L.; Görling, A.; Tykwinski, R. R. Cumulene Rotaxanes: Stabilization and Study of [9]Cumulenes. Angew. Chem. Int. Ed. 2015, 54, 6645-6649. 50. Wendinger, D.; Tykwinski, R. R. Odd [N]Cumulenes (N = 3, 5, 7, 9): Synthesis, Characterization, and Reactivity. Acc. Chem. Res. 2017, 50, 1468-1479. 51. Sun, Q.; Tran, B. V.; Cai, L.; Ma, H.; Yu, X.; Yuan, C.; Stöhr, M.; Xu, W. On‐Surface Formation of Cumulene by Dehalogenative Homocoupling of Alkenyl Gem‐Dibromides. Angew. Chem. Int. Ed. 2017, 56, 12165-12169. 52. Marshall, J. L.; Lehnherr, D.; Lindner, B. D.; Tykwinski, R. R. Reductive Aromatization/Dearomatization and Elimination Reactions to Access Conjugated Polycyclic Hydrocarbons, Heteroacenes, and Cumulenes. ChemPlusChem 2017, 82, 967-1001. 53. Hendon, C. H.; Tiana, D.; Murray, A. T.; Carbery, D. R.; Walsh, A. Helical Frontier Orbitals of Conjugated Linear Molecules. Chem. Sci. 2013, 4, 4278-4284. 54. Imamura, A.; Aoki, Y. Helical Molecular Orbitals Around Straight-Chain Polyyne Oligomers as Models for Molecular Devices. Chem. Phys. Lett. 2013, 590, 136-140. 55. Garner, M. H.; Hoffmann, R.; Rettrup, S.; Solomon, G. C. Coarctate and Möbius: The Helical Orbitals of Allene and Other Cumulenes. ACS Cent. Sci. 2018, 4, 688-700. 56. Szafert, S.; Gladysz, J. A. Carbon in One Dimension:  Structural Analysis of the Higher Conjugated Polyynes. Chemical Reviews 2003, 103, 4175-4206. 57. Szafert, S.; Gladysz, J. A. Update 1 Of: Carbon in One Dimension:  Structural Analysis of the Higher Conjugated Polyynes. Chem. Rev. 2006, 106, PR1-PR33.

30 ACS Paragon Plus Environment

Page 31 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

58. Chalifoux, W. A.; Tykwinski, R. R. Synthesis of Polyynes to Model the Sp-Carbon Allotrope Carbyne. Nat. Chem. 2010, 2, 967-971. 59. Tykwinski, R. R. Carbyne: The Molecular Approach. Chem. Rec. 2015, 15, 1060-1074. 60. Jin, C.; Lan, H.; Peng, L.; Suenaga, K.; Iijima, S. Deriving Carbon Atomic Chains from Graphene. Phys. Rev. Lett. 2009, 102, 205501. 61. Liu, M.; Artyukhov, V. I.; Lee, H.; Xu, F.; Yakobson, B. I. Carbyne from First Principles: Chain of C Atoms, a Nanorod or a Nanorope. ACS Nano 2013, 7, 10075-10082. 62. Artyukhov, V. I.; Liu, M.; Yakobson, B. I. Mechanically Induced Metal–Insulator Transition in Carbyne. Nano Lett. 2014, 14, 4224-4229. 63. La Torre, A.; Botello-Mendez, A.; Baaziz, W.; Charlier, J. C.; Banhart, F. Strain-Induced Metal–Semiconductor Transition Observed in Atomic Carbon Chains. Nat. Comm. 2015, 6, 6636. 64. Ben Romdhane, F.; Adjizian, J.-J.; Charlier, J.-C.; Banhart, F. Electrical Transport Through Atomic Carbon Chains: The Role of Contacts. Carbon 2017, 122, 92-97. 65. Shi, L.; Rohringer, P.; Suenaga, K.; Niimi, Y.; Kotakoski, J.; Meyer, J. C.; Peterlik, H.; Wanko, M.; Cahangirov, S.; Rubio, A. et al., Confined Linear Carbon Chains as a Route to Bulk Carbyne. Nat. Mater. 2016, 15, 634-639. 66. Yang, S.; Kertesz, M. Linear Cn Clusters:  Are They Acetylenic or Cumulenic? J. Phys. Chem. A 2008, 112, 146-151. 67. Fan, X.; Liu, L.; Lin, J.; Shen, Z.; Kuo, J.-L. Density Functional Theory Study of Finite Carbon Chains. ACS Nano 2009, 3, 3788-3794.

31 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 47

68. Casari, C. S.; Milani, A. Carbyne: From the Elusive Allotrope to Stable Carbon Atom Wires. MRS Commun. 2018, 8, 207-219. 69. Aoki, Y.; Imamura, A. An Analytical Hückel‐Type Approach to the Relationship Between Peierls Instability in Polyenes and Interchain Interaction. J. Chem. Phys. 1995, 103, 9726-9737. 70. Tada, T.; Yoshizawa, K. Quantum Transport Effects in Nanosized Graphite Sheets. ChemPhysChem 2002, 3, 1035-1037. 71. Yoshizawa, K.; Tada, T.; Staykov, A. Orbital Views of the Electron Transport in Molecular Devices. J. Am. Chem. Soc. 2008, 130, 9406-9413. 72. Herrmann, C.; Elmisz, J. Electronic Communication Through Molecular Bridges. Chem. Comm. 2013, 49, 10456-10458. 73. Proppe, J.; Herrmann, C. Communication Through Molecular Bridges: Different Bridge Orbital Trends Result in Common Property Trends. J. Comp. Chem. 2015, 36, 201-209. 74. Gu, J.; Wu, W.; Danovich, D.; Hoffmann, R.; Tsuji, Y.; Shaik, S. Valence Bond Theory Reveals Hidden Delocalized Diradical Character of Polyenes. J. Am. Chem. Soc. 2017, 139, 93029316. 75. Voigt, B. A.; Steenbock, T.; Herrmann, C. Structural Diradical Character. arXiv:1802.06695 [physics.comp-ph], 2018. 76. Stuyver, T.; Fias, S.; De Proft, F.; Geerlings, P. Back of the Envelope Selection Rule for Molecular Transmission: A Curly Arrow Approach. J. Phys. Chem. C 2015, 119, 26390-26400.

32 ACS Paragon Plus Environment

Page 33 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

77. Tsuji, Y.; Hoffmann, R.; Strange, M.; Solomon, G. C. Close Relation Between Quantum Interference in Molecular Conductance and Diradical Existence. Proc. Natl. Acad. Sci. 2016, 113, E413-E419. 78. Stuyver, T.; Blotwijk, N.; Fias, S.; Paul, G.; De Proft, F. Exploring Electrical Currents Through Nanographenes: Visualization and Tuning of the Through‐Bond Transmission Paths. ChemPhysChem 2017, 18, 3012-3022. 79. Tsuji, Y.; Estrada, E.; Movassagh, R.; Hoffmann, R. Quantum Interference, Graphs, Walks, and Polynomials. Chem. Rev. 2018. 80. Higashiguchi, K.; Yumoto, K.; Matsuda, K. Evaluation of the Β Value of the Phenylene Unit by Probing Exchange Interaction Between Two Nitroxides. Org. Lett. 2010, 12, 5284-5286. 81. Shinomiya, M.; Higashiguchi, K.; Matsuda, K. Evaluation of the Β Value of the Phenylene Ethynylene Unit by Probing the Exchange Interaction Between Two Nitronyl Nitroxides. J. Org. Chem. 2013, 78, 9282-9290. 82. Nishizawa, S.; Hasegawa, J.-y.; Matsuda, K., Theoretical Investigation of the Β Value of the Π-Conjugated Molecular Wires by Evaluating Exchange Interaction Between Organic Radicals. J. Phys. Chem. C 2013, 117, 26280-26286. 83. Sarbadhikary, P.; Shil, S.; Panda, A.; Misra, A. A Perspective on Designing Chiral Organic Magnetic Molecules with Unusual Behavior in Magnetic Exchange Coupling. J. Org. Chem. 2016, 81, 5623-5630. 84. Sarbadhikary, P.; Shil, S.; Misra, A. Magnetic and Transport Properties of Conjugated and Cumulated Molecules: The [Small Pi]-System Enlightens Part of the Story. Phys. Chem. Chem. Phys. 2018, 20, 9364-9375.

33 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 47

85. Zdetsis, A. D.; Economou, E. N. Interrelation of Aromaticity and Conductivity of Graphene Dots/Antidots and Related Nanostructures. J. Phys. Chem. C 2016, 120, 29463-29475. 86. Martín-Martínez, F. J.; Fias, S.; Hajgató, B.; Van Lier, G.; De Proft, F.; Geerlings, P. Inducing Aromaticity Patterns and Tuning the Electronic Transport of Zigzag Graphene Nanoribbons Via Edge Design. J. Phys. Chem. C 2013, 117, 26371-26384. 87. Coulson, C. A.; Rushbrooke, G. S. Note on the Method of Molecular Orbitals. Math. Proc. Camb. Philos. Soc. 1940, 36, 193-200. 88. Zhao, X.; Geskin, V.; Stadler, R. Destructive Quantum Interference in Electron Transport: A Reconciliation of the Molecular Orbital and the Atomic Orbital Perspective. J. Chem. Phys. 2016, 146, 092308. 89. Kuhn, R.; Schulz, B.; Jochims, J. C. Cis-Trans-Isomerism of 1,6-Diphenyl-1,6-Di-TButylhexapentaene. Angew. Chem. Int. Ed. 1966, 5, 420-420. 90. Auffrant, A.; Jaun, B.; Jarowski Peter, D.; Houk Kendall, N.; Diederich, F. Peralkynylated Buta-1,2,3-Trienes: Exceptionally Low Rotational Barriers of Cumulenic C=C Bonds in the Range of Those of Peptide Cn Bonds. Chem. Eur. J. 2004, 10, 2906-2911. 91. Jarowski, P. D.; Diederich, F.; Houk, K. N. Butatrienes as Extended Alkenes:  Barriers to Internal Rotation and Substitution Effects on the Stabilities of the Ground States and Transition States. J. Phys. Chem. A 2006, 110, 7237-7246. 92. Buehringer, M. U.; Padberg, K.; Phleps, M.; Maid, H.; Placht, C.; Neiss, C.; Ferguson, M.; Goerling, A.; Tykwinski, R. R. Double Bonds? Studies on the Barrier to Rotation About the Cumulenic C=C Bonds of Tetraaryl[N]Cumulenes (N = 3, 5, 7, 9). Angew. Chem. Int. Ed. 2018, 57, 8321-8325.

34 ACS Paragon Plus Environment

Page 35 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

93. Chen, J.; Thygesen, K. S.; Jacobsen, K. W. Ab Initio Nonequilibrium Quantum Transport and Forces with the Real-Space Projector Augmented Wave Method. Phys. Rev. B 2012, 85, 155140. 94. Prasongkit, J.; Grigoriev, A.; Wendin, G.; Ahuja, R. Cumulene Molecular Wire Conductance from First Principles. Phys. Rev. B 2010, 81, 115404. 95. Shen, L.; Zeng, M.; Yang, S.-W.; Zhang, C.; Wang, X.; Feng, Y. Electron Transport Properties of Atomic Carbon Nanowires Between Graphene Electrodes. J. Am. Chem. Soc. 2010, 132, 11481-11486. 96. Cretu, O.; Botello-Mendez, A. R.; Janowska, I.; Pham-Huu, C.; Charlier, J.-C.; Banhart, F. Electrical Transport Measured in Atomic Carbon Chains. Nano Lett. 2013, 13, 3487-3493. 97. D’yachkov, P. N.; Zaluev, V. A.; Kocherga, E. Y.; Sadykov, N. R. Tight Binding Model of Quantum Conductance of Cumulenic and Polyynic Carbynes. J. Phys. Chem. C 2013, 117, 1630616315. 98. Bonardi, P.; Achilli, S.; Tantardini, G. F.; Martinazzo, R. Electron Transport in Carbon Wires in Contact with Ag Electrodes: A Detailed First Principles Investigation. Phys. Chem. Chem. Phys. 2015, 17, 18413-18425. 99. Chen, W.; Widawsky, J. R.; Vázquez, H.; Schneebeli, S. T.; Hybertsen, M. S.; Breslow, R.; Venkataraman, L. Highly Conducting Π-Conjugated Molecular Junctions Covalently Bonded to Gold Electrodes. J. Am. Chem. Soc. 2011, 133, 17160-17163. 100. Hybertsen, M. S.; Venkataraman, L.; Klare, J. E.; Whalley, A. C.; Steigerwald, M. L.; Nuckolls, C. Amine-Linked Single-Molecule Circuits: Systematic Trends across Molecular Families. J. Phys.: Condens. Matter 2008, 20, 374115.

35 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 36 of 47

101. He, J.; Chen, F.; Li, J.; Sankey, O. F.; Terazono, Y.; Herrero, C.; Gust, D.; Moore, T. A.; Moore, A. L.; Lindsay, S. M. Electronic Decay Constant of Carotenoid Polyenes from SingleMolecule Measurements. J. Am. Chem. Soc. 2005, 127, 1384-1385. 102. Visoly-Fisher, I.; Daie, K.; Terazono, Y.; Herrero, C.; Fungo, F.; Otero, L.; Durantini, E.; Silber, J. J.; Sereno, L.; Gust, D. et al., Conductance of a Biomolecular Wire. Proc. Natl. Acad. Sci. 2006, 103, 8686-8690. 103. Wang, C.; Batsanov, A. S.; Bryce, M. R.; Martín, S.; Nichols, R. J.; Higgins, S. J.; GarcíaSuárez, V. M.; Lambert, C. J. Oligoyne Single Molecule Wires. J. Am. Chem. Soc. 2009, 131, 15647-15654. 104. Meisner, J. S.; Kamenetska, M.; Krikorian, M.; Steigerwald, M. L.; Venkataraman, L.; Nuckolls, C. A Single-Molecule Potentiometer. Nano Lett. 2011, 11, 1575-1579. 105. Moreno-García, P.; Gulcur, M.; Manrique, D. Z.; Pope, T.; Hong, W.; Kaliginedi, V.; Huang, C.; Batsanov, A. S.; Bryce, M. R.; Lambert, C. et al., Single-Molecule Conductance of Functionalized Oligoynes: Length Dependence and Junction Evolution. J. Am. Chem. Soc. 2013, 135, 12228-12240. 106. Gulcur, M.; Moreno‐García, P.; Zhao, X.; Baghernejad, M.; Batsanov, A. S.; Hong, W.; Bryce, M. R.; Wandlowski, T. The Synthesis of Functionalised Diaryltetraynes and Their Transport Properties in Single‐Molecule Junctions. Chem. Eur. J. 2014, 20, 4653-4660. 107. Milan, D. C.; Al-Owaedi, O. A.; Oerthel, M.-C.; Marqués-González, S.; Brooke, R. J.; Bryce, M. R.; Cea, P.; Ferrer, J.; Higgins, S. J.; Lambert, C. J. et al., Solvent Dependence of the Single Molecule Conductance of Oligoyne-Based Molecular Wires. J. Phys. Chem. C 2016, 120, 15666-15674.

36 ACS Paragon Plus Environment

Page 37 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

108. Lissau, H.; Frisenda, R.; Olsen, S. T.; Jevric, M.; Parker, C. R.; Kadziola, A.; Hansen, T.; van der Zant, H. S. J.; Brøndsted Nielsen, M.; Mikkelsen, K. V. Tracking Molecular Resonance Forms of Donor–Acceptor Push–Pull Molecules by Single-Molecule Conductance Experiments. Nat. Comm. 2015, 6, 10233. 109. Bouvy, D.; Janousek, Z.; Viehe, H. G.; Tinant, B.; Declercq, J.-P. Cumulogy Supported by X-Ray Analysis of 1,1-Bis(Dimethylamino)-4,4-Dicyanobutatriene. Tetrahedron Lett. 1993, 34, 1779-1782. 110. Tinant, B.; Declercq, J.-P.; Bouvy, D.; Janousek, Z.; Viehe, H. G. Exceptionally Short Central

Bonds

in

the

Cumulogue

and

the

Vinylogue

of

2,2-Dicyano-1,1-

Bis(Dimethylamino)Ethylene. J. Chem. Soc., Perkin Trans. 2 1993, 911-915. 111. Morley, J. O. Theoretical Study of the Electronic Structure and Hyperpolarizabilities of Donor-Acceptor Cumulenes and a Comparison with the Corresponding Polyenes and Polyynes. J. Phys. Chem. 1995, 99, 10166-10174. 112. Valérie, M.; Léo, L.; Aya, H.; Cécile, B.; Alix, S.; Carine, D.; Teruo, S.; Remi, C. 1,4‐ Dialkynylbutatrienes: Synthesis, Stability, and Perspectives in the Chemistry of Carbo‐Benzenes. Chem. Eur. J. 2011, 17, 5086-5100. 113. Gawel, P.; Wu, Y.-L.; Finke, A. D.; Trapp, N.; Zalibera, M.; Boudon, C.; Gisselbrecht, J.P.; Schweizer, W. B.; Gescheidt, G.; Diederich, F. Push–Pull Buta-1,2,3-Trienes: Exceptionally Low Rotational Barriers of Cumulenic C=C Bonds and Proacetylenic Reactivity. Chem. Eur. J. 2015, 21, 6215-6225.

37 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 38 of 47

114. Katagiri, S.; Sudoh, S.; Ozaki, T.; Toda, T.; Shimazaki, N. Dipole Moments of 1Cyclopentadienylidene-2-Cyclo-Heptatrienylidene-Ethylene

and

Its

Related

Compounds.

Tetrahedron Lett. 1992, 33, 6819-6822. 115. Dahlstrand, C.; Rosenberg, M.; Kilså, K.; Ottosson, H. Exploration of the Π-Electronic Structure of Singlet, Triplet, and Quintet States of Fulvenes and Fulvalenes Using the Electron Localization Function. J. Phys. Chem. A 2012, 116, 5008-5017. 116. Toda, T.; Shimazaki, N.; Mukai, T. Synthesis and Properties of 1‐Cycloheptatrienylidene‐ 2‐Cyclopentadienylidene‐Ethylene Derivatives. Angew. Chem. Int. Ed. 1987, 26, 335-336. 117. Möllerstedt, H.; Piqueras, M. C.; Crespo, R.; Ottosson, H. Fulvenes, Fulvalenes, and Azulene:  Are They Aromatic Chameleons? J. Am. Chem. Soc. 2004, 126, 13938-13939. 118. Blake, I. M.; Rees, L. H.; Claridge, T. D. W.; Anderson, H. L. Synthesis and Crystal Structure of a Cumulenic Quinoidal Porphyrin Dimer with Strong Electronic Absorption in the Infrared. Angew. Chem. Int. Ed. 2000, 39, 1818-1821. 119. Ueta, K.; Naoda, K.; Ooi, S.; Tanaka, T.; Osuka, A. Meso‐Cumulenic 2h‐Corroles from Meso‐Ethynyl‐3h‐Corroles. Angew. Chem. Int. Ed. 2017, 56, 7223-7226. 120. Haberland, S.; Finke, A. D.; Kerisit, N.; Katan, C.; Trolez, Y.; Gawel, P.; Leito, I.; Lõkov, M.; Järviste, R.; Kaupmees, K. et al., Enhancement of Push–Pull Properties of Pentafulvene and Pentafulvalene Derivatives by Protonation at Carbon. Eur. J. Org. Chem. 2018, 2018, 739-749. 121. Jimenez, V. G.; Tapia, R.; Medel, M. A.; Mariz, I. F. A.; Ribeiro, T.; Blanco, V.; Cuerva, J. M.; Macoas, E.; Campana, A. G. Aggregation-Induced Emission of [3]Cumulenes Functionalized with Heptagon-Containing Polyphenylenes. Chem. Comm. 2018, 54, 3359-3362.

38 ACS Paragon Plus Environment

Page 39 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

122. Stuyver, T.; Zeng, T.; Tsuji, Y.; Fias, S.; Geerlings, P.; De Proft, F. Captodative Substitution: A Strategy for Enhancing the Conductivity of Molecular Electronic Devices. J. Phys. Chem. C 2018, 122, 3194-3200. 123. Kerisit, N.; Finke, A. D.; Trapp, N.; Leroux, Y. R.; Guillemin, J.-C.; Trolez, Y.; Diederich, F. New Reactivity of 6,6-Bis-Donor-Substituted Pentafulvenes: One-Step Synthesis of Highly Substituted [3]Cumulene and Dihydropentalene. Tetrahedron 2015, 71, 4393-4399. 124. Aviram, A.; Ratner, M. A. Molecular Rectifiers. Chem. Phys. Lett. 1974, 29, 277-283. 125. Metzger, R. M. Unimolecular Electronics. Chem. Rev. 2015, 115, 5056-5115. 126. Capozzi, B.; Xia, J.; Adak, O.; Dell, E. J.; Liu, Z.-F.; Taylor, J. C.; Neaton, J. B.; Campos, L. M.; Venkataraman, L. Single-Molecule Diodes with High Rectification Ratios Through Environmental Control. Nat. Nanotech. 2015, 10, 522. 127. Li, Y.; Zhao, J.; Yin, X.; Yin, G. Ab Initio Investigations of the Electric Field Dependence of the Geometric and Electronic Structures of Molecular Wires. J. Phys. Chem. A 2006, 110, 1113011135. 128. Ye, Y.; Zhang, M.; Zhao, J. Ab Initio Investigations on Three Isomers of Polyacetylene under the Interaction of the External Electric Field. J. Mol. Struct. (Theochem.) 2007, 822, 12-20. 129. Li, Y.; Zhao, J.; Yin, G. Theoretical Investigations of Oligo(Phenylene Ethylene) Molecular Wire: Effects from Substituents and External Electric Field. Comput. Mater. Sci. 2007, 39, 775-781. 130. Movsisyan, L. D.; Kondratuk, D. V.; Franz, M.; Thompson, A. L.; Tykwinski, R. R.; Anderson, H. L. Synthesis of Polyyne Rotaxanes. Org. Lett. 2012, 14, 3424-3426.

39 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 40 of 47

131. Weisbach, N.; Baranova, Z.; Gauthier, S.; Reibenspies, J. H.; Gladysz, J. A. A New Type of Insulated Molecular Wire: A Rotaxane Derived from a Metal-Capped Conjugated Tetrayne. Chem. Comm. 2012, 48, 7562-7564. 132. Movsisyan, L. D.; Franz, M.; Hampel, F.; Thompson, A. L.; Tykwinski, R. R.; Anderson, H. L.Polyyne Rotaxanes: Stabilization by Encapsulation. J. Am. Chem. Soc. 2016, 138, 1366-1376. 133. Milan, D. C.; Krempe, M.; Ismael, A. K.; Movsisyan, L. D.; Franz, M.; Grace, I.; Brooke, R. J.; Schwarzacher, W.; Higgins, S. J.; Anderson, H. L., et al., The Single-Molecule Electrical Conductance of a Rotaxane-Hexayne Supramolecular Assembly. Nanoscale 2017, 9, 355-361. 134. Baisch, B.; Raffa, D.; Jung, U.; Magnussen, O. M.; Nicolas, C.; Lacour, J.; Kubitschke, J.; Herges, R. Mounting Freestanding Molecular Functions Onto Surfaces: The Platform Approach. J. Am. Chem. Soc. 2009, 131, 442-443. 135. Wei, Z.; Wang, X.; Borges, A.; Santella, M.; Li, T.; Sørensen, J. K.; Vanin, M.; Hu, W.; Liu, Y.; Ulstrup, J. et al., Triazatriangulene as Binding Group for Molecular Electronics. Langmuir 2014, 30, 14868-14876. 136. Jasper-Tönnies, T.; Garcia-Lekue, A.; Frederiksen, T.; Ulrich, S.; Herges, R.; Berndt, R., Conductance of a Freestanding Conjugated Molecular Wire. Phys. Rev. Lett. 2017, 119, 066801. 137. Valášek, M.; Mayor, M. Spatial and Lateral Control of Functionality by Rigid Molecular Platforms. Chem. Eur. J. 2017, 23, 13538-13548. 138. Chen, W.; Li, H.; Widawsky, J. R.; Appayee, C.; Venkataraman, L.; Breslow, R. Aromaticity Decreases Single-Molecule Junction Conductance. J. Am. Chem. Soc. 2014, 136, 918920.

40 ACS Paragon Plus Environment

Page 41 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

139. Borges, A.; Solomon, G. C. Effects of Aromaticity and Connectivity on the Conductance of Five-Membered Rings. J. Phys. Chem. C 2017, 121, 8272-8279. 140. Zhang, G.-P.; Xie, Z.; Song, Y.; Wei, M.-Z.; Hu, G.-C.; Wang, C.-K. Is There a Specific Correlation Between Conductance and Molecular Aromaticity in Single-Molecule Junctions? Org. Electron. 2017, 48, 29-34. 141. Fujii, S.; Marqués-González, S.; Shin, J.-Y.; Shinokubo, H.; Masuda, T.; Nishino, T.; Arasu, N. P.; Vázquez, H.; Kiguchi, M. Highly-Conducting Molecular Circuits Based on Antiaromaticity. Nat. Comm. 2017, 8, 15984. 142. Stuyver, T.; Perrin, M.; Geerlings, P.; De Proft, F.; Alonso, M. Conductance Switching in Expanded Porphyrins Through Aromaticity and Topology Changes. J. Am. Chem. Soc. 2018, 140, 1313-1326. 143. Montgomery, L. K.; Huffman, J. C.; Jurczak, E. A.; Grendze, M. P. The Molecular Structures of Thiele's and Chichibabin's Hydrocarbons. J. Am. Chem. Soc. 1986, 108, 6004-6011. 144. Porter, W. W.; Vaid, T. P.; Rheingold, A. L. Synthesis and Characterization of a Highly Reducing Neutral “Extended Viologen” and the Isostructural Hydrocarbon 4,4‘ ‘‘ ‘-Di-N-Octyl-PQuaterphenyl. J. Am. Chem. Soc. 2005, 127, 16559-16566. 145. Zhu, X.; Tsuji, H.; Nakabayashi, K.; Ohkoshi, S.-i.; Nakamura, E. Air- and Heat-Stable Planar Tri-P-Quinodimethane with Distinct Biradical Characteristics. J. Am. Chem. Soc. 2011, 133, 16342-16345. 146. Eaves, S. G.; Hart, S. J.; Whitwood, A. C.; Yufit, D. S.; Low, P. J.; Lynam, J. M. Rapid Markovnikov Addition of Hcl to a Pendant Alkyne: Evidence for a Quinoidal Cumulene. Chem. Comm. 2015, 51, 9362-9365.

41 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 42 of 47

147. Gruber, M.; Padberg, K.; Min, J.; Waterloo, A. R.; Hampel, F.; Maid, H.; Ameri, T.; Brabec, C. J.; Tykwinski, R. R. Acenequinocumulenes: Lateral and Vertical Π‐Extended Analogues of Tetracyanoquinodimethane (Tcnq). Chem. Eur. J. 2017, 23, 17829-17835. 148. Barry, B. M.; Soper, R. G.; Hurmalainen, J.; Mansikkamäki, A.; Robertson, K. N.; McClennan, W. L.; Veinot, A. J.; Roemmele, T. L.; Werner‐Zwanziger, U.; Boeré, R. T. et al., Mono‐ and Bis(Imidazolidinium Ethynyl) Cations and Reduction of the Latter to Give an Extended Bis‐1,4‐([3]Cumulene)‐P‐Carboquinoid System. Angew. Chem. Int. Ed. 2018, 57, 749-754. 149. Hansmann, M. M.; Melaimi, M.; Munz, D.; Bertrand, G. Modular Approach to Kekulé Diradicaloids Derived from Cyclic (Alkyl)(Amino)Carbenes. J. Am. Chem. Soc. 2018, 140, 25462554. 150. Cadierno, V.; Gimeno, J. Allenylidene and Higher Cumulenylidene Complexes. Chem. Rev. 2009, 109, 3512-3560. 151. Coletti, C.; Marrone, A.; Re, N. Metal Complexes Containing Allenylidene and Higher Cumulenylidene Ligands: A Theoretical Perspective. Acc. Chem. Res. 2012, 45, 139-149. 152. Pu, L.; Zhang, Z.; King, R. B.; Allen, W. D. Most Favorable Cumulenic Structures in IronCapped Linear Carbon Chains Are Short Singlet Odd-Carbon Dications: A Theoretical View. Phys. Chem. Chem. Phys. 2018, 20, 15496-15506. 153. Auger, N.; Touchard, D.; Rigaut, S.; Halet, J.-F.; Saillard, J.-Y. Electronic Structure of Ruthenium Cumulene Complexes [Cl(Ph3)4rucnh2]+ (N = 1−8) and of Their Reduced States. Bonding and Properties of the Cationic, Neutral, and Anionic Series with Respect to the Cumulenic Chain Length. Organometallics 2003, 22, 1638-1644.

42 ACS Paragon Plus Environment

Page 43 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

154. Touchard, D.; Haquette, P.; Daridor, A.; Toupet, L.; Dixneuf, P. H. First Isolable Pentatetraenylidene Metal Complex Containing the Ru:C:C:C:C:Cph2 Assembly. A Key Intermediate to Provide Functional Allenylidene Complexes. J. Am. Chem. Soc. 1994, 116, 1115711158. 155. Roth, G.; Fischer, H. Complexes with Diamino-Substituted Unsaturated C3 and C5 Ligands:  First Group 6 Pentatetraenylidenes and New Allenylidene Complexes. Organometallics 1996, 15, 1139-1145. 156. Roth, G.; Fischer, H. On the Way to Heptahexaenylidene Complexes:  Trapping of an Intermediate with the Novel M=C=C=C=C=C=C=Cr2 Moiety. Organometallics 1996, 15, 57665768. 157. Lass, R. W.; Steinert, P.; Wolf, J.; Werner, H. Synthesis and Molecular Structure of the First Neutral Transition‐Metal Complex Containing a Linear M=C=C=C=C=Cr2 Chain. Chem. Eur. J. 1996, 2, 19-23. 158. Kerstin, I.; Helmut, W. The First Structurally Characterized Metal Complex with the Molecular Unit M=C=C=C=Cr2. Angew. Chem. Int. Ed. 2000, 39, 1632-1634. 159. Roth, G.; Fischer, H.; Meyer-Friedrichsen, T.; Heck, J.; Houbrechts, S.; Persoons, A. Synthesis and Nonlinear Optical Properties of New Heptapentaenylidene Complexes:  Study on the Second Harmonic Generation Efficiencies of Amino-Substituted Group 6 Cumulenylidenes. Organometallics 1998, 17, 1511-1516. 160. Zheng, Q.; Gladysz, J. A. A Synthetic Breakthrough Into an Unanticipated Stability Regime:  Readily Isolable Complexes in Which C16−C28 Polyynediyl Chains Span Two Platinum Atoms. J. Am. Chem. Soc. 2005, 127, 10508-10509.

43 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 44 of 47

161. Stahl, J.; Bohling, J. C.; Bauer, E. B.; Peters, T. B.; Mohr, W.; Martín‐Alvarez, J. M.; Hampel, F.; Gladysz, J. A. Sp Carbon Chains Surrounded by Sp3 Carbon Double Helices: A Class of Molecules That Are Accessible by Self‐Assembly and Models for “Insulated” Molecular‐Scale Devices. Angew. Chem. Int. Ed. 2002, 41, 1871-1876. 162. Stahl, J.; Mohr, W.; de Quadras, L.; Peters, T. B.; Bohling, J. C.; Martín-Alvarez, J. M.; Owen, G. R.; Hampel, F.; Gladysz, J. A. Sp Carbon Chains Surrounded by Sp3 Carbon Double Helices:  Coordination-Driven Self-Assembly of Wirelike Pt(C⋮C)Npt Moieties That Are Spanned by Two P(Ch2)Mp Linkages. J. Am. Chem. Soc. 2007, 129, 8282-8295. 163. Schwarz, F.; Kastlunger, G.; Lissel, F.; Riel, H.; Venkatesan, K.; Berke, H.; Stadler, R.; Lörtscher, E. High-Conductive Organometallic Molecular Wires with Delocalized Electron Systems Strongly Coupled to Metal Electrodes. Nano Lett. 2014, 14, 5932-5940. 164. Cao, Z.; Xi, B.; Jodoin, D. S.; Zhang, L.; Cummings, S. P.; Gao, Y.; Tyler, S. F.; Fanwick, P. E.; Crutchley, R. J.; Ren, T. Diruthenium–Polyyn-Diyl–Diruthenium Wires: Electronic Coupling in the Long Distance Regime. J. Am. Chem. Soc. 2014, 136, 12174-12183. 165. Al-Owaedi, O. A.; Bock, S.; Milan, D. C.; Oerthel, M.-C.; Inkpen, M. S.; Yufit, D. S.; Sobolev, A. N.; Long, N. J.; Albrecht, T.; Higgins, S. J. et al., Insulated Molecular Wires: Inhibiting Orthogonal Contacts in Metal Complex Based Molecular Junctions. Nanoscale 2017, 9, 9902-9912. 166. Tanaka, Y.; Kiguchi, M.; Akita, M. Inorganic and Organometallic Molecular Wires for Single‐Molecule Devices. Chem. Eur. J. 2017, 23, 4741-4749. 167. Milan, D. C.; Vezzoli, A.; Planje, I. J.; Low, P. J. Metal Bis(Acetylide) Complex Molecular Wires: Concepts and Design Strategies. Dalton Trans. 2018.

44 ACS Paragon Plus Environment

Page 45 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

TOC Figure

45 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 46 of 47

Author Biographies Marc Hamilton Garner received his MSc in 2016 and is currently a PhD student in Gemma Solomon’s group at the University of Copenhagen. His research interests revolve around theoretical chemistry and electron transport with specific research interests including molecular orbital theory, substituent effects, sigma-conjugation in silanes, and quantum interference effects.

William Bro-Jørgensen is currently a BSc student in Gemma Solomon’s group at the University of Copenhagen. His research interests include molecular orbital theory and its deployment in organometallic systems.

Pernille Dalsgaard Pedersen received her BSc in 2017 and is currently a MSc student at the University of Copenhagen. Her research interests are mainly in the field of computational chemistry including its applications in atmospheric chemistry and chemical kinetics.

Gemma C. Solomon majored in Chemical Physics at the University of Western Australia before moving to the University of Sydney where she completed her BSc(Hons) in 2003, for which she was awarded the University Medal. She received her PhD in Chemistry in 2007 from the University of Sydney where she worked with Prof. Jeffrey Reimers and Prof. Noel Hush. She moved to Northwestern University, USA where she was a postdoctoral fellow with Prof. Mark Ratner until moving to the University of Copenhagen in 2010 to start her independent career as an assistant professor in Chemistry. She is currently an Associate Professor in the Nano-Science Center and Department of Chemistry at the University of Copenhagen where her research focus is

46 ACS Paragon Plus Environment

Page 47 of 47 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

charge and heat transport through molecules. Since 2014 she has been a Senior Editor for the Journal of Physical Chemistry A/B/C.

47 ACS Paragon Plus Environment