REVIEW 4 major objections 4 minor 10 references
Electronic and Thermoelectric Properties of Molecular Junctions Incorporating Organometallic Complexes: Implications for Thermoelectric Energy Conversion
T0 review · 4 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read The paper claims that iron-cyclam molecular junctions reach an experimental thermoelectric figure of merit ZT up to 0.4 at room temperature, among the highest reported for molecular systems.
desk verdict Good multi-technique thermoelectric data on organometallic junctions, but the ZT=0.4 headline is computed with swapped thermal conductances and needs correction. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing objects are two trans metal-bis(arylacetylide) wires of nearly equal length and identical thiol linkers: the Ru-dppe complex and the Fe-cyclam complex, the latter present as both Fe(II) and Fe(III). The experimental machinery combines mechanically controllable break junctions for single-molecule conductance, conductive atomic force microscopy for monolayer conductance, null-point scanning thermal microscopy for heat conductance, and thermovoltage measurements for the Seebeck coefficient, with density-functional-theory-plus-Green's-function transmission calculations. The quantity that carries the argument is the slope of the transmission function $T(E)$ at the Fermi energy (the energy-dependent electron tunneling probability): a steeper slope gives a larger Seebeck coefficient, and the Fe-cyclam transmissions show steep, non-Lorentzian features, including destructive quantum interference (cancellation of transmission paths) and, for the charged state, spin-dependent channels, close to the Fermi level. The figure of merit is assembled as $ZT = S^{2}GT/\kappa$, with the per-molecule thermal conductance $\kappa$ taken from the SThM monolayer value divided by the estimated number of contacted molecules.
What would settle it
Measure the thermal conductance and thermovoltage of a single Fe-cyclam junction, or a well-defined few-molecule junction, directly rather than dividing a monolayer value by a molecule count; if the per-molecule thermal conductance comes out above about 100 pW/K, the room-temperature $ZT$ drops below 0.1.
Extended reading notes
Core claim
The central claim is that Fe-cyclam junctions, contacted through thiol linkers on gold, combine a high Seebeck coefficient with a low heat conductance, giving an experimental $ZT$ of up to 0.4 at room temperature for junctions in the high-conductance tail of the distribution ($G_{\max}\approx 4.2\times 10^{-2}G_0$). The paper states this is among the highest room-temperature ZT values reported for molecular systems. The same molecule shows two oxidation states, [Fe-cyclam]^0 and [Fe-cyclam]^+, which the authors identify from two level alignments in fits with the single-level model (a Lorentzian one-orbital description of junction transmission) and from X-ray photoemission; the two states give different conductance and thermopower fingerprints. The Ru-dppe analogue, with the same length and linkers, shows $S\approx 14$ to 28 μV/K and a much smaller ZT, so the Fe center itself, not just the wire length or the anchor group, is the decisive ingredient.
Load-bearing premise
The room-temperature $ZT$ of 0.4 depends on the assumption that the per-molecule thermal conductance of a Fe-cyclam junction is the SThM monolayer value (about 20 nW/K) divided by the roughly 1000 molecules assumed to sit under the thermal probe; if the real count or the additivity of heat flow is wrong by a factor of two, $ZT$ changes by that same factor.
Editorial extensions
If this is right
- Fe-cyclam junctions become a benchmark room-temperature molecular thermoelectric, with $ZT\approx 0.4$ competitive with the best published molecular systems.
- The Fe(II)/Fe(III) redox pair provides two electrically distinct junction states, so the same molecule can switch between higher-conductance, lower-Seebeck and lower-conductance, higher-Seebeck operation.
- The low per-molecule thermal conductance of about 20 pW/K means heat transport, not electron transport, is what the iron center suppresses.
- Room-temperature operation without external gating or magnetic field is achievable, a practical requirement for energy-conversion devices.
- The computed $ZT$ of about 0.8 for Fe-cyclam, against the estimated experimental 0.4, indicates that cleaner junctions or better contacts could push the experimental value higher.
Reading between the lines
- Beyond the paper, if the 15-50 pW/K per-molecule heat conductance range is generic for heavy-atom organometallic wires, then replacing light organic backbones with metal centers is a design rule for suppressing phonon conduction in molecular thermoelectrics.
- An extension the paper leaves implicit is that, because the two iron oxidation states sit at different level alignments, an electrochemical gate could tune a single Fe-cyclam junction between high-power and high-efficiency operating modes in situ.
- A direct test would be a combined single-molecule thermovoltage and heat-conductance measurement on the same Fe-cyclam junction; if the per-molecule $\kappa$ deviates strongly from 20 pW/K, the claimed $ZT$ would need revision.
- The gap between measured $ZT$ (0.4) and computed $ZT$ (0.8) suggests contact geometry and the fraction of Fe(III) junctions are the main controllable variables.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports a combined experimental and computational study of Au|organometallic|Au junctions based on trans-Ru(dppe)2 bis(arylacetylide) and trans-Fe(cyclam) bis(arylacetylide) complexes. The authors perform MCBJ single-molecule conductance measurements at 4.2 K with IETS verification, C-AFM conductance and Seebeck measurements on SAMs at room temperature, SThM thermal-conductance measurements of the same SAMs, and DFT-NEGF transmission calculations. From these data they estimate thermoelectric power factors and ZT values, concluding that Fe-cyclam junctions can reach an experimental ZT up to 0.4, which they state is among the highest reported for molecular junctions at room temperature.
Significance. The strengths are the breadth and care of the characterization: the same molecules are studied as single molecules and as SAMs, with independent measurements of S, G, and κ on multiple samples, and the computational section includes explicit tests of electrode shape, geometry optimization, exchange-correlation functional, and counterion treatment. If the ZT claim were internally consistent, the observation of S≈130–145 μV/K with low per-molecule thermal conductance in a stable room-temperature organometallic junction would be a significant contribution to molecular thermoelectrics. The MCBJ-IETS identification of molecular junctions, the assignment of two Fe oxidation states using three independent techniques, and the explicit comparison of single-molecule and SAM data are also valuable. However, the headline ZT value is not a directly measured quantity and currently rests on an internal numerical inconsistency in the Discussion; this must be fixed before the central claim can be assessed.
major comments (4)
- [Discussion (ZT estimation paragraph) and Abstract] The stated ZT values are internally inconsistent with the per-molecule thermal conductances given in the same paper. The SThM section (Fig. 7 and surrounding text) and the Discussion assign κ = 50 pW/K to Ru-dppe and κ = 20 pW/K to Fe-cyclam. However, the reported numbers reproduce only with the opposite assignment. Using ZT = S²GT/κ with S = 28 μV/K, G(HC, Ru) = 6.9×10⁻⁴ G0 gives ZT ≈ 2.5×10⁻⁴ for κ = 50 pW/K but ZT ≈ 6.3×10⁻⁴ for κ = 20 pW/K, matching the reported 6.2×10⁻⁴. For Fe-cyclam with S = 145 μV/K and Gmax = 4.2×10⁻² G0, the reported 0.4 corresponds to κ = 50 pW/K, whereas the stated κ = 20 pW/K gives ZT ≈ 1.0. The Ru Gmax value (reported 6.6×10⁻³) likewise uses κ = 20 pW/K. No combination of the stated inputs reproduces the reported Fe HC value of 0.01 with S = 145 μV/K (one obtains ≈0.015 for κ = 50 pW/K and ≈0.037 for κ = 20 pW/K). Because the abstract's central claim of 'experimental ZT ... reaches up to 0.4' and the conclusion's 'experimental value 0.4' rest on these calculations, the authors must correct the assignment of thermal conductances, explicitly state whether electronic thermal conductance is included (at Gmax the Wiedemann–Franz contribution would give total κ ≈ 44 pW/K and ZT ≈ 0.47), or revise the reported ZT values.
- [SI, Null-point scanning thermal microscopy; Discussion] The per-molecule thermal conductances used in the ZT estimate are not directly measured but are obtained by dividing the SAM thermal conductance Gth(SAM) ≈ 20 nW/K by an assumed number of contacted molecules (≈400 for Ru-dppe and ≈10³ for Fe-cyclam). That number is estimated from a Hertzian contact radius of ≈20 nm and an assumed area per molecule. The resulting ZT scales inversely with the assumed N, so a factor-of-two uncertainty in N changes the headline value by the same factor. The estimate also assumes that heat transport through the SAM is additive (independent parallel molecules). The manuscript should propagate this uncertainty and should justify applying the SAM-average per-molecule thermal conductance to the high-conductance single-molecule MCBJ junctions whose G is used in the ZT formula, rather than presenting ZT = 0.4 as a directly measured quantity.
- [Computational study and Figure S40] The computed ZT = 0.8 is not an independent prediction in its current form. The calculation uses the experimentally determined κph = 20 pW/K as an input, and the authors state that absolute values of the Seebeck coefficient are not meaningful because they depend on the arbitrary Fermi-level shift ΔEF, yet they report a single ZT = 0.8 without specifying the corresponding ΔEF and without giving the spread of values shown in Figure S40. The calculated conductance also overestimates several of the MCBJ conductance values by orders of magnitude. The claimed agreement between the computed ZT (0.8) and the experimental estimate (0.4) should be reassessed after correcting the experimental arithmetic identified above.
- [Experimental methods and Discussion] The ZT is presented as a room-temperature figure of merit, but it combines the MCBJ high-conductance G measured at 4.2 K in vacuum with Seebeck and thermal-conductance values measured on SAMs at room temperature. The manuscript provides no evidence that the 4.2-K single-molecule conductance equals the room-temperature value, and the room-temperature C-AFM SAM conductance is orders of magnitude lower than the MCBJ Gmax used in the ZT estimate. If the ZT claim is intended to apply at room temperature, this assumption needs to be stated and justified, or a room-temperature conductance should be used.
minor comments (4)
- [Introduction (page 5)] The phrase 'quantum computing calculations are completing the study' should read 'quantum-chemical calculations complete the study'; the paper uses DFT/NEGF, not quantum computing.
- [SI, Table S2] In Table S2, 'mmersion' is a typo for 'immersion'.
- [SI, Null-point scanning thermal microscopy] In the sentence 'we estimate by the same method that ca. 400 (103, respectively) Ru-dppe (Fe-cyclam, respectively) molecules are contacted', the '103' should be typeset as 10³ for clarity.
- [Throughout] The acronym EDMS is defined as 'ethyl(dimethyl)silane', but the synthetic scheme uses a 2-(trimethylsilyl)ethyl protecting group; please check the nomenclature and use a consistent abbreviation.
Circularity Check
The computed ZT is partly calibrated by a free Fermi-shift choice, and the reported experimental Fe-cyclam ZT=0.4 is internally inconsistent with the stated 20 pW/K thermal conductance, though the underlying S, G, and kappa measurements are independent.
-
fitted input called prediction
[Computational study, paragraph beginning 'The calculated Seebeck coefficients at 300 K...' and Discussion ZT comparison]
"The calculated Seebeck coefficients at 300 K are reaching values in the range of ∼150 μV/K for both [Fe-cyclam]0 and [Fe-cyclam]+ junctions and ∼90 μV/K for Ru-dppe junctions for a moderate shift of the Fermi level (by 0.025 eV). The order of magnitude and ordering of the Seebeck coefficient of Fe-cyclam versus Ru-dppe are in line with the experiments."
The transmission spectrum is shifted by a free parameter ΔEF, which the paper states is needed because 'a quantitative description of the level alignment is not achievable at our level of theory.' The paper also notes that S varies by up to 50% for 0.05 eV shifts, so absolute values are not meaningful. Choosing ΔEF=0.025 eV yields S≈150 μV/K, matching the measured Fe-cyclam S=145 μV/K that is used in the experimental ZT estimate. Since ZT∝S², the computed ZT=0.8 is partly calibrated to the experiment it is compared with; the agreement 'close to the computed value' is therefore not an independent first-principles confirmation of the experimental claim.
-
other
[Discussion, paragraph 'In our study, we perform these evaluations by combining several techniques...']
"we provide an experimental estimation of the ZT by considering ... the thermal conductance per molecule (50 pW/K and 20 pW/K – Ru-dppe and Fe-cyclam) ... If we consider instead the Gmax value measured in MCBJ ..., we can estimate a maximum ZT of 6.6 × 10-3 for Ru-dppe and 0.4 for Fe-cyclam."
Using the standard formula ZT=S²GT/κ at T=300 K, the reported numbers are reproduced only with the opposite per-molecule thermal conductances from those stated in the same paragraph. S=28 μV/K, G=7.3×10⁻³G0, and κ=20 pW/K give 6.6×10⁻³ (the reported Ru value), while S=145 μV/K, G=4.2×10⁻²G0, and κ=50 pW/K give 0.41 (the reported Fe value). The text assigns 50 pW/K to Ru-dppe and 20 pW/K to Fe-cyclam, so the headline Fe ZT=0.4 is not computed from Fe-cyclam's measured 20 pW/K but from Ru-dppe's 50 pW/K. With the stated Fe κ=20 pW/K, the same S and G give ZT≈1.03, so the central numerical claim is not reproducible from the paper's own inputs; this is an internal derivation-chain failure even though it is not a self-referential circularity.
full rationale
The experimental backbone is not self-definitionally circular: the Seebeck coefficient, electrical conductance, and thermal conductance are measured by independent techniques (C-AFM/MCBJ, MCBJ, and SThM), and no equation reduces one observable to another. The per-molecule thermal conductances are obtained by dividing a measured SAM value by an assumed molecule count, which is an assumption rather than a circular definition. No load-bearing self-citation chain is present; previous work by the same groups is used for synthesis, prior conductance values, and background, not to forbid alternatives or force the central claim. The genuine circularity concern is limited to the computational ZT: the Fermi-level shift is a free parameter chosen in a range that reproduces the measured Seebeck coefficient, so the computed ZT=0.8 inherits that calibrated S quadratically. Separately, the reported experimental Fe-cyclam ZT=0.4 is internally inconsistent with the stated 20 pW/K thermal conductance, matching instead Ru-dppe's 50 pW/K; this is a numerical reproducibility flaw rather than a circularity, but it undermines the precise headline value. The qualitative conclusion that Fe-cyclam outperforms Ru-dppe is robust to the swap, but the specific ZT=0.4 and the comparison with computed ZT=0.8 are not reliable as printed. Overall circularity score 4 reflects the partial calibration of the theoretical ZT and the broken derivation chain in the experimental ZT estimate, while acknowledging that the primary measurements themselves are independent and not defined in terms of each other.
Assumptions & free parameters
free parameters (4)
- DFT Fermi-level shift ΔEF =
scanned between -0.10 and +0.10 eV; ZT=0.8 at shift reproducing S≈150 μV/K
- Number of molecules under SThM tip =
≈1000 for Fe-cyclam; ≈400 for Ru-dppe
- Number of molecules under C-AFM tip =
N≈5 (Ru-dppe), N≈15 (Fe-cyclam)
- SLM parameters ε0, Γ1, Γ2 =
e.g., ε0=0.32±0.12 eV (Ru SAM); ε0=0.15/0.55 eV (Fe MCBJ)
assumptions (5)
- standard math Landauer formula and single-level Lorentzian transmission model (SLM, Eq. S4) describe the junction I-V and thermopower.
- domain assumption SAM behaves as independent molecules conducting in parallel for both electrical and thermal transport.
- domain assumption Thermal contact radius rth≈20 nm derived from a water meniscus model with contact angle 30°.
- domain assumption The Seebeck coefficient measured on the SAM is representative of the single high-conductance molecule in the MCBJ when computing ZT.
- ad hoc to paper DFT-GGA with NEGF and a compensating Br- counterion (with -3 eV potential shift and -10 eV Hubbard shift) adequately describes the [Fe-cyclam]+ charged junction.
Cite this review
Pith. "Pith review of Electronic and Thermoelectric Properties of Molecular Junctions Incorporating Organometallic Complexes: Implications for Thermoelectric Energy Conversion." pith.science (2026). https://pith.science/paper/LQHRN2K5
@misc{pith2026250621113,
author = {Pith},
title = {Pith review of: Electronic and Thermoelectric Properties of Molecular Junctions Incorporating Organometallic Complexes: Implications for Thermoelectric Energy Conversion},
year = {2026},
howpublished = {\url{https://pith.science/paper/LQHRN2K5}},
note = {Machine review of arXiv:2506.21113}
}
read the original abstract
The electronic and thermoelectric properties of molecular junctions formed from iron and ruthenium metal-acetylide were studied using complementary experimental techniques and quantum chemical simulations. We performed physical characterizations of single-molecule and self-assembled monolayer junctions of the same molecules that allowed meaningful comparisons between the Ru and Fe adducts. In the case of the Fe-containing junctions, two distinct oxidation states are present. These junctions exhibit one of the highest Seebeck coefficients (S ca. 130 {\mu}V/K) reported to date for similar systems paired with broad electric conductance distribution and limited thermal conductance. As a result, the experimental thermoelectric figure of merit ZT for Fe-containing junctions reaches up to 0.4 for junctions with relatively high conductance. This is one of the highest ZT values reported for molecular systems at room temperature.
Reference graph
Works this paper leans on
-
[1]
S17 Figure S11. XPS spectra of the S2p region: (a) Ru-dppe SAM, (b) Fe-cyclam SAM. The dark lines are the measurements. The blue lines are the fits with the two-doublet decomposition shown in orange. Figure S12-a shows the Ru3d3/2 peak observed for the Ru-dppe SAM. This peak is weak, as also noted in other previous results on Ru-dppe-type SAMs.22,23 In th...
-
[5]
T.; Rubio-Bollinger, G.; Lambert, C
(58) Evangeli, C.; Gillemot, K.; Leary, E.; González, M. T.; Rubio-Bollinger, G.; Lambert, C. J.; Agraït, N. Engineering the Thermopower of C 60 Molecular Junctions. Nano Lett. 2013, 13 (5), 2141–2145. https://doi.org/10.1021/nl400579g. (59) Malen, J. A.; Doak, P.; Baheti, K.; Tilley, T. D.; Majumdar, A.; Segalman, R. A. The Nature of Transport Variations...
-
[50]
Blue curve: transmission at the G point. The electronic transmission depends on the level energy and k-point in the transversal direction A and B of incoming electrons from the electrodes at a given energy (T(E, k)), which implies the possible need for a specific k-sampling for each energy. The time reversal symmetry allows us to do the calculation for on...
work page 2015
-
[858]
https://doi.org/10.1039/dt9760000858. (3) Cao, Z.; Forrest, W. P.; Gao, Y.; Fanwick, P. E.; Ren, T. Trans -[Fe(Cyclam)(C 2 R) 2 ] + : A New Family of Iron(III) Bis-Alkynyl Compounds. Organometallics 2012, 31 (17), 6199–6206. https://doi.org/10.1021/om300515r. (4) Lenfant, S.; Viero, Y.; Krzeminski, C.; Vuillaume, D.; Demeter, D.; Dobra, I.; Oçafrain, M.; ...
-
[1635]
https://doi.org/10.1039/dt9840001635. (2) Chan, P.-K.; Poon, C.-K. Structural and Mechanistic Studies of Co-Ordination Compounds. Part XII. Syntheses and Characterization of Some Dianiono(1,4,8,11-Tetra-Azacyclotetradecane)-Manganese(III), -Iron(III), and -Nickel(III) Salts. J. Chem. Soc., Dalton Trans. 1976, No. 10,
-
[1991]
(19) Rabe, J. P.; Knoll, W. An Ellipsometric Method for the Characterization of Macroscopically Heterogeneous Films. Optics Commun. 1986, 57 (3), 189–192. https://doi.org/10.1016/0030-4018(86)90415-3. (20) Mervinetsky, E.; Alshanski, I.; Lenfant, S.; Guerin, D.; Medrano Sandonas, L.; Dianat, A.; Gutierrez, R.; Cuniberti, G.; Hurevich, M.; Yitzchaik, S.; V...
-
[2009]
Conductance of Molecular Wires Connected or Bonded in Parallel
(35) Magoga, M.; Joachim, C. Conductance of Molecular Wires Connected or Bonded in Parallel. Phys. Rev. B 1999, 59 (24), 16011–16021. https://doi.org/10.1103/PhysRevB.59.16011. S50 (36) Reuter, M. G.; Hersam, M. C.; Seideman, T.; Ratner, M. A. Signatures of Cooperative Effects and Transport Mechanisms in Conductance Histograms. Nano Lett. 2012, 12 (5), 22...
-
[2011]
Nanoscale Thermal Mapping of Few-Layer Organic Crystals
(53) Zhang, Y.; Zhang, C.; Wei, D.; Bai, X.; Xu, X. Nanoscale Thermal Mapping of Few-Layer Organic Crystals. CrystEngComm 2019, 21 (36), 5402–5409. https://doi.org/10.1039/C9CE00827F. (54) Dryden, J. R. The Effect of a Surface Coating on the Constriction Resistance of a Spot on an Infinite Half-Plane. J. Heat Transf. 1983, 105 (2), 408–410. https://doi.or...
Show all 10 references
-
[2013]
(41) Beebe, J
https://doi.org/10.5451/UNIBAS-006119305. (41) Beebe, J. M.; Kim, B.; Gadzuk, J. W.; Daniel Frisbie, C.; Kushmerick, J. G. Transition from Direct Tunneling to Field Emission in Metal-Molecule-Metal Junctions. Phys. Rev. Lett. 2006, 97 (2), 026801. https://doi.org/10.1103/PhysR...
-
[2017]
A.; Chiechi, R
(16) Weiss, E. A.; Chiechi, R. C.; Kaufman, G. K.; Kriebel, J. K.; Li, Z.; Duati, M.; Rampi, M. A.; Whitesides, G. M. Influence of Defects on the Electrical Characteristics of Mercury-Drop Junctions: Self-Assembled Monolayers of n -Alkanethiolates on Rough and Smooth Silver. J...
2007 doi
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.