Contact atomic structure and electron transport through molecules

Using benzene sandwiched between two Au leads as a model system, we investigate from first principles the change in molecular conductance caused by different atomic structures around the metal-molecule contact. Our motivation is the variable situations that may arise in break junction experiments; our approach is a combined density functional theory and Green function technique. We focus on effects caused by (1) the presence of an additional Au atom at the contact and (2) possible changes in the molecule-lead separation. The effects of contact atomic relaxation and two different lead orientations are fully considered. We find that the presence of an additional Au atom at each of the two contacts will increase the equilibrium conductance by up to two orders of magnitude regardless of either the lead orientation or different group-VI anchoring atoms. This is due to a resonance peak near the Fermi energy from the lowest energy unoccupied molecular orbital. In the nonequilibrium properties, the resonance peak manifests itself in a negative differential conductance. We find that the dependence of the equilibrium conductance on the molecule-lead separation can be quite subtle: either very weak or very strong depending on the separation regime.     

Theoretical investigation on electron transport through an organic molecule: effect of the contact structure

Knowing how the contact geometry influences the conductance of a molecular wire junction requires both a precise determination of the molecule/metallic-electrode interface structure and an evaluation of the conductance for different contact geometries with a fair accuracy. With a greatly improved method to solve the Lippmann-Schwinger equation, we are able to include at least one atomic layer of each electrode into the extended molecule. The artificial effect of the jellium model used for the electrodes is therefore significantly reduced. Our first-principles calculations on the transport properties of a single benzene dithiolate molecule sandwiched between Au(111) surfaces show that the transmission of the bridge site contact, which is the most stable adsorption configuration in equilibrium, displays different features from those of other configurations, and that the inclusion of the surface layers of Au electrodes into the extended molecule shifts and broadens the transmission peaks due to a stronger and more realistic S-Au bonding. We discuss the geometry dependence of the transport properties by analyzing the density of states of the molecular orbitals.     

Models of electrodes and contacts in molecular electronics

Bridging the difference in atomic structure between experiments and theoretical calculations and exploring quantum confinement effects in thin electrodes (leads) are both important issues in molecular electronics. To address these issues, we report here, by using Au-benzenedithiol-Au as a model system, systematic investigations of different models for the leads and the lead-molecule contacts: leads with different cross sections, leads consisting of infinite surfaces, and surface leads with a local nanowire or atomic chain of different lengths. The method adopted is a nonequilibrium Green's-function approach combined with density-functional theory calculations for the electronic structure and transport, in which the leads and molecule are treated on the same footing. It is shown that leads with a small cross section will lead to large oscillations in the transmission function T(E), which depend significantly on the lead structure (orientation) because of quantum waveguide effects. This oscillation slowly decays as the lead width increases, with the average approaching the limit given by infinite surface leads. Local nanowire structures around the contacts induce moderate fluctuations in T(E), while a Au atomic chain (including a single Au apex atom) at each contact leads to a significant conductance resonance.     

Bond fluctuation of S/Se anchoring observed in single-molecule conductance measurements using the point contact method with scanning tunneling microscopy

Conductance was measured for the single molecules with S/Se anchoring on a Au surface using the point contact method with scanning tunneling microscopy that enables us to selectively perform a repeated analysis of a chosen target molecule. Apparent conductance changes observed in sequential measurements suggest the existence of bond fluctuation among the adsorption sites.     

Spin-related electronic pathway through single molecule on Au(111)

Spin properties of organic molecules have attracted great interest for their potential applications in spintronic devices and quantum computing. Fe-tetraphenyl porphyrin (FeTPP) is of particular interest for its robust magnetic properties on metallic substrates. FeTPP is prepared in vacuum via on-surface synthesis. Molecular structure and spin-related transport properties are characterized by low-temperature scanning tunneling microscope and spectroscopy at 0.5 K. Density functional theory calculations are performed to understand molecular adsorption and spin distribution on Au(111). The molecular structure of FeTPP is distorted upon adsorption on the substrate. Spin excitations of FeTPP are observed on the Fe atom and high pyrrole groups in differential conductance spectra. The calculated spin density distribution indicates that the electron spin of FeTPP is mainly distributed on the Fe atom. The atomic transmission calculation indicates that electrons transport to substrate is mediated through Fe atom, when the tip is above the high pyrrole group.     

噻吩在M(111)(M=Pd,Pt, Au)表面的吸附

采用密度泛函理论(DFT), 选取DMol³程序模块, 对噻吩在M(111) (M=Pd, Pt, Au)表面上的吸附行为进行了探讨. 通过对噻吩在不同底物金属上的吸附能、吸附构型、Mulliken 电荷布居、差分电荷密度以及态密度的分析发现, 噻吩在Pd(111)面上的吸附能最大, Pt(111)面次之, Au(111)面最小. 吸附后, 噻吩在Au(111)面上的构型几乎保持不变, 最终通过S端倾斜吸附于top 位; 噻吩在Pd(111)及Pt(111)面上发生了折叠与变形, 环中氢原子向上翘起, 最终通过环平面平行吸附于hollow 位. 此外, 噻吩环吸附后芳香性遭到了破坏, 环中碳原子发生sp³杂化, 同时电子逐渐由噻吩向M(111)面发生转移, M(111)面上的部分电子也反馈给了噻吩环中的空轨道, 这种协同作用最终导致了噻吩分子稳定吸附于M(111)面.     

噻吩在M(111)(M=Pd,Pt,Au)表面的吸附(英文)

采用密度泛函理论(DFT),选取DMol3程序模块,对噻吩在M(111)(M=Pd,Pt,Au)表面上的吸附行为进行了探讨.通过对噻吩在不同底物金属上的吸附能、吸附构型、Mulliken电荷布居、差分电荷密度以及态密度的分析发现,噻吩在Pd(111)面上的吸附能最大,Pt(111)面次之,Au(111)面最小.吸附后,噻吩在Au(111)面上的构型几乎保持不变,最终通过S端倾斜吸附于top位;噻吩在Pd(111)及Pt(111)面上发生了折叠与变形,环中氢原子向上翘起,最终通过环平面平行吸附于hollow位.此外,噻吩环吸附后芳香性遭到了破坏,环中碳原子发生sp3杂化,同时电子逐渐由噻吩向M(111)面发生转移,M(111)面上的部分电子也反馈给了噻吩环中的空轨道,这种协同作用最终导致了噻吩分子稳定吸附于M(111)面.     

Physics of artificial nano-structures on surfaces

Quantum electron transport through nano-structures such as metal atomic wires or molecular bridges is investigated with various theoretical approaches. The difference of the quantization feature between Na and Al atom wires is explained based on the eigenchannel analyses combined with the recursion-transfer matrix calculation. The eigenchannels are calculated self-consistently for Au atom wires at finite bias voltage and the nonlinear conductance is explored in relation to the offset energies of d band channels. As for molecular bridges, we find that a remarkable metalization is caused, if the coupling of the molecule with the metal electrode is enhanced. Internal current distribution within the molecular networks is discussed and exotic properties of the quantum transport is found. In particular, a strong induced loop current is revealed circulating the ring part of the molecule. The direction of the loop current is switched sharply when the electron incident energy sweeps a degenerate molecular level.     

Au(111)表面Verdazyl自由基的构象转换

Pure organic radical molecules on metal surfaces are of great significance in exploration of the electron spin behavior. However, only a few of them are investigated in surface studies due to their poor thermal stability. The adsorption and conformational switching of two verdazyl radical molecules, namely, 1, 5-biisopropyl-3-(benzo[b]benzo[4,5]thieno[2, 3-d]thiophen-2-yl)-6-oxoverdazyl (B2P) and 1, 5-biisopropyl-3-(benzo[b]benzo[4,5]thieno[2, 3-d]thiophen-4-yl)-6-oxoverdazyl (B4P), are studied by scanning tunneling microscopy (STM) and density functional theory (DFT). The adsorbed B2P molecules on Au(111) form dimers, trimers and tetramers without any ordered assembly structure in which two distinct appearances of B2P in STM images are observed and assigned to be its "P" and "T" conformations. The "P" conformation molecules appear in the STM image with a large elliptical protrusion and two small ones of equal size, while the "T" ones appear with a large protrusion and two small ones of different size. Likewise, the B4P molecules on Au(111) form dimers at low coverage, strip structure at medium coverage and assembled structure at high coverage which also consists of above-mentioned two conformations. Both B2P molecules and B4P molecules are held together by weak intermolecular interaction rather than chemical bond. STM tip induced conformational switching of both verdayzl radicals is observed at the bias voltage of +2.0 V. The "T" conformation of B2P can be switched to the "P" while the "P" conformation of B4P can be switched to the "T" one. For both molecules, such a conformational switching is irreversible. The DFT calculations with Perdew-Burke-Ernzerhof version exchange-correlation functional are used to optimize the model structure and simulate the STM images. STM images of several possible molecular conformations with different isopropyl orientation and different tilt angle between verdazyl radical and Au(111) surface are simulated. For conformations with different isopropyl orientation, the STM simulated images are similar, while different tilt angles of verdazyl radical lead to significantly different STM simulated images. Combined STM experiments and DFT simulations reveal that the conformational switching originates from the change of tilting angle between the verdazyl radical and Au(111) surface. The tilt angles in "P" and "T" conformations are 0° and 50°, respectively. In this study, two different adsorption conformations of verdazyl radicals on the Au(111) surface are presented and their exact adsorption structures are identified. This study provides a possible way to study the relationship between the electron spin and configuration conversion of pure organic radical molecules and a reference for designing more conformational switchable radical molecules that can be employed as interesting molecular switches.     

Electronic and structural properties of oligophenylene ethynylenes on Au(111) surfaces

The interaction of oligophenylene ethynylene (OPE) on the (111) surface of a gold slab resembling a self-assembled monolayer (SAM) is investigated using ab initio density functional theory calculations. The authors performed a full optimization including all atoms in the OPE and in the slab to better understand OPE adsorption on the surface. It is found that OPE has two energetically favorable adsorption sites on the Au surface with relatively different molecular geometries: the nontop site adsorption greatly modifies the (111) surface structure; however, the extensive electron interactions enable a delocalized electron density distribution, implying an improved conductivity between OPE and Au, and the top site which is 0.9 eV higher in energy than the nontop and features weaker Au-S bonds. Interestingly the on top configuration shows a strong spin imbalance along the molecule and the nontop shows a small spin imbalance on the surface. This feature is of strong interest for the development of resonators for the detection of chemical and biological agents. They have also calculated the frequency spectrum of these SAMs, which yield deformations in the gold surface yielding peak frequency shifts specific to each absorption site.     

Structure and Electron Transport at Metal Atomic Junctions Doped with Dichloroethylene

We have investigated the structure and electron transport at dichloroethylene-doped metal atomic junctions at low temperatures (20 K) in ultra-high vacuum, using Fe, Ni, Pd, Cu, Ag, and Au. The metal atomic junctions were fabricated using the mechanically controllable break junction technique. After introducing the dichloroethylene (DCE), the conductance behavior of Fe, Ni, and Pd junctions was considerably changed, whereas little change was observed for Cu, Ag, and Au. For the Pd and Cu junctions, a clear peak was observed in their conductance histograms, showing that the single-molecule junction was selectively formed. To investigate the structure of the metal atomic junctions further, their plateau lengths were analyzed. The length analysis revealed that the Au atomic wire was elongated, and the metal atomic wires were formed for the other transition metals: those that do not normally form metal atomic wires without DCE doping, as DCE adsorption stabilized the metal atomic states. There is a strong interaction between DCE and the metals, where DCE supports the formation of the metal atomic wire for Fe, Ni, and Pd.     

Effect of Asymmetric Anchoring Groups on Electronic Transport in Hybrid Metal/Molecule/Graphene Single Molecule Junctions

A combined experimental and theoretical study on molecular junctions with asymmetry in both the electrode type and in the anchoring group type is presented. A scanning tunnelling microscope is used to create the “asymmetric” Au-S-(CH₂)n-COOH-graphene molecular junctions and determine their conductance. The measurements are combined with electron transport calculations based on density functional theory (DFT) to analyze the electrical conductance and its length attenuation factor from a series of junctions of different molecular length (n). These results show an unexpected trend with a rather high conductance and a smaller attenuation factor for the Au-S-(CH₂)_n-COOH-graphene configuration compared to the equivalent junction with the “symmetrical” COOH contacting using the HOOC-(CH₂)_n-COOH series. Owing to the effect of the graphene electrode, the attenuation factor is also smaller than the one obtained for Au/Au electrodes. These results are interpreted through the relative molecule/electrode couplings and molecular level alignments as determined with DFT calculations. In an asymmetric junction, the electrical current flows through the less resistive conductance channel, similarly to what is observed in the macroscopic regime.     

Site- and orientation-selective anchoring of a prototypical molecular building block

The controlled anchoring of molecular building blocks on appropriate templates is a major prerequisite for the rational design and fabrication of supramolecular architectures on surfaces. We report on a particularly selective adsorption process of hexa-peri-hexabenzocoronene on Au(111), which leads to well-controlled adsorption position and orientation of the polycyclic aromatic hydrocarbons. Scanning tunneling microscopy reveals selective adsorption on monatomic steps in the fcc stacking regions with a specific orientation of 18 degrees between the molecular axis and the step normal. Ab initio calculations for various adsorption sites reveal the lowest total energy for adsorption on a kink site. Energy considerations and the excellent agreement between experimental and simulated images show that adsorption on kink sites is responsible for the specific adsorption angle.     

Self-assembled monolayers of alkaneselenolates on (111) gold and silver

Self-assembled monolayers (SAMs) formed from didodecyl diselenide (C12SeSeC12) and didodecyl selenide (C12SeC12) on (111) Au and Ag substrates were extensively characterized by several complementary techniques. C12SeSeC12 was found to form contamination-free, densely packed, and well-ordered C12Se SAMs on both substrates, whereas the adsorption of C12SeC12 occurred only on Au and resulted in the formation of a SAM-like C12SeC12 film with a low packing density and a conformational disorder. The properties of the C12Se SAMs were compared with those of dodecanethiolate (C12S) SAMs. The packing density, orientational order, and molecular inclination in C12Se/Au and C12S/Au were found to be very similar. In contrast, C12Se/Ag exhibited significantly lower packing density, a lower degree of the conformational and orientational order, and a larger molecular inclination than C12S/Ag. The results suggest a sp3 bonding configuration for the selenium atom on Au and Ag and indicate a larger corrugation of the headgroup-substrate binding energy surface in C12Se/Ag than in C12S/Ag.     

Theoretical study of CO adsorption on gold/alumina substrates

Aiming to understand the role of the substrate in the adsorption of carbon monoxide on gold clusters supported on metal-oxides, we have started a study of that process on two different alumina substrates: an amorphous-like fully relaxed stoichiometric (Al2O3)20 cluster and the Al terminated (0001) surface of alpha-(Al2O3) crystal. In this paper, we present first principles calculations for the adsorption of one Au atom on both alumina substrate and the adsorption of Au8 on (Al2O3)20. Then, we study the CO adsorption on the minimum energy structure of these three different gold/alumina systems. A single Au adsorbs preferably on top of an Al atom with low coordination, the binding energy being higher in the case of Au/(Al2O3)20. CO absorbs preferably on top of the Au atom, but in the case of Au/(Al2O3)20, Au forms a bridge with the Al and O substrate atoms after CO adsorption. We find other stable sites for CO adsorption on the cluster but not on the surface. This result suggests that the Au activity toward CO may be larger for the amorphous cluster than for the crystal surface substrate. For the most stable Au8/(Al2O3)20 configuration, two Au atoms bind to Al and a O atoms respectively and CO adsorbs on top of the Au which binds to the Al atom. We find other CO adsorption sites on supported Au8 which are not stable for the free Au8 cluster.     

An atomistic model and key parameters for devising single molecular nanowire sensors

Based on one impurity model Hamiltonian describing a nanowire upon adsorption of a molecule, we obtain an analytical formula of the conductance which is governed clearly by modulating key parameters. The formula shows that the conductance change in nanowire upon adsorption of a molecule is mainly controlled by three factors, electron hopping between adsorbed molecule and nanowire, chemical potential, and the change of atomic configurations of the nanowires near the adsorption site. Conductance is very sensitive to the choice of these key parameters; therefore, a proper nanowire system that renders matched chemical potential as well as hopping strength between the nanowire and the adsorbed molecule should be devised for the sensor applications. Our model calculations give similar conductance features to the conductance obtained by the first principle calculations for a singe-molecule-adsorbed molecular wire. It is worthy of note that the system can be in antiresonance, which is characterized by a quick drop in conductance when a molecule is adsorbed on the nanowires.     

Robust conductance of dumbbell molecular junctions with fullerene anchoring groups

The conductance of a molecular wire connected to metallic electrodes is known to be sensitive to the atomic structure of the molecule-metal contact. This contact is to a large extent determined by the anchoring group linking the molecular wire to the metal. It has been found experimentally that a dumbbell construction with C(60) molecules acting as anchors yields more well-defined conductances as compared to the widely used thiol anchoring groups. Here, we use density functional theory to investigate the electronic properties of this dumbbell construction. The conductance is found to be stable against variations in the detailed bonding geometry and in good agreement with the experimental value of G=3×10(-4) G(0). Electron tunneling across the molecular bridge occurs via the lowest unoccupied orbitals of C(60) which are pinned close to the Fermi energy due to partial charge transfer. Our findings support the original motivation to achieve conductance values more stable towards changes in the structure of the molecule-metal contact leading to larger reproducibility in experiments.     

Nonequilibrium electronic transport of 4,4'-bipyridine molecular junction

The electronic transport properties of a 4,4'-bipyridine molecule sandwiched between two Au(111) surfaces are studied with a fully self-consistent nonequilibrium Green's-function method combined with the density-functional theory. The 4,4'-bipyridine molecule prefers to adsorb near the hollow site of the Au(111) surface and distorts slightly. The modifications on the electronic structure of the molecule due to the presence of the electrodes are described by the renormalized molecular orbitals, which correspond well to the calculated transmission peaks. The average Fermi level lies close to the lowest unoccupied renormalized molecular orbital, which determines the electronic transport property of the molecular junction under a small bias voltage. The total transmission is contributed by a single channel. The transmission peaks shift with the applied bias voltage, and this behavior depends on the spatial distribution of the renormalized molecular orbitals and the voltage drop along the molecular junction. The shape of the calculated conductance curve of the equilibrium geometric configuration reproduces the main feature of the experimental results, but the value is larger than the measured data by about 6 times. Good agreement with the experimental measurements can be obtained by elongating the molecular junction. The electronic transport behaviors depend strongly on the interface configuration.     

Effect of anchoring groups on single-molecule conductance: comparative study of thiol-, amine-, and carboxylic-acid-terminated molecules

We studied the effect of anchoring groups on the conductance of single molecules using alkanes terminated with dithiol, diamine, and dicarboxylic-acid groups as a model system. We created a large number of molecular junctions mechanically and analyzed the statistical distributions of the conductance values of the molecular junctions. Multiple sets of conductance values were found in each case. The I-V characteristics, temperature independence, and exponential decay of the conductance with the molecular length all indicate tunneling as the conduction mechanism for these molecules. The prefactor of the exponential decay function, which reflects the contact resistance, is highly sensitive to the anchoring group, and the decay constant is weakly dependent on the anchoring group. These observations are attributed to different electronic couplings between the molecules and the electrodes and alignments of the molecular energy levels relative to the Fermi energy level of the electrodes introduced by different anchoring groups. For diamine and dicarboxylic-acid groups, the conductance values are sensitive to pH due to protonation and deprotonation of the anchoring groups. Further insight into the binding strengths of these anchoring groups to gold electrodes is obtained by statistically analyzing the stretching length of molecular junctions.     

Conductivity in alkylamine/gold and alkanethiol/gold molecular junctions measured in molecule/nanoparticle/molecule bridges and conducting probe structures

Charge transport through alkane monolayers on gold is measured as a function of molecule length in a controlled ambient using a metal/molecule/nanoparticle bridge structure and compared for both thiol and amine molecular end groups. The current through molecules with an amine/gold junction is observed to be more than a factor of 10 larger than that measured in similar molecules with thiol/gold linkages. Conducting probe atomic force microscopy is also used to characterize the same monolayer systems, and the results are quantitatively consistent with those found in the nanoparticle bridge geometry. Scaling of the current with contact area is used to estimate that approximately 100 molecules are probed in the nanoparticle bridge geometry. For both molecular end groups, the room-temperature conductivity at low bias as a function of molecule length shows a reasonable fit to models of coherent nonresonant charge tunneling. The different conductivity is ascribed to differences in charge transfer and wave function mixing at the metal/molecule contact, including possible effects of amine group oxidation and molecular conformation. For the amine/Au contact, the nitrogen lone pair interaction with the gold results in a hybrid wave function directed along the molecule bond axis, whereas the thiol/Au contact leads to a more localized wave function.