Robert M. Metzger
Laboratory of Molecular Electronics, Department of Chemistry, University of Alabama, Tuscaloosa, AL 35487-0336, USA, rmetzger@bama.ua.edu
First published on UnassignedUnassigned22nd December 1999
We found unimolecular rectification in a single molecule, (N-hexadecylquinolin-4-ium-1-yl)methylidenetricyanoquinomethanide, 1, in which the ground state is zwitterionic: D+–π–A−, and has a large dipole moment (43 D) while the first excited state is undissociated: D0–π–A0 (and has a smaller moment, between 3 and 9 D). The intervalence optical absorption band connecting these two states is strongly hypsochromic. Spectroscopic measurements (FTIR, NMR, VIS–UV, XPS) all confirm this assignment. This significant ‘first step' to unimolecular electronics, i.e. a rectifier only 2.3 nm thick, confirms a 1974 proposal by Aviram and Ratner. However, the ‘real proof' that the large current under forward bias actually travels through the molecule, has not yet been obtained.
By now some important milestones towards that goal have been reached:
(1) tunneling currents through aliphatic chains are larger than through aromatic chains,5
(2) the resistance of a single molecule (benzene-1,4-dithiol) bonded to two Au electrodes was measured as several MΩ (it would have been smaller, but the work function of Au and the LUMO of the molecule were mismatched, and the thiolate bonds are partially ionic),6
(3) the quantum of electrical resistance (12 kΩ) was measured at room temperature between a single carbon nanotube, glued to a conducting AFM tip, and liquid Hg,7
(4) the Aviram–Ratner mechanism,1 slightly modified, was confirmed in conductivity measurements through a monolayer of (Z)-[β-(N-hexadecylquinolin-4-ium-1-yl)-α-cyano-4-styryl]dicyanomethanide or [(N-hexadecylquinolin-4-ium-1-yl)methylidene]tricyanoquinodimethanide, C16H33–Q3CNQ, 1: this is the first proven two-terminal unimolecular electronic device.8
We review here what we have learned about the chemistry and physics of 1.8–16
![]() | ||
Fig. 1 The Aviram–Ratner Ansatz, showing the through-molecule electron flow from the excited zwitterion state D+–σ–A− to the undissociated ground state D0–σ–A0. |
T–D+–π–A− → T–D0–π–A0
In 1, the quinolinium and tricyanoquinodimethanide rings may be twisted by some unknown angle θ. If θ = 0°, then the two states would be degenerate resonance states. If θ = 90° then IVT vanishes, and the blue color should disappear.
When LB multilayers of 1 were sandwiched between a noble metal electrode (Pt or Ag) and a Mg electrode (protected by an overlayer of Ag), then asymmetric currents through the laterally macroscopic multilayer indicated electrical rectification through the molecule:27
Pt|LB multilayer of 1|Mg|Ag
and even in the sandwich:
Pt|LB monolayer of 1|Mg|Ag
Doubts raised about a possible Schottky barrier at the Mg electrode|molecule interface20 prompted the insertion of an insulating monolayer of a fatty acid between the metal electrodes and three monolayers of C16H33–Q3CNQ;28I–V asymmetries were seen again, albeit between metal electrodes of very different work functions, and were ascribed to molecular rectification.28,29
The NMR spectrum of the compound8 is very diagnostic of a dipolar ground state. Since it was opined that 1 may undergo a twisted internal charge transfer (TICT) transition, maybe as a function of temperature,30 a temperature-dependent NMR study was undertaken, but no dramatic changes in the NMR spectrum were found.12
The dipole moment of a solution of 1 in CH2Cl2 was studied as a function of both temperature and concentration. After correcting for the polarization of the solvent, two estimates emerged (one from the temperature dependence of the relative permittivity, which extrapolated to about 47 D at zero concentration, the other using the Kirkwood–Fröhlich equation, which extrapolated to 40 D), whose average is 43 ± 8 D for 1 at infinite dilution in CH2Cl2: this result strongly supports the zwitterionic nature of the ground state of 1.8
There are two sets of fluorescent emission lines: emissions below 600 nm, and emissions in the 1000–1400 nm region. The emissions in the near IR region are very weakly hypsochromic. That two fluorescent emissions exist for 1 is an apparent violation of Kasha's rule, which dictates that, within a set of states of the same spin multiplicity, a fluorescent emission should occur only from the lowest excited state. Other, higher excited states usually decay internally and without radiative emission, into the lowest excited state, which then emits the photon. However, if the higher-energy emission is weak (here ϕ ≈ 0.01), and is not re-absorbed, then emission from a lower-energy state becomes possible. Alternatively, molecule 1 may have different rotameric states: molecules of 1 with θ ≈ 90°, with no IVT, would emit in the visible, while molecules with a lower θ and a ‘blue' IVT will emit in the near infrared region.
An estimate of the excited state dipole moment µe can be made by comparing the absorption and fluorescent emission energies, assuming no change in the orientation of the dipole: the ratio of the energy shift between any two solvents for absorption ΔEa and for fluorescence ΔEf is ΔEa/ΔEf = μg/μe. If μg = 43 D,1 one obtains μe = 5.53 D (CHCl3–CH3CN) or μe = 2.97 D (CH2Cl2–CH3CN).12 The solvatochromism, explained by an ellipsoidal cavity model,31,32 yields an estimated excited-state dipole moment of 8.7 D.12
Further calculations by density functional theory, AM1, and ab initio methods15 found new conformers with high dipole moments (which increase if a solvent is included). The twist angles θ were larger than expected (Fig. 2): the ‘true' energy minimum that explains all experimental data may not yet have been found.15
![]() | ||
Fig. 2 Conformational search (with and without CH3CN solvent) for high-moment energy minimum of CH3–Q3CNQ.15 |
![]() | ||
Fig. 3 Cyclic voltammogram of a 10−4 M solution of 1 in CH2Cl2, measured using a Pt wire working electrode, an SCE reference electrode, 0.1 M NBu4ClO4 electrolyte, N2 gas bubbled through the solution, and a scan rate of 20 mV s−1.8 |
A saturated solution of 1 in a 0.1 M Bu4NPF6–dimethylformamide solution, held in an electrochemical EPR cell at room temperature under Ar at −0.645 V vs. SCE, where the molecule should be reduced to its radical anion, shows the EPR spectrum shown in Fig. 4.12 This spectrum was simulated, using g = 2.0027, and isotropic hyperfine constants aN = 1.20 G (2 ×
14N, I = 1) and 2.80G (14N), aH = 1.44 G (2 ×
1H, I = 1/2), 1.80 G (2 ×
1H), and 3.0 G (1H) (Fig. 4). The coupling constants are very close to the data for the radical anion of TCNQ (aN = 1.02, aH = 1.44 G), and thus testify that most of the spin density in the radical anion of 1 is localized on the 3CNQ part of the molecule and the bridgehead proton.12
![]() | ||
Fig. 4 Electron paramagnetic resonance spectrum of the radical anion of 1, electrogenerated in DMF at −0.645 V vs. SCE, with 0.1 M Bu4NPF6. The simulated spectrum is shown as dashed lines.12 |
![]() | ||
Fig. 5 Pressure–area isotherms for 1 at 14![]() ![]() ![]() |
Monolayers and multilayers of 1 were transferred as LB films at 15 mm min−1 and 20 mN m−1. Molecule 1 first transfers on the downstroke onto hydrophobic HOPG (transfer ratio 58%), but on the upstroke onto a freshly cleaned hydrophilic Au(111) single crystal and on an Al electrode (transfer ratio 90%).8 At the lower transfer speed of 15 mm min−1 the LB multilayers (past the first monolayer) are Z-type, while at the faster transfer speed of 45 mm min−1, the LB multilayers become Y-type.8 The thickness per monolayer, measured by ellipsometry (λ = 670 nm) and by surface plasmon resonance, is 22 ± 2 Å; thus, the approximate molecular axis may be tilted 48 ± 5° from the normal to the monolayer plane.8 X-Ray diffraction spectra, taken several days after film transfer onto a glass substrate, show that the Y-type films reorder to Z-type after transfer, with a thickness of 23 Å per monolayer.8
The geometry of LB multilayers of 1 on HOPG is:
HOPG|→ ← ← ← ← ← ← ← etc.
On Al it is Z-type throughout:
Al|← ← ← ← ← ← ← etc.
Here the arrow → indicates the rough molecular orientation; it points from the hexadecyl tail T and the quinolinium ring D+ (tail) towards the dicyanomethylene head of A− (arrowhead).8
A sharp and narrow absorption band peaked at 570 nm is found for LB films transferred onto a quartz substrate at 14°C in the dark, and measured at 25
°C (Fig. 6). This narrow band disappears if the films are protonated, but returns when they are exposed to ammonia vapor,8 it disappears irreversibly if the films, or the solution, are exposed to bright light and air for some time.
![]() | ||
Fig. 6 Visible spectrum of an 11-layer LB film of 1 transferred onto quartz; the narrow absorbance has a peak at 570 nm (2.17 eV; spectral width at half-maximum = 0.19 eV); the absorbance per monolayer at 633 nm is 0.0021.8 |
An LB monolayer of 1, transferred onto quartz, shows weak absorption maxima at 340 nm and at 565 nm. The same monolayer, excited at 350 nm, produces a clear emission peak at 492 nm, but no emission when excited at 570 nm. The absorption of 15 blue-violet LB monolayers of 1 on quartz did not change between room temperature and 77 K: this argues for no large change in the twist angle θ between these two temperatures.12
Grazing-angle Fourier transform infrared (FTIR) spectra, using p-polarized light for a monolayer of 1 on a glass substrate covered by Al (Fig. 7), show absorbances at 2850 cm−1 (symmetric C–H), 2920 cm−1 (asymmetric –CH2–), 2139 cm−1 (–CN) and 2175 cm−1 (–C
N), indicating that the cyano groups and the methylene groups do not lie in the plane of the monolayer, but are somewhat tilted from it. The CN stretches show evidence of shifts due to partial change: negatively charged C
N groups (2175 cm−1) and a shifted C
N group (2135 cm−1). For a monolayer of 1 on either Au(111) or HOPG, the C–H and –CH2– signals can be seen clearly; a C
N signal at 2175 cm−1 is barely visible for the monolayer on Au(111).8
![]() | ||
Fig. 7 p-Polarized grazing-angle FTIR spectrum (θ = 85°) of LB monolayer of 1 on Al.8 |
Fig. 8 shows the N 1s region of the XPS spectrum of 1 adsorbed on an oxide-covered Si substrate. There are three peaks, which have been fitted using standard techniques. The dominant feature at 400.5 eV is in the range for N in an organic matrix. However, there are features at higher (402.3 eV) and lower (398.8 eV) binding energies, which indicate the charge transfer state of the molecule. The lowest energy feature (398.8 eV) can be associated with nitrogen in the C≡N moiety; the highest binding energy feature (402.3 eV) is in the range associated with ammonium salts.12 This confirms that 1 is zwitterionic.12
![]() | ||
Fig. 8 Nitrogen N 1s region of the X-ray photoelectron spectrum of an LB multilayer of 1 on Si. Also shown are a sloping linear background fit, Gaussian peak fits to simulate the spectrum, and the error in the fit (small circles). N 1s peaks are seen at 402.3 eV, 400.3 eV, and 398.8 eV.12 |
A valence band scan of a multilayer of 1 on Si (Fig. 9) reveals several peaks at binding energies 3 eV to 25 eV. These peaks can be compared with PM3/RHF calculations. The onset of the valence-band XPS at 3.7 eV on Si (−7.8 eV versus vacuum)12 agrees with the calculated HOMO at −7.98 eV.8
![]() | ||
Fig. 9 Valence-band XPS spectrum of a LB multilayer of 1 on Si, with Gaussian line-fits.12 |
![]() | ||
Fig. 10 Orientation of the LB monolayer (a) or multilayer (b) on the glass, quartz, or Si substrate; the electrode (+) for positive bias, and the direction of ‘easy' electron flow for V > 0 are marked. In some cases Ag paste was used instead of the Ga/In eutectic.8 |
![]() | ||
Fig. 11 Conductivity of a sandwich ‘Ga/In eutectic|Al (100 nm)|4 LB monolayers of 1 (Z-type)|Al (100 nm)|Ga/In eutectic'. Scan rate 100 mV s−1. Easier electron flow (I > 0 for V > 0) is from the bottom electrode through the film to the top electrode.8 |
![]() | ||
Fig. 12 Rectification through a single monolayer of 1 sandwiched between Al electrodes (top Al pad area 4.5 mm2, thickness 100 nm), using Ga/In eutectic and Au wires.8 |
The results on both monolayers and multilayers show that the electrons preferentially flow by IVT from 3CNQ (A−) to quinolinium (D+); this is explained by modifying the AR Ansatz for the case that the ground state is zwitterionic, D+–π–A−, and the excited state is undissociated, D0–π–A0.
The junctions that do not rectify have much lower conductivities.16
Many samples, which rectified at room temperature, were measured at decreasing temperature.13 Results could easily be obtained above 250 K. Below 250 K, either open circuits or permanent or intermittent short circuits happened often: only a few reliable results were obtained below 100 K. Fig. 13 shows the I–V plot for a sample at 105 K.
![]() | ||
Fig. 13 DC current I versus DC applied voltage V for a monolayer of 1 at 105 K.13 |
Although the barrier resistances are small for the Au|Ga/In contact, an ohmic contact resistance of about 100 kΩ was measured in a separate experiment for a non-wetting Ga/In|Al2O3|Al|Al2O3|Ga/In couple.13 Happily, this contact resistance is many orders of magnitude smaller than the resistance across the organic monolayer.13
![]() | ||
Fig. 14 STM Micrograph of Langmuir–Schaefer monolayer of 1 on HOPG, viewed with a Pt/Ir nanotip (Nanoscope II, modified to increase current 100 fold). Scan size = 2 nm × 2 nm. For details, see ref. 8. |
![]() | ||
Fig. 15 STM Image of a LB monolayer of 1 on HOPG, with Pt/Ir tip (Nanoscope III). Scan size = 4.5 nm × 4.5 nm.8 |
![]() | ||
Fig. 16 STS I–V curve for an LB film of 15 Z-type monolayers of 1 on HOPG, with Pt/Ir tip. The higher current for V < −1.35 V corresponds to electron flow from HOPG through molecules to the Pt tip.8,9 |
Asymmetric STS spectra are even seen for a monolayer of 1 on HOPG, and for saturated solution of 1 in DMSO over HOPG.8
Using the bias voltages ΔV measured before asymmetric conduction occurs through monolayers and multilayers on various electrodes8ELUMO = 3.0 ± 0.3 eV is estimated. This value is in rough agreement with the electron affinity of p-benzoquinone (2.4 eV), whose reduction potential in solution matches that of 1.8
The onset of valence-band XPS current,12 and a simple AM1 calculation8 both suggest EHOMO = 8 ± 0.3 eV.12 Clearly EHOMO is below the work function of most metals. Then ELUMO − EHOMO = 8.0 − 3.0 = 5 eV, which is probably large. The IVT energy is measured as only 2.17 eV in films8 and estimated to be 1.1 eV in the gas phase:12 these energies are much smaller than 5 eV, but Coulomb correlation in the D+–π–A− state and Mulliken charge-transfer mixing of states may lower the gap between the D+–π–A− and D0–π–A0 states. Clearly more calculations are needed.
We have confirmed that 1 is not only a multilayer rectifier, but also a monolayer rectifier, by an improvement of Sambles and co-workers' original technique,27–29and in accord with a plausible modification of the Aviram–Ratner Ansatz.1 We have eliminated the uncertainty due to the previous use27,28 of two metal electrodes with different work functions.
The observed hypsochromism and large Stokes shift of the IVT fluorescence emission are consistent with the large dipole moment (43 ± 8 D) for the zwitterionic ground D+–π–A− state, and a much smaller moment (between 3 and 9 D) for the undissociated excited state D+–π–A−.
The core-level XPS spectrum confirms several valences for N in the molecule; the valence-band XPS spectra can be correlated very well with the calculated PM3 (AM1) occupied molecular orbitals.
The conspicuous asymmetry seen here for monolayers of 1 between symmetrical electrodes between 1 V and 1.5 V in the temperature range 370 K to 105 K, is most likely due to a molecular process. Since rectification by a single monolayer between symmetrical electrodes has now been observed as a function of temperature, the possibility that rectification is caused by some thermally activated mechanism can be discounted.
An oriented monolayer with a high dipole moment may induce at zero bias polarization charges in the metal electrodes M1 and M2 adjacent to the monolayer, as follows: M1−|D+–π–A−|M2+. However, the observed rectification8 is in a direction opposite to what would be expected from such a polarization model: the modified Aviram–Ratner model pushes electrons (at resonance) from metal to molecule, i.e. the direction of easy electron transfer is when M1 is biased +, and M2 is biased −.8
An earlier22 critique that there may be Schottky barriers in the ‘Pt|monolayer of 1|Mg|Ag' sandwich at the magnesium–molecule interface,27 or, less likely, in ‘Al|1|Al', is probably not valid: we have found that 1 is a weak one-electron acceptor in solution (E1/2 = −0.513 V vs. SCE), on a par with benzoquinone. The gas-phase electron affinity of benzoquinone (1.95 eV) is much smaller than the work function of Mg (3.66 eV). The energetics speak against a salt formation at the ‘Mg|benzoquinone' or, worse, ‘Al|benzoquinone' interface. If the analogy holds, then molecule 1 will not form a Schottky barrier with Al, or even with Mg. The strong acceptor TCNQ, or 7,7,8,8-tetracyanoquinodimethane, whose electron affinity is 3.3 eV, should form such a Schottky barrier quite easily.
There are still some unsolved puzzles:
(1) The sandwiches using Al or Mg electrodes bear an inevitable oxide layer. Al is a ‘valve' metal, and its thin covering with oxide is not defect-free, unless it is anodized. Control experiments using arachidic acid8 reduce the problem, but may not eliminate it. A centrosymmetric Y-type LB multilayer of arachidic acid exhibits no rectification.8 Does an acentric monolayer of 1 somehow change the properties of the Al2O3 between it and the base Al electrode, or does it change the Al2O3 between it and the Al pad on top of the monolayer? This puzzle is difficult to answer at present. It seems unlikely that a physisorbed organic monolayer could interfere with the oxide formation on an electrochemically active metal such as Al. By crafting a symmetrical Al|monolayer|Al geometry, we think that we avoided the problem, but cannot be sure. It would be nice if one could work with an oxide-free substrate. Unfortunately, adhesion of LB films to hydrophilic Au is poor, and depositing oxide-free Au pads on an LB monolayer destroys it by heating, despite cryocooling the sample holder.
(2) The Ga/In eutectic has, typically, a 100 kΩ contact resistance with the Al pads,13 which is several orders of magnitude less than the resistance of the LB monolayer. When the eutectic wets the Al by piercing through the oxide layer, then the pad lifts off the monolayer. Ag paste has similar problems.
(3) The measured current of 0.33 electrons per molecule per second (0.053 aA) through a single monolayer8 is many orders of magnitude lower than the currents measured in an STM experiment (10 pA to 1 na):8 maybe only one molecule in a million is ‘at work'. Is this due to the oxide limiting the conductivity to the defect sites?
(4) The rectification ratio drops upon cycling:8 physisorption clearly is a weak force, compared to the intense electrical fields encountered during measurement (1.5 V/25 Å). This could be eliminated by chemisorbing a suitably modified version of 1 onto Si or Al. An acidic thiol termination is incompatible with 1. A silanized version of molecule 1 was prepared, but did not form a uniform layer on Si, but experiments are continuing.
(5) The Volta potential of about 0.5 V for a monolayer of 1 at the air/water interface8 or for a dry monolayer of 2 on Al is one order of magnitude smaller than expected for a zwitterionic monolayer.
(6) A theoretical calculation of the I–V asymmetry for 1 would be welcome.
Footnote |
† Basis of a presentation given at Materials Chemistry Discussion No. 2, 12–15 September 1999, University of Nottingham, UK. |
This journal is © The Royal Society of Chemistry 2000 |