Polarizability as a Molecular Descriptor for ... - ACS Publications

Oct 21, 2016 - †School of Molecular Sciences and ‡School of Electrical Computer and Energy Engineering, Arizona State University, Tempe, Arizona 8...
0 downloads 5 Views 5MB Size
Subscriber access provided by Vanderbilt Libraries

Article

Polarizability as a Molecular Descriptor for Conductance in Organic Molecular Circuits Shobeir K.S. Mazinani, Reza Vatan Meidanshahi, Julio L Palma, Pilarisetty Tarakeshwar, Thorsten Hansen, Mark A. Ratner, and Vladimiro Mujica J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b06241 • Publication Date (Web): 21 Oct 2016 Downloaded from http://pubs.acs.org on October 25, 2016

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 free 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 accessible to all readers and 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.

The Journal of Physical Chemistry C 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 24

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

Polarizability As A Molecular Descriptor For Conductance In Organic Molecular Circuits. Shobeir K. S. Mazinani,† Reza Vatan Meidanshahi,‡ Julio L. Palma,¶ Pilarisetty Tarakeshwar,† Thorsten Hansen,§ Mark A. Ratner,k and Vladimiro Mujica∗,† School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1604, USA, School of Electrical Computer and Energy Engineering, Arizona State University, Tempe, Arizona 85287-8806, USA, Department of Chemistry, The Pennsylvania State University, Fayette, The Eberly Campus. 2201 University Drive, Lemont Furnace, Pennsylvania 15456, USA, Department of Chemistry, University of Copenhagen, Universitetsparken 5, 2100, Copenhagen,Ø, Denmark, and Department of Chemistry, Northwestern University, 1145 Sheridan Road, Evanston, Illinois 60208-3113, USA E-mail: [email protected]

To whom correspondence should be addressed School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1604, USA ‡ School of Electrical Computer and Energy Engineering, Arizona State University, Tempe, Arizona 852878806, USA ¶ Department of Chemistry, The Pennsylvania State University, Fayette, The Eberly Campus. 2201 University Drive, Lemont Furnace, Pennsylvania 15456, USA § Department of Chemistry, University of Copenhagen, Universitetsparken 5, 2100, Copenhagen,Ø, Denmark k Department of Chemistry, Northwestern University, 1145 Sheridan Road, Evanston, Illinois 60208-3113, USA ∗ †

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

Abstract We explore a connection between the static molecular polarizability and the molecular conductance that arises naturally in the description of electrified molecular interfaces and that has recently been explored experimentally. We have tested this idea by using measured conductance of few different experimental design motifs for molecular junctions and relating them to the molecular polarizability. Our results show that, for a family of structurally connected molecules, the conductance decreases as the molecular polarizability increases. Within the limitations of our model, this striking result is consistent with the physically intuitive picture that a molecule in a junction behaves as a dielectric that is polarized by the applied bias, hence creating an interfacial barrier that hinders tunneling. The use of the polarizability as a descriptor of molecular conductance offers significant conceptual and practical advantages over a picture based on molecular orbitals. To further illustrate the plausibility of this idea, we have used Simmons' tunneling model that incorporates image charge and dielectric effects to describe transport through a barrier that represents the molecular junction. In such a model, the barrier height depends on the effective dielectric constant of the electrode-molecule-electrode junction, which in turn can be approximately expressed in terms of the molecular polarizability via the classical Clausius-Mossotti relation. Despite the simplicity of our model, it sheds light on a hitherto neglected connection between molecular polarizability and conductance, and paves the way for further experimental, conceptual and theoretical developments.

INTRODUCTION Organic molecular electronics have garnered tremendous interest 1–5 due to the possible ability to operate as verstaile circuit elements. 5–7 This is outshone by the recently explored interface of molecular electronic devices with biochemistry 8 and the study of various phenomena such as charge transfer at the bio-nano interface. 9 One of the key properties for an 2 ACS Paragon Plus Environment

Page 2 of 24

Page 3 of 24

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

efficient design is the ability of these molecular structures to transport electrons effectively. There are endless numbers of possibilities for the molecular structures. 10 While having a plethora of options for design gives scientists and engineers many opportunities, it can become overwhelming in the absence of a physical parameter that can be used as a screening factor that can a priori discern between desirable and unattractive candidates. 7,11 For example, such a factor can be very important in finding an optimized reader molecule in a recognition junction since synthesis and measuring or computing the conductance of these systems is very time consuming. 12,13 Traditionally, some rules of thumb have been used to make fast predictions about the conductance in a single molecule junction, and parameters such as length, 14 temperature and the energy gap between HOMO and LUMO have been considered in the literature. 15 But due to the intrinsic complexity of transport and its relation to molecular electronic structure, these rules do not always hold. 16,17 In this context, molecular polarizability can be used as a guideline for designing systems with desired transport properties.As important as it might be, this type of relatively simple correlation must be taken with caution because conductance depends on some other factors, e.g. the detailed coupling of the molecule to the electrodes and non-linear voltage effects that cannot possibly be captured in a molecule-only description. The physical plausibility of the connection between conductance and polarizability can be inferred from Mujica and coworkers’ work that described the spatial profile of the electrostatic potential in a junction by solving Schrödinger and Poisson equations self-consistently by connecting the quantum electronic density to the electrostatic potential. 18 The picture that emerges from this model is that if many-body and inelastic effects, such as charging and electron-phonon coupling, are neglected molecules behave to a large extent as a dielectric whose polarization response counteracts the driving field (See figure 1). 18

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

Figure 1: The figure displays the spatial voltage profile across a junction. For a vacuum junction (left panel) the voltage V drops linearly between the contacts, as predicted by Poisson equation. For a molecular junction, assuming that the molecule behaves as a dielectric, which polarizes in opposite direction to the applied voltage causing the onset of molecular charges (δ − δ + ) at the interface. This is accompanied by a voltage drop essentially at the interface.

This leads to a spatial profile that differs substantially from that of a vacuum junction between two electrodes, which for a one dimensional system, is a linear function of the inter-electrode separation as found by solving Poisson equation for zero charge density, and corresponds rather to a S-shape function associated with a spatial profile characterized by the fact that the potential drop occurs essentially at the interfaces between the molecule and the electrodes. 18 Once the crucial role of the molecular bridge in determining the local dielectric properties of a junction has been accepted, the next conceptual step is to connect the dielectric constant to the molecular polarizability, which we introduce by assuming the validity of the ClausiusMosotti relation from electromagnetism, a very important aspect of our model that will be discussed below. Once this is achieved, the missing link is the connection between the conductance itself and the molecular polarizability. A strong theoretical support for this connection for quasi-metallic junctions has been addressed previously. 19,20 . In this work, we approach this problem in its simplest form, assuming that a molecular junction can be described as a simple tunneling junction and use Simmons’ model, which connects the barrier properties, that in turn determine the current and the conductance, to the junction’s dielectric constant. 4 ACS Paragon Plus Environment

Page 4 of 24

Page 5 of 24

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 main goal of this study is to explore the robustness of this connection for different design motifs and discuss the domain in which it is applicable. To this end, we have investigated the correlation between calculated molecular polarizabilities and the experimental zero voltage conductance of different families of molecules used in junctions, and our results are remarkably consistent with the qualitative and even quantitative predictions of our simple model. We finish by discussing the limitations and possible applications of the use of polarizability as a molecular descriptor for conductance in helping the design of next generation nano-electronics components.

METHODS Geometry optimizations and polarizability calculations were carried out at the Density Functional Theory level by using the Becke gradient-corrected exchange functional and Lee-YangParr correlation functional with three parameters (B3LYP) and the 6-31G* basis set by ORCA electronic structure package. 21–26

RESULTS AND DISCUSSION As mentioned in the Introduction, we will approach the rational design of molecular electronics materials by studying the correlation between polarizability and conductance. It is worth mentioning that since polarizability is a tensor property and conductance is a scalar quantity, we use as a molecular descriptor the isotropic molecular polarizability which is defined as: 1 1 α ¯ = tr(α) = (αxx + αyy + αzz ) 3 3 We have investigated systems whose experimental conductance has been reported and tested the correlation between calculated polarizabilities and such measured conductance. A first important result is that an increase in the static polarizability of the molecules results in a

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

decrease in the conductance values, as should be expected from the simple model of molecules as dielectrics. Although the functional mathematical representation of this correlation is fairly general, each family of molecules has its own fitting parameters. We classify our results based on different molecular groups and show that the correlation holds within each group, thus stressing the potential of polarizability as a molecular descriptor that could provide valuable guidance for the design of new molecular electronic devices. It should be noted that the isotropic polarizability is clearly the simplest metric to use and other important properties of the polarizability tensor, particularly anisotrpic effects in the direction of orientation of the molecules can be of interest.

Halobenzene, The Effect of Substituents As an initial step, we investigated the effect of halogen substitution on benzene’s conductance. Due to their electronegative nature, halogens are electron withdrawing groups. Hence they reduce the electron density on the benzene’s π-system which results in a decrease in the energy of the HOMO orbital in comparison to HOMO energy of benzene. 27 The conductance values reported by Venkataraman and coworkers are plotted against calculated molecular polarizability; Figure 2 indicates that the conductance of the molecule decreases, although not significantly (< 5%), as its static polarizability increases due to the effect of the substituents. Further systems and design motifs are studied in order to find more conclusive results regarding the connection between polarizability and conductance in molecular junctions.

Amine-Gold Linked Molecular Motif The thiol group has been used extensively as an anchoring group between the molecular bridge and gold electrodes. 4,14,28,29 The strong bond between gold and sulfur is the reason for this general interest towards the Au-S anchoring groups. The main issue with sulfur-gold motif is the dramatic dependence of the conductance on the anchoring groups’ geometrical 6 ACS Paragon Plus Environment

Page 6 of 24

Page 7 of 24

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: The correlation between conductance (milig0 ) and the static molecular polariz3 ability (Å ). In these systems the halogens affect the conductance, not very significantly, through modifying the HOMO of the molecule. It can be inferred that polarizability can be used to predict the effect of this substituents on conductance of benzene. 27 parameters. 27,30,31 Nitrogen-based anchoring groups have been used as an alternative to S-Au bonding interfaces. 32 Due to the reproducibility and smaller effect of the geometry of the linker on conductance; using this amine-based motif in designing the molecular circuits as well as thermoelectric systems seems to be a sensible path. 27,33,34 To investigate the robustness of our conjecture, we calculated the polarizabilities of the diamine family and the bipyridine family and plotted them against the measured conductance values. 27,33,34 Figure 3 summarizes the results.

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

(a)

(b)

(c)

Figure 3: The structures of corresponding compounds are presented. Carbons are represented by black, hydrogens with grey and nitrogens with blue spheres (a) Structures 1, 2 and 3 belong to the diamine family whereas 4, 5 and 6 are from the bipyridine family of molecules (b). (c) Measured conductance (milig0 ), obtained from Venkataraman and 3 coworkers work, are plotted with respect to calculated polarizability (Å ). The figure suggests that for a family of molecules we can use polarizability as a descriptor of conductance to predict the electron transport properties of a molecule.

We attribute the difference between the conductances of the two families to the different nature of their individual anchoring group.

8 ACS Paragon Plus Environment

Page 8 of 24

Page 9 of 24

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

Thiophene and Furan Oligomers Five-membered heterocycles, in particular thiophene and furan, have been used extensively in organic electronic devices such as organic field effect transistors (OFET) and organic photovoltaics (OPV). 35–40 The oligomers can be considered as a cis π-conjugated (CH)x system that are stabilized by a heteroatom. 41 These systems are aromatic. Breslow and coworkers showed that an increase in the (negative) aromatic stabilization energy correlates with the conductance values in systems with thiophene and furan motifs. 42 We have selected the aromatic systems in Breslow’s work, to put the polarizability-based design guideline through a strict test. Compound 7, 8, 9 and 10 (named in accord with Breslow’s work) are 2,5-bis-(4-aminophenylethinyl)furan,. 2,5-bis-(4-aminophenylethinyl)thiophene, 2,5-bis(4-aminophenyl) furan and 2,5-bis-(4-aminophenylethinyl) thiophene respectively (Figure 4-a) . The calculated polarizability and measured conductance are presented in Figure 4-b. There are two sets of correlations for these systems. First, the correlation between 7 and 9 (8 and 10) that are furan (thiophene) based compounds that are tweaked by the addition of an alkyne group. Second, the furan versus thiophene (7, 8 and 9, 10). It can be seen that for both sets of molecular wires, an increase in polarizability is followed by a decrease in conductance of the heterocyclic molecules. It is worth mentioning that polarizability has been used in other studies as an indicator of aromaticity, which establishes a connection between electron delocalization and conductance for this family of compounds. 43,44 Increasing the length of a molecular wire by increasing furan/ thiophene units mainly modifies the LUMO and it has been shown to have a negligible effect on the energy of the HOMO, which again points out to the limitations of schemes based on molecular orbitals to predict the behavior of the conductance . 42,43,45

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

(a)

(b)

(c)

Figure 4: (a) and (b) show the structure of the compounds 7, 8 and 9, 10. In the structures Nitrogen, oxygen, sulfur, carbon and hydrogen are blue, red, yellow,black and grey. (C) The calculated polarizabilities versus experimental conductance shows that for the families of molecules, we observe how the change in the polarizability correlates with the changes in the conductance of the junction.

Hydrogen Bond Motif A ubiquitous motif in organic molecular electronics, in particular bio-inspired designs, whose use permits taking advantage of the unique properties of hydrogen bonding to build low cost, bottom-up organic materials with desirable properties. 46–52 Nishino et al have measured the conductance through hydrogen bond and showed that at short distances hydrogen-bonded wires have higher conductance than alkane chains. 53 Based on their results, we investigate the connection between polarizability and conductance for the family of carboxylic acid alkane thiols. Figure 5 shows the robustness of the correlation between conductance and

10 ACS Paragon Plus Environment

Page 10 of 24

Page 11 of 24

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

static molecular polarizability. We have previously shown for a set of biologically-relevant hydrogen-bonded systems, that polarizability can be used as a guideline to qualitatively predict conductance through the hydrogen bonded system. 54

Figure 5: The correlation between measured conductance and calculated static polarizability of three hydrogen bonded carboxyl groups is shown in this plot. The labels show the total number of carbons in the molecular junction in a manner consistent with the work of Nishino et al . 53

Barrier Model of Conductance To explore the connection between polarizability and conductance, we consider a model of molecular conductance as a tunneling process. This model essentially ignores all the manybody, vibronic interactions, and inelastic aspects of the transport process and assumes that the molecule acts as a one-dimensional tunneling barrier where the barrier is specified by two parameters: the height and the width of the barrier. Simmons’ model, which includes image charges and dielectric effects, has been extensively used for the description of tunneling through metal-molecule interfaces with remarkable success. 55 We are particularly interested ¯ which has been shown in the in the expression for the height of the tunneling barrier, φ, Simmons’ model to be related to the dielectric constant as: 55 s1 + s2 1.15λx s2 (x − s1 ) 1 φ¯ = φ0 − qV − ln( ). 2x ∆s s1 (x − s2 ) r

11 ACS Paragon Plus Environment

(1)

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 24

with λ defined by: λ=

q 2 ln2 8π0 x

(2)

Where q is the electron charge, V is the bias voltage, φ0 is the barrier height, s1 and s2 are the turning points in the barrier shape, ∆s = s2 − s1 , x is the nominal width of the barrier, 0 is the vacuum permittivity and r is the dielectric constant. Since we are not using an atomistic representation of the contacts, the specific values of the turning points have to be taken as parameters of the model. The relation between the current and the voltage in the tunneling junction can be recast in the following form: J = J0 (φ¯ e−A J0 =



φ¯

− (φ¯ + qV ) e−A



¯ φ+qV

q 2πh(σ∆s)2

) (3)

√ A = ( 4π∆s 2me ) h where me is mass of electron, h is the Plank’s constant, and σ is a correction factor as described by Simmons and is usually around 1. 55

The simplest connection between the dielectric constant r of an electrified interface and the molecular polarizability α of the intervening medium is given by Clausius-Mossotti equation: r = γ=

0 +2γα 0 −γα

(4)

NA d 3M

with NA being Avogadro’s number, M is the molar mass of the material and d is its density. It applies to the dielectric constant of a bulk dielectric material that is homogeneous and isotropic, and it connects the static polarizability of a single molecule with the susceptibility of a three-dimensional molecular material. 56,57 The basic microscopic premise of this relation is that in a uniform electric field, each molecule is represented as a polarizable point dipole

12 ACS Paragon Plus Environment

Page 13 of 24

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

that experiences the external field inducing a polarization response. Although these conditions might seem too simplistic to describe a molecular junction, it has been shown that the Clausius- Mossotti equation can accurately represent the relation between polarizability and the dielectric function even for covalently-bonded bulk semiconductors. 58 More interestingly for our current purposes, the nanoscopic validity of the Clausius- Mossotti equation has been systematically explored by Natan et al. 57 and has been shown to correspond to a well-defined limit that provides us with a physically reasonable starting point, that is only applicable to simple junctions where tunneling is the dominant transport mechanism. The inclusion of hopping or hybrid transport mechanisms, like those present in DNA, is well beyond the validity of our simple model. The combined use of Equations 4 and 1 results in the desired connection between the effective barrier’s height and the polarizability: 0 − γα s1 + s2 − B( ) φ¯ = φ0 − qV 2x 0 + 2γα where B =

1.15λx 1) ln( ss21 (x−s ). ∆s (x−s2 )

(5)

The differential conductance g is defined as:

g(V ) =

∂J ∂V

(6)

And the differential conductance in the limit of zero voltage can be obtained straightforwardly from Equation (3). √ lim g(V ) = −qJ0 e−A

V →0



q 

A φ¯ φ¯   1− 2

(7)

In deriving Equation 7, we assumed that the molecular polarizability is not dependent on the bias voltage. Expanding the function representing the barrier height in terms of polarizability and inserting it into equation 7 gives:

g = g1 e

−β1 α

β1 C − α + ... 2

!

13 ACS Paragon Plus Environment

(8)

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

where β1 =

ABγ √ 20

φ0 −B

, g1 = −J0 qe−A



φ0 −B

and C = 1 −

A 2



Page 14 of 24

φ0 − B. It can be seen from

equation 8, polarizability adds more corrections to a normal exponential decay behavior. Also we have shown that other underlying effects such as the role of substituents or the aromaticity that cannot be justified by length dependence can be predicted by polarizability as a molecular descriptor. Inspired by the correlation between polarizability and the classical analogue of the molecular "volume", we write α ' L3 , where L is a typical molecular linear dimension, obtained adding the values of bond lengths. Figure 6 plots Ln of the experimental conductance versus α1/3 and compares it with the conductance values calculated via equation 8 for the diamine and bipyridine systems 3. As it can be seen in the Figure 6, this correlation holds for families of structurally related compounds. For instance, comparison of the polarizability of a substitued-cyclohexane with benzene derivatives will not follow the aforementioned discussion since the modeling parameters for benzene is expected to be different for these two systems.

(a)

(b)

Figure 6: Plot of Ln of conductance versus third root of polarizability for diamine (a) and bipyridine (b). We used the following values and fitted our equation 8. We used the following values to fitting. (a) For diamine we used g1 = 11.26, C = 0.01 and β1 = 3.76. (b) For bipyridine family, we used g1 = 2.75, C = 0.01 and β1 = 4.31 A theoretical formulation beyond a barrier model can be cast in terms of Green’s functions, as is often done for the Landauer theory for molecular conductance. Response properties, such as the molecular polarizability, have been connected to single particle Green’s functions in the literature 59 , also they can be derived invoking the Keldysh contour techniques of non-equilibrium Green’s functions. 60 14 ACS Paragon Plus Environment

Page 15 of 24

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 field-induced dipole moment as a function of time is given by

hµ(t)i = i

X

µij µkl

Z

n

o

< < adv dτ Gret jk (t − τ )Gli (τ − t) + Gjk (t − τ )Gli (τ − t) E(τ ),

(9)

ijkl

where Gret (t), Gadv (t), and G< (t), are the retarded, advanced, and lesser molecular Green’s functions respectively. This implies the familiar expression for molecular polarizability:

α(ω) = −

X

n

o

adv µij µkl Gret jk (ω)%li (−ω) + %li (ω)Gjk (−ω)

(10)

ijkl

Here % = −iG< is the single-particle reduced density matrix. Using this equation, together with the Landauer expression for the conductance, n

g = g0 T r ΓL Gret ΓR Gadv

o

(11)

where g0 is the quantum of conductance and Γ is the spectral density describing the coupling between the molecule and the L,R electrodes, it is possible to derive an expression equivalent to 8. A detailed exposition of this derivation will be reported elsewhere.

CONCLUSIONS In this work we have explored the rather substantial evidence of an existing correlation between the static isotropic molecular polarizability and the molecular contribution to the zero-voltage conductance of a molecular junction for certain families of molecules. We have also examined the physical origin of such a correlation via a model that connects the local dielectric properties of the junction to its transport behavior via changes in the effective shape of the associated tunneling barrier. While our model is clearly not a first-principle theory, it does incorporate strong physical plausibility arguments and provides a rationale for the apparent correlation between

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

experimental conductances and calculated molecular polarizabilites. The use of a molecular property such as polarizability, has clear advantages over correlations based on an orbital picture of conductance and information about the molecular energy spectrum, e.g. the HOMO-LUMO gap, because these models are based on a strong assumption about the nature of conductance channels that breaks down in many important cases. It also relates to recent efforts by Jackson et al to use the Kirchhoff constant as a single variable characterizing the extent of electric network connectivity at the molecular level. 61 Finally, our findings strongly suggest, together with the already demonstrated connection between conductance and electron transfer rate 62–64 , that it should be possible to reformulate Marcus theory of electron transfer in terms of response functions associated to the frequency-dependent polarizabilities, a subject we are currently working on and that is related to the original approach taken by Marcus to this subject.

Acknowledgement V.M. thanks Prof. Abraham Nitzan for interesting discussions about the subject. The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. The authors declare no competing financial interest.

16 ACS Paragon Plus Environment

Page 16 of 24

Page 17 of 24

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

References (1) Lörtscher, E. Wiring Molecules into Circuits. Nat. Nanotechnol. 2013, 8, 381–384. (2) Heath, J. R. Molecular Electronics. Annu. Rev. Mater. Res. 2009, 39, 1–23. (3) Huang, C.; Rudnev, A. V.; Hong, W.; Wandlowski, T. Break Junction under Electrochemical Gating: Testbed for Single-Molecule Electronics. Chem. Soc. Rev. 2015, 44, 889–901. (4) Ratner, M. A Brief History of Molecular Electronics. Nat. Nanotechnol. 2013, 8, 378– 381. (5) Coskun, A.; Spruell, J. M.; Barin, G.; Dichtel, W. R.; Flood, A. H.; Botros, Y. Y.; Stoddart, J. F. High Hopes: Can Molecular Electronics Realise its Potential? Chem. Soc. Rev. 2012, 41, 4827–4859. (6) Asadi, K.; Katsouras, I.; Harkema, J.; Gholamrezaie, F.; Smits, E. C.; Biscarini, F.; Blom, P. W.; de Leeuw, D. M. Organic Field-effect Transistors as a Test-bed for Molecu17 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

lar Electronics: A Combined Study with Large-area Molecular Junctions. Org. Electron. 2012, 13, 2502–2507. (7) Savoie, B. M.; Kohlstedt, K. L.; Jackson, N. E.; Chen, L. X.; Olvera de la Cruz, M.; Schatz, G. C.; Marks, T. J.; Ratner, M. a. Mesoscale Molecular Network Formation in Amorphous Organic Materials. Proc. Natl. Acad. Sci. U. S. A. 2014, 111, 10055–10060. (8) Lindsay, S. Biochemistry and Semiconductor Electronics- The Next Big Hit for Silicon? J. Phys. Condens. Matter 2012, 24, 164201–164208. (9) Tarakeshwar, P.; Palma, J. L.; Holland, G. P.; Fromme, P.; Yarger, J. L.; Mujica, V.; Yarger, L. Probing the Nature of Charge Transfer at Nano-Bio Interfaces: Peptides on Metal Oxide Nanoparticles. J. Phys. Chem. Lett. 2014, 5, 3555–3559. (10) Pathem, B. K.; Claridge, S. a.; Zheng, Y. B.; Weiss, P. S. Molecular Switches and Motors on Surfaces. Annu. Rev. Phys. Chem. 2013, 64, 605–630. (11) Jackson, N. E.; Savoie, B. M.; Chen, L. X.; Ratner, M. a. A Simple Index for Characterizing Charge Transport in Molecular Materials. J. Phys. Chem. Lett. 2015, 1021, 1018–1021. (12) Krsti, P.; Ashcroft, B.; Lindsay, S. Physical Model for Recognition Tunneling. Nanotechnology 2015, 26, 084001–084010. (13) Zhao, Y.; Ashcroft, B.; Zhang, P.; Liu, H.; Sen, S.; Song, W.; Im, J.; Gyarfas, B.; Manna, S.; Biswas, S. et al. Single-Molecule Spectroscopy of Amino Acids and Peptides by Recognition Tunnelling. Nat. Nanotechnol. 2014, 9, 466–473. (14) Lafferentz, L.; Ample, F.; Yu, H.; Hecht, S.; Joachim, C.; Grill, L. Conductance of a Single Conjugated Polymer as a Continuous Function of its Length. Science 2009, 323, 1193–1197.

18 ACS Paragon Plus Environment

Page 18 of 24

Page 19 of 24

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

(15) Nitzan, A. Chemical Dynamics in Condensed Phases; Oxford University Press: New York, 2006; p 741. (16) Solomon, G. C.; Andrews, D. Q.; Van Duyne, R. P.; Ratner, M. a. When Things Are Not as They Seem: Quantum Interference Turns Molecular Electron Transfer "Rules" Upside Down. J. Am. Chem. Soc. 2008, 130, 7788–7789. (17) Xu, B. Q.; Li, X. L.; Xiao, X. Y.; Sakaguchi, H.; Tao, N. J. Electromechanical and Conductance Switching Properties of Single Oligothiophene Molecules. Nano Lett. 2005, 5, 1491–1495. (18) Mujica, V.; Roitberg, A. E.; Ratner, M. Molecular Wire Conductance: Electrostatic Potential Spatial Profile. J. Chem. Phys. 2000, 112, 6834–6839. (19) Verstraete, M. J.; Bokes, P.; Godby, R. W. First-Principles Conductance of Nanoscale Junctions from the Polarizability of Finite Systems. J. Chem. Phys. 2009, 130, 124715– 1–124715–8. (20) Bokes, P.; Godby, R. W. Conductance and Polarization in Quantum Junctions. Phys. Rev. B 2004, 69, 245420–1–245420–8. (21) Neese, F. The ORCA Program System. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2012, 2, 73–78. (22) Devlin, F. J.; Finley, J. W.; Stephens, P. J.; Frisch, M. J. Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields: A Comparison of Local, Nonlocal, and Hybrid Density Functionals. J. Phys. Chem. 1995, 99, 16883–16902. (23) Rassolov, V. A.; Pople, J. A.; Ratner, M. A.; Windus, T. L. 6-31G* Basis Set for Atoms K through Zn. J. Chem. Phys. 1998, 109, 1223.

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

(24) Becke, A. D. Density-functional Thermochemistry.III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648. (25) Lee, C.; Yang, W.; Parr, R. G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B 1988, 37, 785–789. (26) Francl, M. M.; Pietro, W. J.; Hehre, W. J.; Binkley, J. S.; Gordon, M. S.; DeFrees, D. J.; Pople, J. a. Self-consistent Molecular Orbital Methods. XXIII. A Polarization-type Basis Set for Second-row Elements. J. Chem. Phys. 1982, 77, 3654–3665. (27) Quek, S. Y.; Venkataraman, L.; Choi, H. J.; Louie, S. G.; Hybertsen, M. S.; Neaton, J. B. Amine - Gold linked Single-Molecule Circuits: Experiment and Theory. Nano Lett. 2007, 7, 3477–3482. (28) Tao, N. J. Electron Transport in Molecular Junctions. Nat. Nanotechnol. 2006, 1, 173– 181. (29) Aradhya, S. V.; Venkataraman, L. Single-Molecule Junctions Beyond Electronic Transport. Nat. Nanotechnol. 2013, 8, 399–410. (30) Dell, E. J.; Capozzi, B.; Dubay, K. H.; Berkelbach, T. C.; Moreno, J. R.; Reichman, D. R.; Venkataraman, L.; Campos, L. M. Impact of Molecular Symmetry on Single-Molecule Conductance. J. Am. Chem. Soc. 2013, 135, 11724–11727. (31) Quek, S. Y.; Khoo, K. H. Predictive DFT-Based Approaches to Charge and Spin Transport in Single-Molecule Junctions and Two-Dimensional Materials: Successes and Challenges. Acc. Chem. Res. 2014, 47, 3250–3257. (32) Zhou, X. S.; Chen, Z. B.; Liu, S. H.; Jin, S.; Liu, L.; Zhang, H. M.; Xie, Z. X.; Jiang, Y. B.; Mao, B. W. Single Molecule Conductance of Dipyridines with Conjugated Ethene and Nonconjugated Ethane Bridging Group. J. Phys. Chem. C 2008, 112, 3935–3940. 20 ACS Paragon Plus Environment

Page 20 of 24

Page 21 of 24

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

(33) Venkataraman, L.; Klare, J. E.; Nuckolls, C.; Hybertsen, M. S.; Steigerwald, M. L. Dependence of Single-Molecule Junction Conductance on Molecular Conformation. Nature 2006, 442, 904–907. (34) Quek, S. Y.; Choi, H. J.; Louie, S. G.; Neaton, J. B. Thermopower of Amine - GoldLinked, Aromatic Molecular Junctions from First Principles. ACS Nano 2011, 5, 551– 557. (35) Li, Y.; Sonar, P.; Singh, S. P.; Zeng, W.; Soh, M. S. 3,6-Di(furan-2-yl)pyrrolo[3,4c]pyrrole-1,4(2H,5H)-dione and Bithiophene Copolymer with Rather Disordered Chain Orientation Showing High Mobility in Organic Thin Film Transistors. J. Mater. Chem. 2011, 21, 10829–10835. (36) Roncali, J.; Leriche, P.; Blanchard, P. Molecular Materials for Organic Photovoltaics: Small is Beautiful. Adv. Mater. 2014, 26, 3821–3838. (37) Scharber, M. C.; Mühlbacher, D.; Koppe, M.; Denk, P.; Waldauf, C.; Heeger, A. J.; Brabec, C. J. Design Rules for Donors in Bulk-heterojunction Solar Cells - Towards 10 % Energy-conversion Efficiency. Adv. Mater. 2006, 18, 789–794. (38) Woo, C. H.; Beaujuge, P. M.; Holcombe, T. W.; Lee, O. P.; Fréchet, J. M. J. Incorporation of Furan into Low Band-gap Polymers for Efficient Solar Cells. J. Am. Chem. Soc. 2010, 132, 15547–15549. (39) Wang, C.; Dong, H.; Hu, W.; Liu, Y.; Zhu, D. Semiconducting Pi-conjugated Systems in Field-Effect Transistors: A Material Odyssey of Organic Electronics. Chem. Rev. 2012, 112, 2208–2267. (40) Niimi, K.; Mori, H.; Miyazaki, E.; Osaka, I.; Kakizoe, H.; Takimiya, K.; Adachi, C. [2,2âĂš]Bi[naphtho[2,3-b]furanyl]: A Versatile Organic Semiconductor with a Furanfuran Junction. Chem. Commun. 2012, 48, 5892–5894.

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

(41) Roncali, J. Conjugated Poly(thiophenes) - Synthesis, Functionalization, and Applications. Chem. Rev. 1992, 92, 711–738. (42) 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, 918–920. (43) Delaere, D.; Nguyen, M. T.; Vanquickenborne, L. G. Influence of Building Block Aromaticity in the Determination of Electronic Properties of Five-membered Heterocyclic Oligomers. Phys. Chem. Chem. Phys. 2002, 4, 1522–1530. (44) De Proft, F.; Geerlings, P. Conceptual and Computational DFT in the Study of Aromaticity. Chem. Rev. 2001, 101, 1451–1464. (45) Bell, S.; Dines, T. J. Advanced Chemistry Classroom and Laboratory Using Modern Electronic Structure Methods. J. Chem. Educ. 2007, 84, 1364–1370. (46) Dujardin, E.; Mann, S. Bio-Inspired Materials Chemistry. Adv. Mater. 2002, 14, 775– 788. (47) Etter, M. Hydrogen Bonds as Design Elements in Organic Chemistry. J. Phys. Chem. 1991, 95, 4601–4610. (48) Babu, S. S.; Prasanthkumar, S.; Ajayaghosh, A. Self-assembled Gelators for Organic Electronics. Angew. Chemie - Int. Ed. 2012, 51, 1766–1776. (49) Forrest, S. R. The Path to Ubiquitous and Low-cost Organic Electronic Appliances on Plastic. Nature 2004, 428, 911–918. (50) Horiuchi, S.; Tokura, Y. Organic Ferroelectrics. Nat. Mater. 2008, 7, 357–366. (51) Zhang, M. X.; Zhao, G. J. Modification of N-type Organic Semiconductor Performance of Perylene Diimides by Substitution in Different Positions: Two-dimensional Pi-stacking and Hydrogen bonding. ChemSusChem 2012, 5, 879–887. 22 ACS Paragon Plus Environment

Page 22 of 24

Page 23 of 24

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

(52) Xia, F.; Jiang, L. Bio-inspired, Smart, Multiscale Interfacial Materials. Adv. Mater. 2008, 20, 2842–2858. (53) Nishino, T.; Hayashi, N.; Bui, P. T. Direct Measurement of Electron Transfer through a Hydrogen Bond between Single Molecules. J Am Chem Soc 2013, 135, 4592–4595. (54) Meidanshahi, R. V.; Mazinani, S. K. S.; Mujica, V.; Tarakeshwar, P. Electronic Transport across Hydrogen Bonds in Organic Electronics. Int. J. Nanotechnol. 2015, 12, 297–312. (55) Simmons, J. G. Generalized Formula for the Electric Tunnel Effect between Similar Electrodes Separated by a Thin Insulating Film. J. Appl. Phys. 1963, 34, 2581–2590. (56) Feynman, R. P.; Leighton, R. B.; Sands, M. The Feynman Lectures on Physics; Basic books, 2013; Vol. 1. (57) Natan, A.; Kuritz, N.; Kronik, L. Polarizability, Susceptibility, and Dielectric Constant of Nanometer-Scale Molecular Films: A Microscopic View. Adv. Funct. Mater. 2010, 20, 2077–2084. (58) Zhang, G.-X.; Tkatchenko, A.; Paier, J.; Appel, H.; Scheffler, M. Van der Waals Interactions in Ionic and Semiconductor Solids. Phys. Rev. Lett. 2011, 107, 245501. (59) Deleuze, M. S.; Pickup, B. T.; Wilton, D. J. Calculation of Molecular Response Properties with the Second-order Coupled Perturbed Electron Propagator. International Journal of Quantum Chemistry 2000, 77, 625–640. (60) Hansen, T.; Pullerits, T. Nonlinear Response Theory on the Keldysh Contour. Journal of Physics B: Atomic, Molecular and Optical Physics 2012, 45, 154014. (61) Jackson, N. E.; Savoie, B. M.; Chen, L. X.; Ratner, M. A. A Simple Index for Characterizing Charge Transport in Molecular Materials. J. Phys. Chem. Lett. 2015, 6, 1018–1021. 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

(62) Nitzan, A. A Relationship between Electron-Transfer Rates and Molecular Conduction. J. Phys. Chem. A 2001, 105, 2677–2679. (63) Mujica, V.; Kemp, M.; Ratner, M. Electron Conduction in Molecular Wires. II. Application to Scanning Tunneling Microscopy. J.Chem. Phys. 1994, 101, 6856–6864. (64) Wierzbinski, E.; Venkatramani, R.; Davis, K. L.; Bezer, S.; Kong, J.; Xing, Y.; Borguet, E.; Achim, C.; Beratan, D. N.; Waldeck, D. H. The Single-Molecule Conductance and Electrochemical Electron-Transfer Rate are Related by a Power Law. ACS Nano 2013, 7, 5391–5401.

24 ACS Paragon Plus Environment

Page 24 of 24