1. Introduction
The study of molecular conduction and molecule-scale devices is and has been an active area at the boundary between chemistry and physics for at least half a century [
1], and by now, it has accumulated a substantial number of textbooks in the literature, e.g., [
2,
3,
4,
5,
6]. Over this period, different theoretical approaches have emerged. Techniques for calculation have mainly favoured Green’s Function (GF) methods, both equilibrium and non-equilibrium GFs [
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19]. Especially in chemistry, there has been a sustained development of scattering models for the ballistic conduction of electrons though molecular systems, with an emphasis on conduction through channels governed by energies and symmetry characteristics of orbitals [
20,
21,
22,
23,
24,
25,
26,
27,
28,
29,
30]. An advantage of such pictures is that they lend themselves to even simpler qualitative modelling of the effects of electron energy, relative interaction strengths, contact placement, and qualitative features of molecular electronic structures, as treated in e.g., [
31,
32,
33,
34] and as illustrated by many examples in our own work, as cited below.
One successful model of this type is the Source-and-Sink-potential (SSP) model proposed and developed by Ernzerhof and co-workers in tandem with their DFT approach to an a priori calculation of conduction [
35,
36,
37,
38,
39,
40,
41,
42,
43,
44,
45,
46,
47,
48,
49,
50,
51,
52,
53,
54]. The SSP approach replaces the doubly infinite system of molecule and leads by the finite system of a molecule dressed with a source and a sink pseudo atom equipped with complex potentials, hence in the tight-binding model replacing an
matrix problem by an
problem [
40,
55]. These ideas have their roots in earlier work [
56,
57] and have parallels in approaches such as SSM [
58] and CAP [
59].
A special feature of the SSP model in its graph theoretical incarnation [
55] is that it gives a useful qualitative account of the two-lead device, giving selection rules for conduction, rationalisation of the sensitivity of current to placing of leads; models for composite devices; and a general classification of types of conduction behaviour, such as equiconduction, omniconduction, omni-insulation, and perfect reflection [
60,
61,
62,
63,
64,
65,
66,
67,
68,
69]. Although the graph theoretical version of SSP is based on the Hückel (tight-binding) model of electronic structure, many of its most useful features persist at higher levels of theory. Conduction of the Hückel device can be analysed in terms of internal channels based on molecular orbitals of the central molecule, which are either
active or
inert [
63]. It can be shown [
70] that, as the level of theory is raised, a similar analysis applies to Hartree–Fock MO channels and, then at the second-order (or higher) Green’s Function level, to internal channels associated with Dyson orbitals. This evolution of the channel picture is consistent with the standard Meir–Wingreen [
19] and Landauer–Büttiker [
71,
72] expressions for elastic conduction of electrons through molecules.
Green’s Function approaches have been applied to various multi-lead configurations at various levels of theory. For instance, Reference [
59] quotes a dozen examples, and public domain codes such as TRANSIESTA [
73] allow for calculations with many leads. Explicit high-level calculations of transmission for individual systems are undoubtedly valuable, but there is still a role for simple models that can identify broad trends for families of systems, can give a framework for interpretation, and can identify cases for more detailed analysis by a priori/ab initio methods.
In line with this philosophy, our aim in the present paper is to generalise the two-lead SSP formalism to deal with multi-lead devices, retaining its many advantages. To be concrete, we call this version of SSP the SMSP (source-and-multiple-sink-potential) model. The paradigm system for this model has a single-source lead connected to one atom of the central molecule, and several drain leads attached to other specific atoms. Ernzerhof already considered an extension of SSP with a three-lead device that has a gating arrangement that includes two sources and one sink [
51]. We defer treatment of the many possible multi-source variations of SMSP to future work.
As with SSP, we are able to obtain an analytical solution at the Hückel/Tight-Binding level for non-trivial molecular systems in terms of a set of characteristic polynomials. General expressions for the currents in all the leads can be developed and used to identify global selection rules. They suggest new classes of conductors, by analogy with developments in the theory for the original two-lead device. The model shares with SSP the same opportunities for seamless transition from Hückel to self-consistent and correlated treatments.
The structure of the paper is as follows:
Section 2 describes the basic SSP model;
Section 3 introduces the SMSP model equations, their formal solution and simplification for devices with symmetric (i.e., chemically identical) leads, and the analysis of molecular orbital conduction channels and lead-lead interference effects;
Section 4 treats two classes of ‘complete’ devices, in which either all vertices of the molecular graph (CSD) or all vertices of a given partite set (CBSD) are connected to leads and which turn out to have minimal interference effects;
Section 5 shows how to break down conduction into contributions from internal molecular channels;
Section 6 reports illustrative results,
Section 7 outlines the connections of the present model with the Meir–Wingreen approach; and
Section 8 reports our main conclusions.
2. Two-Lead Devices
In the standard two-lead SSP model, the device consists of a central molecule that is attached by single-atom contacts to the leads (
Figure 1). In the steady state, the source lead supports an incoming electron wave that undergoes partial reflection at the interface with the molecule, and the transmitted part of the wave scatters through the molecule and exits into the sink lead. The essential observation behind the model [
40] is that the two semi-infinite leads may be replaced by source and sink pseudo-atoms equipped with appropriate complex potentials, which play the roles of delivering the transmitted fraction of the electron to the molecule and removing it downstream, respectively. In the purely graph theoretical version of the model, this step allows replacement of a doubly infinite system by one in which the
Hermitian (in fact, real symmetric) molecular Hamiltonian matrix is expanded to an
non-Hermitian matrix. It is convenient to recast this generalised eigenvalue problem in the form of an inhomogeneous linear equation [
63]:
where
is the SSP matrix, the solution vector
gives the device wavefunction specified by complex entries on the
vertices, and
is the inhomogeneity term that embodies the boundary conditions. The solution to this equation for each possible value of the energy of the incoming electron yields a compact expression for the zero-voltage transmission as a function of energy,
, and as a bonus, gives a useful analysis of current in terms of bond or molecular-orbital (MO) contributions.
The transmission function can be cast in terms of the characteristic polynomials of four graphs,
G,
,
, and
, where
and
are the centres within the graph that are in direct contact with the leads. These polynomials are denoted
s,
t,
u, and
v in [
55]. Details of the algebra are given in the reference, but the result is that the device transmission has the form [
63]:
where
is a band-pass function ensuring that the electron energy is within the conduction bands of the leads, and the denominator is
where the polynomials
s,
t,
u, and
v are to be defined below in (
7), The incoming and outgoing wavevectors
and
satisfy the Hückel dispersion relations for the leads:
and we take the symmetric case where
, and
, and refer the zero of all energies to
. The polynomial
satisfies [
74]:
where the four structural polynomials (SPs),
s,
t,
u, and
v, which are each the characteristic polynomial of a graph, are defined with respect to the characteristic matrix
, where
is the adjacency matrix of the molecular graph:
and
j can be calculated directly [
55,
74] as
In these equations, superscripts indicate the sets of rows and columns to be deleted from the characteristic matrix.
The transmission
is therefore a ratio of polynomial functions of
E, taking values between zero and one and depending implicitly on the molecular graph and the placement of connections. Questions such as transparency, opacity, selection rules, and linkage of transmission to MO channels can all be explored either analytically or numerically using this function [
55,
60,
61,
62,
63,
64,
65,
66,
67,
68,
69].
In going to the more general device with multiple leads, the conceptual framework remains intact (
Section 3), but the notation needs to be overhauled. The first point of difference concerns the naming of leads. We replace the
and
nomenclature for leads with numerical labelling. For the construction of derivations and proofs, it is convenient to adopt the convention that the source is lead 1, and the sink leads have labels
to
l (
). Centres in the molecule are then numbered 1 to
n such that, for the first
l centres,
in the molecule is connected to lead
, and the numbers
to
n are then used for any centres that are not connected to leads. Formulas for transmission use this convention so that
describes the transmission from the source to the molecule and
is the transmission from the molecule to sink lead,
. Contiguous numbering of this type is advantageous for proofs that use the block structure of the characteristic matrix of bipartite molecular graphs, for example. (In this paper, we consider only devices where all leads are connected to distinct vertices of the molecular graphs. Degenerate ‘ipso’ [
60] devices would need a modified numbering convention.)
A second adjustment is to the nomenclature for SPs. For larger numbers of leads, the labelling of every SP with a different letter becomes unwieldy, so we adopt a single symbol that can be extended to arbitrarily large cases. It is convenient to define
where the SPs of increasing order are now defined as cofactors (cf. Equation (
8)) of the characteristic matrix and where the indices indicate the rows (left of the comma) and columns (right of the comma) to be struck out. We also chose to scale each SP by multiplying by the connection parameters between each deleted vertex and its connecting lead, a feature that makes the ensuing derivations simpler. Diagonal SPs, those for which the two deletion lists are equal, are now scaled versions of the four structural polynomials from the earlier two-lead papers, e.g.,
We can also define some related quantities by dividing by
s, i.e.,
These ‘hatted’ SPs are important in two respects. First, they emerge naturally in the solution of the SMSP equations, as we shall see in the next section. Second, they are intimately connected with the molecular GF, which in the Hückel tight-binding formalism is the inverse
An extension of a theorem by Jacobi (cf. Equation (
12), p. 774, [
75]) establishes that all our hatted polynomials are determinants:
It is also convenient to use an extra piece of notation for the
traces of principal minors of SPs. We define the traced quantities,
where the sums are over the set of molecular vertices attached to leads. The transmission for the two-lead device, expressed in this revised formalism, is
The reader may object that this expression could diverge when the value of
E is equal to a molecule eigenvalue,
. The divergences in (
15) are easily removed. The zero-rank polynomial,
s, is expressed in terms of molecular eigenvalues through
and if we multiply numerator and denominator of (
15) by
, we obtain
which with the equivalences
and
is identical with the expression for
in [
55].
Expression (
17) has no divergences at molecular eigenvalues and has been used in our previous work [
60,
63,
70] to derive selection rules for transmission. The numerical computation of SSP transmission does not use expressions such as (
15) directly, relying instead upon a standard linear numerical solution of the SSP matrix equations that does not suffer from numerical instabilities. The purpose of equations such as (
15) is for analysis to reveal the physics involved in the transmission process. Our aim in this paper is to generalise these expressions for multi-lead devices.
3. Multi-Lead Devices
We consider a device comprising a molecule with
n atoms and a set of
l leads (
Figure 2). The leads are simple ‘unstructured’ infinite atomic chains described in a Hückel tight-binding model in which a given lead,
, has edge weights
and has a simple connection to atom
(using an appropriate molecular numbering scheme) in the molecule with edge weight
. Without loss of generality, the first lead is the single source lead, and the remaining
leads are sinks. Each lead has a dispersion relation
that relates the energy,
E, of a stream of electrons in the lead to the parameter
and the associated wavevector
(
in an infinite wire. The description of the
molecule may, in principle, go beyond the Hückel model in the manner we have described previously [
70], but, for simplicity, here, we derive the Hückel version of the theory.
3.1. The SMSP Equations
The Source-Multiple-Sink potential (SMSP) equations are modelled on our previous work with SSP [
63], so that we can write them in the form (
1) as
where the
l-dimensional (diagonal) lead block of the S(M)SP matrix,
, is
Each nonzero element in is the matrix element for the terminal atom of the inverse Green’s function of the semi-infinite linear chain, .
The molecule block is the characteristic matrix constructed from the molecular adjacency matrix,
,
The diagonal matrix,
, is the matrix of
n molecular-orbital eigenvalues, and
is the associated
n-dimensional matrix of molecular-orbital eigenvectors. The unitarity of the eigenvectors implies that
The lead-molecule block,
, is an
-dimensional matrix with elements
On the right-hand side of Equation (
19), the inhomogeneity vector,
, has a single non-zero element
which corresponds to the fact that lead 1 is the single source lead. The vector
ensures that the system of equations satisfies the correct boundary conditions for a semi-infinite source lead with current travelling in the forward direction with pseudo-momentum,
. The factor
is the normalising factor (cf. [
63]), ensuring that the forward current is equal to unity.
The SMSP solution vector
in Equation (
19) is the (complex) device wavefunction. Consequently, we can write the transmission from the molecule to lead
using the standard formula for the electron current [
2] as
This formula has the correct sign for transmission into sink leads (p ≠ 1), but the sign must of course be reversed to obtain the transmission from the source lead (p = 1).
3.2. The Formal Solution to the SMSP Equations
The SMSP equations in Equation (
19) can be written out in terms of the separate lead and molecule blocks as
We can solve for
in terms of
using Equation (
21):
We can now substitute our formal solution for the molecular coefficients into the lead equation to give
where the reduced lead matrix,
, is
Equation (
30) can be solved by inversion to give
where the cofactor matrix and the determinant are defined as
and from Equation (
29)
Substitution of Equations (
32) and (
34) into (
26) gives
where we moved real quantities to the left of the imaginary part operator. Summing Equation (
35) over all leads,
, and using the molecular GF for real energies gives
since the terms under the summation on the right-hand side are all real, thus verifying Kirchhoff’s law for the conservation of current in the device.
From the definition of the reduced lead matrix in Equation (
31),
We can now use the Laplace expansion of the determinant
to show that
For a
sink lead (
), we deduce that
It follows from Kirchhoff’s conservation rule that
Equations (
35) and (
40) give the general SMSP expressions for the currents that flow through the single-source molecular device. These expressions can be used directly for the calculation of ballistic currents in general SMSP devices. However, some far-reaching simplifications are available when we restrict consideration to the usual case where a device consists of a molecule connected to chemically identical leads, i.e., the
symmetric device.
3.3. Symmetric Multi-Lead Devices
The leads in symmetric devices have identical lead parameters,
, and identical connection parameters,
. This enables a relatively easy simplification of the transmission formula (
39) into a cluster expansion, involving pairs, triples, etc. of leads. All leads satisfy a single dispersion relation
and furthermore, the reduced lead matrix of Equation (
31) is now
where the bar over the matrix
indicates that we require rows of the matrix only for those atoms connected to leads. The rectangular matrix
is obviously not unitary and indeed gives different products
but the first of these and the diagonal nature of the first term in the
l-dimensional matrix
, allow us to write
where the
diagonal n-dimensional matrix,
, which involves all eigenstates, is
Equation (
44) expresses the reduced lead matrix,
, as a triple product involving rectangular matrices. The Cauchy–Binet theorem [
76,
77] can be used to expand the determinant of a product of two rectangular matrices. Hence, we obtain a further reduction in the quantities in Equation (
33) as a sum over a set of products of
l-dimensional principal minors
where
is an ordered subset of
, with
, and
. By deriving this equation, we are able to separate the roles of eigenvectors,
, and the state energies contributing to
. Further manipulation of Equation (
46) uses algebra familiar in the reduction of matrix elements over Slater determinants in quantum chemistry. We can replace the restricted summation over the ordered sets,
, by unrestricted summations from 1 to
n over the state indices (i.e., Hückel MO or GF pole indices),
. Equalities between the
indices are allowed because they give rise to zero contributions in Equation (
46), since they produce determinants,
, with identical columns. We can also use the antisymmetry of the determinants to replace the first determinant by a simple product. Hence, we can conclude that
The next step is to expand the product of the
using Equation (
45) to generate a cluster expansion in powers of
and
. Whilst it is easy to pursue this generally, it is probably more enlightening to consider the example of a three-lead case explicitly.
The Symmetric Three-Lead Device
In the case of the three-lead device, the denominator (
47) is
The central triple product can be expanded exactly in the form
where
and
are in the set
, and we have defined energy denominators
Our aim is to obtain an expression in terms of generalised SPs. We proceed by substitution of Equation (
49) into Equation (
48) and then by examining the resultant coefficient of each power of
. The first term in (
49) generates a coefficient
where
is a permutation of parity
belonging to the permutational group
of the lead labels 1, 2 and 3, and we used the orthonormality of the
matrix rows.
Using the same methodology, the
term in Equation (
49) gives three contributions to the coefficient, the first of which is
Likewise, the other two denominators give contributions
and
. Reduction in the remaining terms proceeds in a similar manner, but the increased number of denominators is analogous to going from one-electron to two- or three-electron matrix elements over Slater determinants, and it is necessary to include more of the permutational symmetry of the included determinant. Collecting terms, we reach an explicit formula for
or simply
where we have used the Jacobi relation from Equation (
13) and the generalised SP notation in Equations (
11) and (
14). We derived this theory in terms of the Hückel tight-binding approximation, but it also holds for correlated molecular GFs, where the MO coefficients,
, in the definitions of SPs are replaced by the GF Dyson orbital coefficients [
70]. Furthermore, the cluster expansion in Equation (
53) holds for devices with any number of leads. The alternating sign expansion contains terms with determinants of GF matrix elements with increasing dimension up to
l, the total number of leads in the device.
Now that we have the denominator for the sink lead transmission, we also need to derive an expression for the cofactor matrix element in Equation (
39). As an example, we take the
matrix element for the three-lead case. We have the same structure as before, but lead 2 matrix elements are omitted from the left-hand part of the expression and lead 1 matrix elements are omitted from the right so that the determinants are of order two:
We write out the determinant in full. The simplification of this expression is carried out exactly as before. The first term,
, in the expansion of the
factors, however, gives no contribution because the off-diagonal nature of the cofactor creates matches only between orthogonal
matrix elements. The final result for the cofactor is
The difference between the cluster expansion of the determinant in Equation (
53) and the cofactor in Equation (
55) is significant. The determinant expression has traces over principal minors. These satisfy the interlacing theorems we previously utilised in our work on two-lead devices to deduce selection rules [
60].
It is straightforward to extend (
53) and (56) to more leads. For example, the equivalent formulas for four-lead SMSP devices are as follows:
3.4. The Wide-Band Limit
The wide-band limit (WBL) is a commonly used approximation for calculation of electron transport [
41]. It is assumed that the wave-vector of the ballistic electron can be taken to be
, and all other quantities in the expression for the transmission tend to their values for
. In the SSP model of the symmetric two-lead device, the wide-band limiting transmission,
, for non-ipso devices reduces to a formula in terms of the tail coefficients of the characteristic polynomials of the four subgraphs induced by deletion of 0, 1, and 2 distinct connection vertices of the molecular graph
G [
55]. This expression has been used to develop selection rules for Fermi-level conduction of two-lead devices [
60,
61].
A similar approach can be used for many-lead devices. As an example, the WBL version of the sink transmission formula for lead 2 in the symmetric three-lead device is
This can be expressed entirely in terms of the characteristic polynomials of the seven subgraphs induced by deletion of 0, 1, 2, and 3 distinct connection vertices as follows:
where
, and we understand the lower-case symbols
as the Fermi limits of the characteristic polynomials of the graphs
G,
,
,
, and the upper-case symbols
as their traces. This expression reduces to the corresponding two-lead equation (Equation (21) in [
55]) when all SPs involving lead 3 are set to zero. It is straightforward to extend Equation (
59) to arbitrary numbers of leads. As the seven polynomials in the three-lead case are linked by an identity for the principal minors of the
hyperdeterminant (see, e.g., Equation (2) in [
78]),
can be calculated from all seven or just six of them. In either case, the WBL approximation can be applied to this instance of a many-lead device to give the transmission from a knowledge of the graph itself and its connections. Only the tail coefficients of the characteristic polynomials are required; this will have particular significance for future development of multi-case selection rules based on the numbers of zero roots of the various polynomials.
3.5. Analysis of Interference in Multi-Lead Transmission
In order to appreciate the nature of the transmission for multi-lead devices, we need to examine the numerator and denominator quantities
and
. Again taking the three-lead case as the example, we have
The first term in the expansion involves the same rank 1 off-diagonal SP that appears in the two-lead device constructed from leads 1 and
, albeit with a different
denominator. We shall refer to this as the
direct term for the
constituent device . The remaining two terms represent the effect of the third lead (i.e.,
) on the conduction through the second lead (
). We shall call these
mixed interference and
pure interference terms, respectively. Each term has a different dependence on the lead parameter,
, and the lead-molecule parameter,
(cf. Equation (
9)). The direct term varies as
, whilst mixed and pure interference components vary as
and
, respectively. The different scaling implies that the relative importance of direct and interference effects may change markedly for different choices of the relative sizes of these parameters.
The denominator term in the transmission is expressed using our trace notation as
Equation (
61) is derived for the symmetric three-lead device with
, and so the two-lead device cannot be recovered simply by taking
, but in the limit that
, Equation (
61) tends to a product of a factor
multiplied by the denominator in Equation (
15). The transmission formula for sink lead
in multi-lead symmetric devices has the pre-factor
, which gives the quadratic rise and fall in transmission as a function of energy near the lower and upper band edges, as in the two-lead device.
The analysis of interference terms can be extended to the case where the device has more leads. In the case of
l leads, the leading term in the numerator is of the form
It has the same SP as the cognate constituent two-lead device, with a dependence on , and this can be expected to be the most important term in the transmission spectrum. The other terms in represent sink-lead interference terms. They form a cluster expansion with the inclusion of higher rank SPs up to the number of leads, l, and traces over all sink leads other than . An analysis in terms of a hierarchy of constituent two-lead, three-lead, etc. devices could be envisaged.
The equations for symmetric devices, exemplified by (
60) and (
61) can be used to derive general features of transmission spectra, selection rules for Fermi transmission, and systematic trends in interference effects. Even further simplifications are possible with additional assumptions about the choice of leads, in particular when the number of leads is large.
6. Results
The SMSP formalism described in
Section 3 has been implemented in our suite of routines that use the Maple 2019 package [
81] to provide analytical and numerical calculations of the zero-voltage transmission curves
for multi-lead systems.
The Hückel tight-binding calculations that are used for the central molecule use a single
-orbital basis function on each atom with a hopping parameter,
, to represent the interaction between
-bonded carbon centres. Hartree–Fock (HF), and second-order Green’s function (GF2) calculations require two-electron
-interactions. We used a Hamiltonian with a Hubbard single-centre interaction parameter,
, for all centres and a single two-centre interaction parameter,
, between atoms that are
-bonded. The methodology is identical to that in our earlier work [
70] for the two-lead SSP model, where it is described in more detail. No two-electron interactions were used to describe the source and sink atoms that represent the leads in the SMSP method. For most calculations, we used
for the lead hopping parameter and
for the lead molecule connection parameter. All energies are shown in units of the negative quantity
.
The labelling schemes for the molecules under study are shown in
Figure 3.
Figure 4,
Figure 5,
Figure 6 and
Figure 7 show transmission curves plotted against electron energy. In each case, the black outer envelope is the total transmission from the source lead into the molecule, and coloured curves show (symmetry distinct) transmissions from the molecule into sink leads. A colour code is used to distinguish the leads in clockwise order from the source. The plots also indicate poles of the molecular Green’s function, using an encoding of blue for attachment and red for ionisation poles.
6.1. Devices Based on the Benzene Ring
As a first example, we take the six-membered ring of carbon atoms, which can be connected as a two-lead device in three ways (‘ortho’, ‘meta’, ‘para’) and as a three-lead device in six ways, leading to the nine transmission curves shown in
Figure 4. All curves are symmetric about the Fermi energy in the Hückel model, as the six-ring is a bipartite graph. Transmission curves for the three two-lead devices are qualitatively different. The ortho and para devices have non-zero transmission at the Fermi energy, with different patterns in the wings, as predicted by the analytical expression for
of a cycle [
55] and selection rules [
61,
63] for conduction at degenerate and non-degenerate eigenvalues (here,
). In contrast, the meta device is insulating at the Fermi level.
Plots for the various three-lead devices show a variety of patterns, but it is striking that the overall shapes of the curves can be interpreted in terms of constituent two-lead components: to a first approximation, the curve
for the three-lead device with source 1, and sinks
a and
b has lead contributions that follow the patterns for devices
and
in the region of the Fermi level but are damped in the higher and lower energy wings. The implication is that the interference terms (see Equation (
60)) have smooth, predictable effects on
for the composite device. This pattern is especially clear for the symmetric
device (
Figure 4i), where the source is in a meta relations to both sinks, and the Fermi-level transmission is zero. In fact, this device is a CBSD: interference terms between leads do not appear in the numerators of the lead currents, and the device is interference-free. In other devices such as
and
, a vanishing Fermi-level contribution from the meta constituent device is masked by the contribution from the other constituent. The counting of ortho, meta and para constituent devices works well as a rough guide to
, though we may expect this simplicity to be diluted in larger systems and devices with more leads.
The final row of
Figure 4 shows the transmission curves calculated at the GF2 level for
,
constituent devices, flanking the plot for the composite three-lead device
. The plots show strong resemblances to those calculated at the Hückel level, and interpretation of the lead contributions in terms of constituent devices survives intact in the correlated treatment.
6.2. Devices Based on Anthracene
This section presents results on devices based on anthracene, which, similar to benzene, has a bipartite molecular graph. The devices considered here all have three leads, with the source on atom 4 in the numbering of
Figure 3. The devices for which transmission curves are shown in
Figure 5a–c are (4,8,14), (4,9,13) and (4,8,11), which span the types (
), (
), and (
), where the colours of the circles represent the partite sets (
Figure 3b). All three necessarily have single-vertex deleted and triple-vertex deleted graphs that are non-Kekulean and have at least one non-Kekulean constituent two-lead device. Devices (4,8,14) and (4,9,13) are C
symmetric, and their total tramsission is simply twice the contribution of each sink lead. The second row of the Figure shows the three two-lead constituent devices (4,8), (4,9), and (4,11). In each case, there is a clear qualitative correspondence between the three-lead device and its constituent devices, i.e., (a) with (d); (b) with (e); and in the mirror-symmetric (c), the distinct lead contributions strongly resemble (f) and (d), generally as expected. Exceptions to this agreement are in the features where the transmission falls sharply to zero at eigenvalues
and
in the two-lead (4,8)- and (4,9) constituent devices; these are suppressed in the mirror symmetric three-wire devices as an effect of specific cancellations with the denominator.
The final row of panels
Figure 5g–i refers to interference terms in lead currents in the three-lead devices. Explicit calculation of interference contributions to the numerators of the expressions for the lead currents shows that the mixed and pure interference contributions are identically zero for the mirror-symmetric devices (4,8,14) and (4,9,13) but not for the non-symmetric one (4,8,11). As shown in
Figure 5i, interference terms of both mixed and pure types appear. It is risky to generalise too much from one example, but it is at least interesting to see that these terms exert little effect in the central region of the spectrum (where they are both small) but are more influential in the wings (where they are more intense and have regions of opposite sign).
6.3. A Non-Bipartite Case: Devices Based on Pentalene
Figure 6 shows a cascade of calculated curves
for non-bipartite (i.e., non-alternant) pentalene. We take one example of a four-lead device and compare it with all of the constituent two- and three-lead devices. The chosen four-lead case has connections at vertices 1,3,5,7 (
Figure 6a), which are the four core vertices for the non-degenerate non-bonding LUMO of pentalene (i.e., they carry the non-zero charge/spin densities for this orbital).
The plots are no longer symmetric, as the spectrum of a non-bipartite graph is not paired, but there are still discernible qualitative features of the transmission curve that are associated with eigenvalues: broad or compound maxima at most eigenvalues, and a sharp spike in the transmission at the HOMO eigenvalue ( in Hückel theory). The plot also shows a broad dip in transmission at antibonding energies above the LUMO. The origins of the lead contributions to devices in the middle row are plausibly traced from the appropriate two-lead devices ((e) and (f) to (b); (e) and (g) to (d); (f) and (g) to (d)), and the constituent device curves for (e), (f), and (g) appear in the plot for the full device (a), with some blurring of features by interference effects. The family relationships of devices with progressively increasing numbers of leads are yet again helpful as a rough guide to the features of the transmission spectrum of a complex device.
6.4. Complete Devices of Types CSD and CBSD
Figure 7 shows calculated curves
for molecular devices of CSD and CBSD types based on cycles and linear polyenes. The first three panels show the transmission curves for CSD based on 6-, 8-, and 10-membered rings (
Figure 7a–c). They show qualitatively similar behaviour, with broad maxima, shoulders, and minima associated with the eigenvalues of the respective graphs. The pattern of a central minimum for
-rings and a central maximum for
-rings persists to larger ring sizes. Traces of the constituent devices are apparent in the lead contributions, e.g., compare ortho, meta, and para benzene devices (
Figure 4a–c) with the CSD in
Figure 7a. As noted in
Section 4.1, transmission is non-zero across the whole band window. However, the summed current can include some very low sink-lead contributions. In the case of the 10-ring, for example, the lead transmissions at the Fermi level for leads 2, 3, 4, 5, and 6 are 0.1581, 0.0445, 0.0853, 0.0042, and 0.0674, respectively.
The next three panels (
Figure 7d–f) show the transmission curves for CBSD based on 8-, 10-, and 12-membered rings (
Figure 7a–c). Again, a pattern of alternating behaviour at the Fermi level is apparent and is confirmed by extension of the calculations to larger rings: CBSD based on
-rings show a central maximum in total transmission, consistent with superposition of sink-lead curves with central maxima, whereas
-rings have an insulating minimum at the Fermi level, indicating vanishing of all lead contributions. As noted in
Section 4.3, this is part of a general pattern based on the number of non-bonding orbitals (NBOs), or in other words, the nullity of the graph. In the
-ring, the nullity is 2 and all centres carry non-zero density arising from any occupation of the non-bonding shell (i.e., a
-ring is a core graph [
82]). Cancellation of the
factor between numerator and denominator in the expression for lead currents therefore yields conduction at the Fermi level for these rings. The
-rings have no such cancellation. This convex/concave pattern for antiaromatic/aromatic rings is one example of a nullity-based selection rule for a family of devices. There will be many more to be found.
Finally, the third row (
Figure 7g–i) deals with complete bipartite devices based on the heptratriene linear polyene. As shown in
Figure 3, the molecular graph has two partite sets of vertices. The graph has nullity one and the unique Hückel NBO has non-zero entries only on the larger (black) partite set. This leads to a qualitative distinction between devices of CBSD type based on the two sets. For the (1,3,5,7) device, there is non-zero density in the NBO at all lead positions, and cancellation of the
factor in the numerator leads to Fermi conduction. On the other hand, the two devices based on the smaller (white) partite set show insulation at the Fermi level, even though the underlying graph is still of course singular, as all lead positions are at nodes in the NBO.
A comparison of
Figure 7h–i shows that this central feature of the transmission spectrum is independent of the choice of source lead position within the same partite set. There is, however, a significant difference between the (2,4,6)- and (4,2,6)-devices. The molecule has a mirror plane through the central vertex, so that MOs that are antisymmetric with respect to the mirror have
. The (2,4,6)-device in panel (h) displays weak inertness for lead 2 attached to vertex 4 (in red) but not for lead 3 attached to vertex 6 (the green curve). Hence, no maximum is visible in the red curve near the second eigenvalue from the right or the second from the left. The (4,2,6)-device, in panel (i), exhibits strong inertness for these same eigenvalues because the source lead is attached to vertex 4. These observations are again indicative of the wealth of underlying selection rules.
7. Connections with the Meir–Wingreen Formula
The Meir–Wingreen (MW) formula gives an exact expression for the current through a molecular device described using a correlated molecule but with leads and lead–molecule interactions treated in the one-electron approximation. In previous work [
70], we showed that the SSP method is consistent with the MW formula, in the elastic-scattering limit when the electron–phonon interactions are neglected (i.e., in the Born–Oppenheimer approximation). SSP produces identical formulae for the transmission, including the case in which the description of the central molecule includes electron correlation.
The transmission in the SMSP approach is determined by boundary conditions on a device wave function with conceptual advantages in the chemical interpretability of the internal molecular channels for conduction. These channels are determined by the attached and ionised states appearing in the Lehmann representation of the equilibrium molecular GF. The properties of these channels are defined in terms of the characteristics of the DOs associated with each of the GF poles. The sets of orbitals and poles define spectral expansions of structural polynomials used to formulate the expressions we derived in this paper. The relevant formal properties of the polynomials are retained regardless of the level of theory, from Hückel, through HF, to sophisticated GF formalisms. The molecular part of the problem may therefore be treated with empirical, semi-empirical, or ab initio methods, with or without the inclusion of correlation, thus giving a seamless single formalism to describe conduction through molecular devices.
At any of these levels of approximation, the inverse
of the SSP/SMSP device matrix is the device GF (Equation (
85) of [
70]). The poles of this matrix are complex-valued and represent device resonances corresponding to the presence of an extra electron (or hole) in the sea of (possibly correlated) electrons. The imaginary parts of these poles are inversely related to the lifetimes of the resonances. These lifetimes represent transit times for the particles or holes crossing the device, times that are particularly short for Koopmans poles but long for shake poles. In the SMSP context, it can be seen that the poles of the device GF can be obtained by solving the equation
. Each complex conjugate pair of poles is then associated with a Dyson orbital and (real) pole energy of the molecular GF. Hence, the SMSP model gives a direct route to calculation and interpretation of all the properties that can be accessed through GF methods.
8. Conclusions
Working equations have been derived for an extended version of the SSP model for calculation of ballistic currents flowing through molecular conductors under a potential difference. They constitute the new Source-and-Multiple-Sink Potential model. The new model describes systems with an incoming electron travelling along a single source lead that drains out through the molecule into multiple sinks. As with its SSP predecessor, SMSP lends itself to a graph-theoretical approach, which leads to a formulation that is compatible with, but not confined to the Hückel/tight-binding description of electronic structure of
systems. SMSP therefore retains the advantages of SSP. It yields closed-form expressions for device transmission as a function of electron energy: for general lead parameters (
Section 3.2); for chemically symmetric leads (
Section 3.3); and for leads with coverage of all
-centres (
Section 4.1), or of all starred or unstarred centres of an alternate
-system (
Section 4.3).
As its SSP predecessor did, the SMSP model gives predictions for qualitative selection rules and generic patterns of conduction for families of molecular systems, such as the opposite behaviour of aromatic and anti-aromatic rings at the Fermi level (
Section 6.4). In addition to a calculus for the prediction of total transmission with energy, the model allows for an analysis of molecular current at various levels of detail: by bond current, by orbital-based molecular channels, and by destination lead (
Section 5). Patterns of current in multi-lead devices can be interpreted as built up from notional constituent two-lead devices, where SMSP gives a natural breakdown into direct contributions, interference terms (again governed by selection rules), and mean-field corrections (
Section 3.5). The wide-band limit yields particularly simple expressions for transmission involving only the characteristic polynomials of the set of vertex-deleted graphs (
Section 3.4). In the hypothetical limit of full or half coverage of a molecule by leads, the specific interference effects are predicted to vanish (
Section 4.1 and
Section 4.3). The symmetry of the device can also cause the vanishing of interference effects (
Section 5.3) when the connection point for the sink leads are equivalent.
SMSP, similar to SSP, is a parameterised model, with two parameters
, and
(defined in units of
, the resonance integral for the molecule). Adjusting the parameters affects the detailed appearance of the transmission spectrum and the relative importance of interference effects. The effective lead parameter
is likely to be larger in magnitude than
, even in all-carbon devices since the co-ordination number is higher with structured leads. The molecular energy levels, therefore, lie inside the band as assumed in our sample calculations. The role of
is to describe the perturbation of the molecule by the leads, the effect of which is to move poles of the molecular GF off the real axis, and hence, lower
values give sharper transmission peaks. Interference effects become more important as the ratio
increases, as shown in (
Section 3.5). Perhaps the main advantage of models such as SMSP is that many interpretative features persist when more realistic models of electronic structure are used. In particular, the analysis by molecular-orbital-based channels survives the introduction of electron–electron interaction into the theory.
Thus far, we are not aware of experimental work on molecular systems with dense attachment of leads to single molecules. Mechanically controlled break-junction (MCB) techniques [
83] provide IV characteristics that give insight into averaged conduction for two-lead molecular devices rather than their detailed transmission spectra. Even so, an impressive range of data on the factors influencing molecular conductivity can be inferred from painstaking experimentation based on careful molecular design. For example, encouraging agreement with experimentally observed regiospecificity of conduction [
84] can already be achieved using non-empirical and graph-theoretical SSP calculations [
43,
60]. The long-term aim of the present work is to carry forward this modelling of trends to multi-lead devices, in the hope that it may be of use in the design of future experiments on these intriguing systems.