1 Introduction
Since 1974 when Aviram and Ratner firstly proposed a theoretical model of molecular rectifiers [
1], molecular electronics have gradually become a cause for concern. Single-molecule junctions show the numerous innovative properties that may be applicable in new forms of electronic devices, such as negative differential resistance, shot noise, and Fano resonances [
2-
4]. Experimentally, the current–voltage (
I–
V) characteristics of individual molecules are measured within the mechanically-controllable break junctions or the scanning tunneling microscopes [
5,
6]. However, the critical aspects about the electronic structure generally remains unknown by the
I–
V characteristics alone. Until recently, with the advancement of experimental measurements, not only the electrical conductance (
) but also the thermopower can be obtained in a modified scanning tunneling microscope [
7]. The thermopower can indicate whether single-molecule junctions are p-type or n-type [
8].
The thermopower is also called the Seebeck coefficient (S) which, in a linear response regime, can be defined as
is the voltage counterbalancing the electronic current initiated by a temperature difference (
) across single-molecule junctions. For instance, the measured thermopower of C
60-based single-molecule junctions are with a mean value of -18 μV/K by using a scanning tunneling microscope and ~-14.5 μV/K within the mechanically-controlled break junction [
9,
10]. In addition, the “-” sign of the Seebeck coefficient implies that the position of the chemical potential of the metal electrodes is closer to the lowest unoccupied molecular orbital (LUMO), namely, the electron is the majority carrier.
Yet, single-molecule junctions offer the promising thermoelectric energy conversion owing to the tunability of their Seebeck coefficients [
11-
13]. There seems to be a linear dependence of the thermopower with the molecular lengthes for 1,4-benzenedithiol (8.7 μV/K), 4,4-dibenzenedithiol (12.9 μV/K), and 4,4-tribenzenedithiol (14.2 μV/K) molecules [
7]. When a C60 dimer bridges the gap between Au electrodes, it exhibits negative thermopower values with an average value of -33 μV/K which doubles the magnitude observed in the monomer C
60-based junctions [
9]. Other than increasing the molecular length, there may be other ways of tuning the thermopower, such as by strengthening the distance between the left/right electrodes or by controlling the chemical potential within a gate [
14,
15]. For example, a pure spin thermopower or/and a pure spin current can be obtained by only tuning the gate voltage in a single-molecule-magnet junction [
16].
In this manuscript, we choose the B80-based single-molecule junction as a typical model and investigate the effect of the strain as well as the chemical potential and their co-effects on its thermoelectricity. The investigation of thermoelectricity is of crucial significance in the design of the heat-to-electricity generation whose efficiencies are dependent on a dimensionless figure of merit (ZT),
Here,
κ is the total thermal conductance, which should include all possible contributions (electron (
), phonon (
) and other heat carriers (
) like photon). To account of the strain-induced modifications of thermoelectric properties of B
80-based single-molecule junctions, we carried out
, as well as
computations by systematically adjusting the closest approaching distance between electrodes. On the other hand, the chemical potential by a gate is an efficient way to adjust the electron and thermal transport in single-molecule junctions [
17-
19]. Then, in this work we hope to map out the trends by investigating the influence of the strain as well as the chemical potential on the thermoelectric transport through B
80 molecule. However, a huge number of degrees of freedom in single-molecule junctions cause a particular challeng to establish the general trend.
2 Methodology
The transmission function [] is an effective parameter to designate the coherent transport of single-molecule junctions. First of all, a set of functions [] are introduced as
with
h and
as the Planck constant and Fermi-Dirac distribution, respectively.
μ is the chemical potential regulated by a gate in a double-gate experiment [
20].
Additionally,
σ,
S, and
are expressed as [
21]
and
respectively. Obviously,
S,
σ, and
all depend on the electron transport behavior in single-molecule junctions. Actually, they are highly correlated with each other in bulk materials, thus leading to a dire challenge to enhance the thermoelectric figure of merit [
22-
24].
Lastly, ZT can be obtained by
where
is completely determined by the electron transmission function. For single-molecule junctions, it is reasonable to assume a vanishing phonon-derived thermal conductance. On one hand,
varies in a wide range of magnitude and is very small in single-molecule junctions [
25-
27]. On the other hand, the phonon transport is found to play a minor role in the thermal conductance in C
60-based single-molecule junctions [
28]. In comparison with C
60-based single-molecule junctions, B
80-based single-molecule junctions have a better electron transport behavior due to the strong coupling between B
80 molecules and metal electrodes. Therefore,
can be used to qualitatively describe the thermoelectric performance of B
80-based single-molecule junctions. In addition, electron-vibration coupling is not taken into consideration, which typically tend to reduce the thermoelectric efficiency [
29,
30].
was calculated by using the non-equilibrium Green’s function method as realized in the SMEAGOL package [
31,
32]. The norm-conserving Troullier-Martins pseudopotential is applied to describe the core electrons [
33,
34]. The exchange-correlation functional is described by the Perdew-Burke-Ernzerhof parametrization of generalised gradient approximation [
35]. The valence electrons are expanded in double zeta plus the polarized basis sets. The Brillouin zone has been sampled with
points within the Monkhorst–Pack sampling scheme [
36]. The electrostatic potentials are determined on a real-space grid with a mesh cutoff energy of 300 Ry.
3 Result and discussion
As sketched in Fig. 1(a), a B80 molecule is rigidly connected to two semi-infinite Au(100) nano-electrodes. In addition, both the planar surfaces of electrodes are taken into consideration. To start with, the total energy-distance calculations were executed and are shown in Fig. 1(b). It is found that the connection between the total energy and the distance between the left/right electrodes is parabolic when the distance is in the range from 11.0 Å to 12.2 Å. That is, in this range the deformation remains as an elastic deformation. When the distance is less than 11.0 Å, the single-molecule junction is in the plastic deformation, thus resulting in a decreased total energy. Consequently, it is assumed that the equilibrium distance between the left and right electrodes is 11.4 Å.
Experimentally, the ambient temperature (
T) of single-molecule junctions is changeable. Consequently, we investigate
σ,
S,
, and
as a function of the temperature (
T) from 50 to 1000 K, which are displayed in Fig. 2. Note that the distance between the left/right electrode is fixed to 11.4 Å. As presented in Fig. 2(a) and (c),
σ is insensitive to
T while
boosts linearly with the rising of
T. Particularly,
increases from 0.14 nW/K to 3.01 nW/K when
T ascends from 50 K to 1000 K. Actually,
σ and
correlate to each other and are constrained by a Lorenz number (
L) within the Wiedemann–Franz law (
). The obtained
L is approximately
W·Ω/K
2 and is very close to the famous value (
W·Ω/K
2) in metals or degenerate semiconductors. This is because B
80-based single-molecule junction has a well-behaved electrical conductivity which can be seen from
, as appeared in the inset of Fig. 2(a). From Fig. 2(b), the sign of
S is negative. That is, the electrons are the majority carriers. In addition, the absolute
S increases with
T increasing and reach its maximum (
μV/K) at 750 K, which is typically rather low and comparable to the measured
S in C
60-based single-molecule junctions [
9,
10]. Additionally,
in B
80-based single-molecule junctions reaches the maximum at 700 K and is about 0.035. In other words, the thermoelectric performance of B
80-based single-molecule junction is very poor, because the relationship between
σ and
is bound by the Wiedemann–Franz law [
23]. Furthermore, the Wiedemann–Franz law is proved to be invalid in a single-electron and Fano-resonant junctions, further leading to a better thermoelectric performance [
37-
40].
In B
80-based single-molecule junctions, the strain can be termed as the difference between the electrode-electrode distance and its equilibrium distance (~11.4 Å). The negative sign suggests a compressive strain and the positive value implies the tensile strain. In this work, the strain is in the range of
Å. Though there is a plastic deformation in case that the compressive strain surpasses -0.6 Å, the junctions are stable in this region. As depicted in Figs. 3(a) and (c),
σ and
are also severely dependent on the tensile/compressive strain in B
80-based single-molecule junction as in C
60-based single-molecule junction [
41]. They goes up at the same time when the strain declines from the tensile (0.8 Å) to compressive (-0.8 Å) strain, because the strain provokes the enhancement of the contact coupling between B
80 molecules and Au electrodes.
σ changes from 66.74 μA/V (0.8 Å) to 171.26 μA/V (-0.8 Å) and
increases from 0.52 nW/K (0.8 Å) to 1.12 nW/K (-0.8 Å).
In general, the effect of the strain on the thermopower is very complicated in B
80-based single-molecule junction, as shown in Fig. 3(b). However, both the experimental measurements and theoretical calculations show that
S will escalate considerably with increasing the electrode-electrode distance in a benzene-1,4-dithiolate single-molecule junction [
14,
15]. In B
80-based single-molecule junction, the absolute
S decreases when the strain increases from 0.2 Å to 0.8 Å and decreases from 0.2 Å to -0.4 Å, respectively. Unexpectedly, the sign of
S changes into “+” when the compressive strain is -0.6 Å. And then, the sign of
S changes back into “-” under -0.8 Å compressive strain. That is, B
80-based single-molecule junction can be changed from n-type to p-type and then to n-type by utilizing a compressive strain. The absolute
S has a maximum value (~-24.21 μV/K) at 0.2 Å tensile strain and its corresponding maximum
values are 0.023, as shown in Fig. 3(d). As a result, the strain can be used to modulate the thermoelectric transport properties of single-molecule junctions, but its effect is more complex.
In order to resolve the influence of the strain in detail, the transmission spectra
via the strain are shown in Fig. 4. As the strain diminishes from 0.8 to -0.4 Å, there are a few common characteristics in the variation of transmission spectra: the LUMO shifts to the chemical potential of metal electrodes and becomes more and more broadened as a result of the enhanced hybridization between B
80 molecular orbitals and the surface states of left/right electrodes [
42]. Therefore,
σ and
increase in B
80-based single-molecule junctions when the strain shifts from 0.8 Å tensile strain to -0.4 Å compressive strain. In addition,
S can be further simplified using the Sommerfeld expansion as
Obviously, S is related not only to the magnitude but also to the slope of in the vicinity of the Fermi level. Here, the maximum absolute S locates at 0.2 Å tensile strain instead of -0.4 Å compressive strain where a transmission peak locates at the chemical potential of metal electrodes. When the strain is in the range Å, S is “-” because . On the contrary, S is “+” because at -0.6 Å compressive strain.
From the above-mentioned equations,
σ,
S, and
are governed by the chemical potential (
μ), which can be controllably adjusted by a gate in the dual-gate experiments [
17,
18,
20]. Therefore, for the equilibrium B
80-based single-molecule junction we illustrate the changes of the thermoelectric properties as a function of
μ in Fig. 5. Obviously, as shown in the (a) and (c) panels,
μ induces
σ and
change along evenly together. When
μ shifts to -0.16 eV, the minimum
σ and
are 39.11 μA/V and 0.32 nW/K, respectively. Accordingly, the effect of
μ on the thermoelectric performance in B
80-based single-molecule junctions chiefly rely on the changes of
S along with
μ. As shown in Fig. 5(b), the results show that the characteristics of
S are sensitive to
μ. The most striking feature is that the B
80-based single-molecule junctions can be transformed from n-type to p-type by shifting
μ.
S is close to zero when
μ shifts to -0.16 eV. This can be noticeably comprehended from the transmission spectrum, as revealed in the inside panel of Fig. 2(a). A valley appears at -0.16 eV in the transmission function. As
μ further decreases, the sign of
S becomes “+” (p-type) because
. The maximum
S are 23.56 μV/K and -22.77 μV/K in case that
μ shifts to -0.35 eV and -0.05 eV, respectively. The corresponding maximum
[Fig. 5(d)] are 0.023 and 0.021, respectively.
Finally, we consider the combined action of the strain and the chemical potential on the thermopower and thermoelectric figure of merit, as shown in Fig. 6. From Fig. 6(a), when μ shifts to about -0.4 eV, the maximum positive S are obtained, irrespective of the strain applied to B80-based single-molecule junction. On the other hand, as the strain increases from -0.4 Å to 0.8 Å, the minimum negative S locates from -0.14 eV to 0.16 eV. In general, the maximum positive and minimum negative thermopower are 36.28 μV/K and -26.07 μV/K under 0.8 Å tensile strain when μ shifts to -0.38 eV and 0.16 eV, respectively. The corresponding are 0.054 and 0.027, as displayed in Fig. 6(b). In comparison with the equilibrium case, it can be increased a two-fold by the co-action of the strain and the chemical potential.
4 Summary
By utilizing density functional theory together with non-equilibrium Green’s function, the thermoelectric properties of B
80-based single-molecule junctions are explored by taking into consideration the effects of the temperature, the strain and the chemical potential. Mostly, the calculated thermoelectric figure of merit (
) is rather small. There are two reasons for this. First, the coupling between B
80 molecule and Au electrodes is robust and the transmission spectrum is expanded, thus leading to a smaller Seebeck coefficient. Second, the relationships between
σ and
obey the Wiedemann–Franz law very well and the Lorenz number is close to the famous value in metals or degenerate semiconductors. In response, some works recommend to improve the thermoelectric performance by violating the Wiedemann–Franz law due to Fano effect [
38-
40] or reducing the Lorenz number due to pudding mold type band structures [
43].
Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature 2019