Mohammed
Noori‡
ab,
Hatef
Sadeghi
a and
Colin J.
Lambert
*a
aQuantum Technology Centre, Department of Physics, Lancaster University, Lancaster LA1 4YB, UK. E-mail: h.sadeghi@lancaster.ac.uk; c.lambert@lancaster.ac.uk
bDepartment of Physics, College of Science, Thi-Qar University, Thi-Qar, Iraq
First published on 31st March 2017
If high efficiency organic thermoelectric materials could be identified, then these would open the way to a range of energy harvesting technologies and Peltier coolers using flexible and transparent thin-film materials. We have compared the thermoelectric properties of three zinc porphyrin (ZnP) dimers and a ZnP monomer and found that the “edge-over-edge” dimer formed from stacked ZnP rings possesses a high electrical conductance, negligible phonon thermal conductance and a high Seebeck coefficient of the order of 300 μV K−1. These combine to yield a predicted room-temperature figure of merit of ZT ≈ 4, which is the highest room-temperature ZT ever reported for a single organic molecule. This high value of ZT is a consequence of the low phonon thermal conductance arising from the stacked nature of the porphyrin rings, which hinders phonon transport through the edge-over-edge molecule and enhances the Seebeck coefficient.
The efficiency of a thermoelectric (TE) material or device is determined by the dimensionless thermoelectric figure of merit ZT = GS2T/κ, where G is the electrical conductance, T is temperature, S is the thermopower (Seebeck coefficient) and κ = κel + κph is the thermal conductance due to electrons (κel) and phonons (κph). The Seebeck coefficient characterizes the ability of a thermoelectric material to convert heat to electricity and is defined as S = −ΔV/ΔT, where ΔV is the voltage difference generated between the two ends of the junction when a temperature difference ΔT is established between them.7,16–19 Enhancing the efficiency of TE materials is not easy, because all parameters are correlated. For example at the fundamental level, the electronic properties G, S and κel are related, because they are all derived from the transmission coefficient Tel(E) describing electrons of energy E passing from one electrode to other through a molecule (see Methods). In particular the Seebeck coefficient S is approximately proportional to the slope of lnTel(E), evaluated at the Fermi energy EF, whereas the electrical conductance is proportional to Tel(EF). Therefore, if the Fermi energy lies in a region of high slope, close to a transmission resonance then both G and S are enhanced.20 On the other hand, to decrease the thermal conductance κ, which appears in the denominator of ZT, both electron and phonon transport must be engineered. Therefore, simultaneous consideration of both electron and phonon transport is needed to develop new materials for thermoelectricity.
Since only a few groups worldwide are able to measure the thermal conductance of single molecules, theoretical investigation is needed to identify new strategies to simultaneously suppress phonons and enhance S and G. Recent proposals to reduce phonon transport in molecular junctions include weakening the overlap between the continuum of vibrational states in the electrodes and discrete vibrational modes of the molecules,21 taking advantage of the weak interaction between different parts of the molecules, as in π–π stacked structures19 and using the low Debye frequency of electrodes to filter high-frequency phonons.20 The challenge is to identify new materials and device structures in which such strategies can be realized in the laboratory. In this paper, we present a comparative theoretical study of the thermoelectric properties of four different zinc porphyrin structures and elucidate a new strategy for simultaneously increasing their thermopower and reducing their thermal conductance leading to a high value of ZT.
To calculate the electronic properties of the molecules in the junction, from the converged DFT calculation, the underlying mean-field Hamiltonian H was combined with our quantum transport code, GOLLUM29 to calculate the transmission coefficient Tel(E) for electrons of energy E passing from the source to the drain. The electrical conductance Gel(T) = G0L0, the electronic contribution of the thermal conductance κel(T) = (L0L2 − L12)/hTL0 = and the thermopower S(T) = −L1/eTL0 of the junction are calculated from the electron transmission coefficient Tel(E) where and f(E,T) is the Fermi–Dirac probability distribution function f(E,T) = (e(E−EF)/KBT + 1)−1, T is the temperature, EF is the Fermi energy, G0 = 2e2/h is the conductance quantum, e is electron charge, and h is the Planck's constant.
For the structures of Fig. 1, Fig. 2 show the transmission coefficients for electrons with energy E and phonons of energy ħω, passing through a molecule from the left electrode to the right electrode, calculated using the method described in ref. 20. We first carry out geometrical optimization of each molecule placed between two gold electrodes using the SIESTA28 implementation of density functional theory (DFT) to find the ground state optimized positions of the atoms relative to each other (see Methods). From the ground state geometry, we obtain the mean-field Hamiltonian of each system comprising both electrodes and molecules and use our transport code GOLLUM29 to calculate the transmission coefficients Tel(E) (see Methods). In each case the optimal angle between the porphyrins is zero, which corresponds to the maximum conductance that could be obtained.22 The electronic transport properties of 1 and 3 have been studied experimentally in the literature,22 so we used these to benchmark our calculations. As shown in Table S1 of the ESI,† our calculated conductances for these molecules are in good agreement with experiment. The electron transmission of 4 is much smaller than 1, 2 and 3, whereas the transmission of 2 is either equal to that of 3 near the HOMO resonance or lower in the vicinity of the middle of the HOMO–LUMO gap. As shown in Fig. (SI2†), this is reflected in the electrical conductance as a function of temperature.
To calculate the vibrational properties of each structure, we use the harmonic approximation to construct the dynamical matrix D. Each atom is displaced from its ground-state equilibrium position by δq′ and −δq′ in the x, y, and z directions and the forces on all atoms calculated in each case. For 3n degrees of freedom (n = number of atoms), the 3n × 3n dynamical matrix Dij = (Fqi(δq′j) − Fqj(−δq′j))/2Mijδq′j is constructed, where F and M are the force and mass matrices (see ref. 20 for more details). For an isolated molecule, the square root of the given values of D determines the frequencies ω associated with the vibrational modes of the molecule in the junction (see the ESI†). For a molecule within a junction, the dynamical matrix describes an open system composed of the molecule and two semi-infinite electrodes and is used to calculate the transmission coefficient Tph(ω) for phonons with energy ħω passing through the molecule from the right to the left electrodes.20
Fig. 2b–e shows Tph(ω) for the four structures of Fig. 1. It is apparent that the widths of the resonances in the edge-over-edge ZnP-dimer 2 are smaller than those of the other structures and the low energy phonons (in the range 2–5 meV) are either suppressed or pushed to higher frequencies. This can be demonstrated using the participation ratio of the dimer molecular cores 2, 3 and 4 and comparing the integrated density of states N(ħω) of 2, 3 and 4. As shown in Fig. S1 of the ESI,† the participation ratio of the molecule core (ZnPs and linkers) connected to the gold surface is mostly due to the in-plane (PRy) and out of plane (PRx) transverse modes in structures 3 and 4, whereas out-of plane transverse modes are mainly suppressed or converted to in-plane transverse modes and moved to the higher frequency, reflecting the higher rigidity of the edge-over-edge structure. In addition, the integrated density of states is almost the same for 3 and 4, whereas for low frequencies, the integrated density of states of 2 is smaller than 3 and 4. This means the thermal conductance is reduced significantly in 2, because transmission of the low energy modes is suppressed due to the scattering from in-plane modes to cross-plane transverse modes. In addition, some modes are pushed to higher frequency, although this is smaller effect compared with the suppression of low frequency transmission. Overall, these two effects combine to yield a lower phonon thermal conductance in 2.
The thermal conductance of the junction (κ = κph + κel) is obtained by summing the contributions from both electrons (κel) and phonons (κph). The electronic (phonon) thermal conductances are calculated from the electron (phonon) transmission coefficients shown in Fig. 2a–e. Fig. 3a shows that the ZnP monomer 1 has the lowest value of κph while 4 has the highest. This is counter-intuitive, because one would expect a higher thermal conductance for shorter molecules. However, due to the more rigid nature of the monomer, its vibrational modes are pushed to higher frequencies and therefore their contribution to the room temperature conductance is suppressed. In addition, Fig. 3b and c show that the thermal conductance due to the electrons κel of the dimer ZnP 3 is higher than those of the edge-over-edge ZnP and structures 1 and 4 for a wide range of energy in the vicinity of DFT predicted Fermi energy. The crucial point is that almost for all Fermi energies, the electronic contribution to the thermal conductance is higher than the phonon contribution. This is significant, because to achieve a high-ZT material, one needs to only focus on engineering the electronic properties of structure 2.
To examine the thermoelectric properties of 1–4, we obtained the Seebeck coefficient of all structures from the electron transmission coefficient Tel(E), as described in the Methods. Fig. 4a shows the Seebeck coefficients as a function of Fermi energy EF (and also as a function of temperature in Fig. S3 of the ESI†) and reveals that the edge-over-edge ZnP dimer 2 has a higher Seebeck coefficient than 1, 3 and 4 due to the higher slope of lnTel(EF) over a wide range of Fermi energies between the HOMO and LUMO. Since the electronic contribution to the thermal conductance is higher in 1, 2 and 3, the contribution of the phonons is negligible. Furthermore the electrical conductance is proportional to the electronic thermal conductance, so they cancel each other in ZT. Consequently as shown in Fig. 4b, due to the high Seebeck coefficient of the edge-over-edge dimer, a ZT as high as ≈4 is obtained when EF lies in a wide energy window in the vicinity of the DFT-predicted Fermi energy. Fig. 4b also shows that the less-rigid structure 4 is not promising for efficient conversion of the heat to electricity. Although all of these structures are made from ZnP, this study shows the importance of the molecular design. The more rigid edge-over-edge ZnP dimer 2 shows a very high ZT, whereas the less conductive structure 3 is unattractive for thermoelectricity.
To our knowledge, there currently exist no measurements of the Seebeck coefficient of porphyrins and no measurements of their ZT. The Seebeck coefficient of n-alkanedithiol (ADT) (n = 2, 3, 4, 5, 6, 8) was found to range from 6.8 to 2.4 μV K−1, depending on the length n,32,33 whereas 1,4-benzenedithiol (BDT), 4,4′-dibenzenedithiol (DBDT) and 4,4′′-tribenzenedithiol (TBDT) attached to Au electrodes have been measured by several groups34–41 and found to increase with the number of phenyl rings, ranging from 7 μV K−1 to 16 μV K−1, by changing the terminal groups, measurements on 1,4-bis((trimethylstannyl)methyl)-n-phenyl (n = 1, 2, 3, 4) revealed that the Seebeck coefficient increased to 24 μV K−1 for the longest molecules. On the other hand 1,4-n-phenylenediamine (PDA) (n = 1, 2, 3) was found to have lower values in the range of 2.1–10.4 μV K−1.32,42,43 These values are comparable to the Seebeck coefficient of oligothiophenes, which can be as high as 14.8 μV K−1 for ethanethioate-terminated terthiophenes on gold.44
Changing the terminal groups to pyridine generally switches transport towards the LUMO and changes the sign of the Seebeck coefficient. For example the Seebeck coefficients of 4,4′-bipyridine and 1,2-di(4-pyridyl)ethylene attached to gold electrodes were found to be approximately −9 μV K−1 and −10 μV K−1 respectively.43,45
Fullerenes tend to have higher Seebeck coefficients, ranging from −10 to −30 μV K−1 for C60 with different metallic electrodes46 to −33 μV K−1 for C60 dimers47 and up to −31.6 μV K−1 C82 endohedral fullerenes.48 The sign of the endohedral fullerene Sc3N@C80 was shown to be sensitive to pressure, ranging from −25 μV K−1 to +25 μV K−1, depending on the orientation of the molecule on a gold substrate.49
The length dependent thermal conductance in Au-alkanethiol-SiO2 junctions by the scanning thermal microscopy (SThM) technique has been measured to vary from 25 pW K−1 for the shortest chains to 10 pW K−1 for chains of 18 carbon atoms.50 The thermal conductance of Au-alkanedithol (n = 8, 9, 10)-GaAs junctions at room temperature51 was measured to be approximately 27 MW m−2 K−1, assuming a 1 nm2 footprint per molecule, this equates to 27 pW K−1 per molecule. Values ranging from 35 to 65 MW m−2 K−1 were measured for different terminal groups.52
Single-molecule thermoelectricity is a rapidly expanding field of fundamental research and our paper is aimed primarily at further developing our understanding at the single-molecule level. Acquisition of the fundamental understanding of single-molecule thermoelectricity will underpin the design and synthesis of high-performance molecular films, but will not be sufficient, because a range of additional issues arise when single molecules are placed in parallel to form a molecular film. Ref. 53 discusses the question of how a statistical ensemble of single-molecule conductances and Seebeck coefficients combine in parallel to yield the thermoelectric properties of a dilute molecular film and concludes that a conductance-weighted single-molecule Seebeck coefficient is the relevant quantity. For more dense films, thermoelectric properties may be modified by intermolecular interactions, although as mentioned above, measurements of single-molecule thermal conductances of alkanes50 are comparable to the thermal conductance per molecule obtained from the measurements on molecular films51,52 and therefore at least for alkanes, such interactions do not appear to be significant.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c6nr09598d |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2017 |