{"id":"1bbb2a08-c4a6-4ab2-ac43-b62d0ed8ca93","arxiv_id":"1908.02665","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A one-level model with both electron-vibration and electron-electron interactions reproduces the measured conductance and Seebeck coefficient of Au-C60-Au junctions at 100 K, where the bare one-level model fails.","lead":"This paper models a single C60 molecule between gold electrodes, including the vibration of the molecule and the repulsion between electrons, to explain measured conductance and thermopower. It shows that including both interactions together reproduces the experimental gate-voltage dependence better than a simple one-level model.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's 'only combined' claim is under-supported because the U-only control is absent; Eq. 24's two-peak ansatz is secondary and does not change the conditional verdict.","rationale":"The reader's weakest_assumption centers on Eq. 24; I agree that approximation deserves scrutiny, but I do not think it is the most load-bearing. The paper is claiming necessity, and necessity claims require an explicit ablation. The theory section states the validity condition for Eq. 24 as Gamma approximately EP << U; with U=0.3 eV, the Coulomb peak separation is an order of magnitude larger than the broadening, so the ansatz is at least plausible. The real gap is logical: no U-only baseline is shown, and the text's 'demonstrate' outruns the evidence in Fig. 3. This is not a charge of inconsistency; it is a request for a control experiment that the model itself supplies. The paper has genuine strengths: the self-consistent adiabatic approach is non-perturbative, and the thermal-conductance prediction in Fig. 4 is an independent falsifiable output. But the headline claim's force depends on the missing control. I therefore keep the reader's CONDITIONAL verdict and recommend requiring the U-only calculation, or an explicit citation to a prior calculation, before the 'only combined' claim is accepted.","tokens_in":14690,"tokens_out":10258,"duration_ms":125560,"concrete_test":"Compute G(VG) and S(VG) with U=0.3 eV and EP=0 (lambda=0 in Eq. 24), using exactly the Fig. 3 parameters: E0-mu=0.065 eV, Gamma=0.032 eV, alpha=0.006 eV/V, T=100 K. Compare the U-only peak position, height, and Seebeck curve with the experimental data in Fig. 3. Then scan E0-mu over at least 0.045-0.075 eV and, if needed, Gamma over 0.02-0.05 eV to check whether any U-only parameter set fits both G and S. If none does, the combined claim is supported; if one does, the paper's central assertion fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing concern is not the two-peak ansatz in Eq. (24) but the missing ablation control for the paper's own central claim, stated in the abstract and Sec. IV: 'only the combined effect of local electron-vibration and electron-electron interactions'. The comparison shown in Fig. 3 is coherent (E0-mu=0.065 eV), EP-only (EP=0.018 eV), and EP+U (U=0.3 eV). No U-only curve is presented anywhere. Since E0 was shifted from the experimental fit value 0.057 to 0.065 eV and EP and U are order-of-magnitude estimates, the 'only combined' claim requires showing that U alone fails over the plausible parameter range. The absence of this control leaves open a simpler alternative: a U-only (or U plus adjusted E0/Gamma) model might already place the conductance peak near VG about 5 V and produce an acceptable Seebeck curve, in which case the combined interaction is sufficient but not necessary. The Eq. (24) concern is secondary because U=0.3 eV is an order of magnitude larger than both Gamma=0.032 eV and EP=0.018 eV, so the two Coulomb peaks remain well resolved despite EP/Gamma being comparable.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a self-consistent adiabatic treatment of a one-level Anderson-Holstein model to describe charge conductance and Seebeck coefficient of Au-C60-Au molecular junctions at T = 100 K. The model includes a low-frequency center-of-mass vibrational mode with polaron energy EP and a local Hubbard repulsion U, and averages transport quantities over the oscillator distribution. The authors show that the coherent one-level model fitted to the Seebeck data of Ref. [17] fails to reproduce the conductance peak position, that adding electron-vibration coupling shifts and reduces the conductance peak, and that adding electron-electron interaction further reshapes the curves so that both G and S appear to match the experimental gate-voltage dependence. The paper also presents predictions for the electronic thermal conductance and the electronic figure of merit ZT_el, and concludes that only the combined effect of electron-vibration and electron-electron interactions can explain the experimental data.","tokens_in":14927,"tokens_out":4140,"duration_ms":47817,"significance":"If fully established, the result would be a valuable demonstration that a non-perturbative adiabatic method incorporating both electron-vibration coupling and Coulomb blockade can quantitatively capture two complementary transport observables in a molecular junction. The paper is commendable for using a well-defined one-level model, for clearly separating coherent, EP-only, and EP+U contributions, and for providing concrete predictions for Gel_K and ZT_el that are falsifiable by future experiments. However, the central 'only combined effect' claim is currently under-supported because no U-only control calculation is shown and because the quantitative agreement with the experimental data is asserted by visual inspection rather than by a goodness-of-fit measure or a parameter-sensitivity analysis. The contribution is therefore significant if the missing control calculations are supplied, but as it stands the conclusion is stronger than the evidence.","major_comments":[{"comment":"The abstract and Sec. IV claim that 'only the combined effect of local electron-vibration and electron-electron interactions' can reproduce both G and S, but the paper never presents a U-only calculation. The comparison in Fig. 3 contains only the coherent result, the EP-only result, and the EP+U result. Since the level alignment E0-mu is shifted from the experimental fit value of 0.057 eV to 0.065 eV, and since a finite U alone creates a second spectral peak at epsilon = -U and transfers spectral weight from the main peak, it is possible that a U-only model with an adjusted E0 or Gamma already places the conductance peak near VG ~ 5 V and produces an acceptable Seebeck curve. Showing that the U-only model fails over the full plausible parameter range is essential to support the 'only combined' necessity claim; without it, the combined model is merely sufficient, not necessary.","section":"Sec. IV, Fig. 3"},{"comment":"The quantitative basis of the 'very good agreement' statement is not established. The parameters EP = 0.018 eV, U = 0.3 eV, and the shifted E0 - mu = 0.065 eV are order-of-magnitude estimates or ad hoc adjustments, and no error metric, confidence interval, or sensitivity study is provided. The authors should report, for example, the root-mean-square deviation between theory and experiment for G and S in the relevant VG range, and show how the agreement degrades when EP, U, Gamma, and alpha are varied within physically motivated ranges. Without such an analysis, the reader cannot distinguish a robust physical description from an overfit to a single dataset.","section":"Sec. IV, Fig. 3 and Sec. III parameter values"},{"comment":"The spectral function in Eq. (24) is written as a weighted sum of two independent Lorentzians of equal width Gamma, and the text states that this approximation is valid when Gamma ~ EP << U. The chosen parameters give EP/Gamma = 0.018/0.032 ~ 0.56, so the stated scale separation is not well satisfied, even though U >> Gamma ensures that the two Coulomb peaks are well separated. The authors should assess the error introduced by the two-peak ansatz at EP/Gamma ~ 0.56, either by benchmarking against a numerically exact or established approximate method in this parameter regime, or by estimating the neglected interference and non-Lorentzian corrections. If those corrections are non-negligible, the computed G and S curves would not reliably reflect the underlying many-body model.","section":"Eq. (24) and parameter regime discussed after it"}],"minor_comments":[{"comment":"There is a typographical error in the title/abstract: 'thermoele ctric' should be 'thermoelectric'.","section":"Title and abstract"},{"comment":"The sentence 'only very the thermal conductance of single-molecule junctions has been fully characterized' appears to be missing a word; it should probably read 'only very recently has the thermal conductance...'.","section":"Sec. I, fourth paragraph"},{"comment":"The definition of the hybridization width matrix contains a redundant repetition: 'Γ m,n = ∑α Γ m,n α = ∑α Γ m,n α' should be a single sum with the same symbol.","section":"Sec. II, after Eq. (6)"},{"comment":"The axis labels 'SC60' and 'GC60' should be typeset as 'S_{C60}' and 'G_{C60}', and the unit of S is written as 'K/V' instead of 'V/K' in the caption of Fig. 1.","section":"Fig. 1 and Fig. 3 captions"},{"comment":"The text estimates EP ~ 0.030 eV from experimental parameters but then uses EP = 0.018 eV in the calculations; the reason for choosing the lower value should be stated explicitly, beyond the statement that increasing EP shifts the peak too much.","section":"Sec. III"},{"comment":"The phrase 'in the unities chosen in Figure 4' should be 'in the units chosen in Figure 4'.","section":"Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The core concern is that the paper's headline claim of necessity is not established by the evidence presented. The missing U-only control is a straightforward calculation that the authors can add without changing the model, and a sensitivity analysis would substantially strengthen the paper. The Eq. (24) concern is secondary given U >> Gamma, but it should be addressed for rigor. I would be willing to look at a revised version with these additions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"K.,\n\nWorth a look. The paper takes the adiabatic approach that these authors have been developing (electron-vibration coupling plus Coulomb blockade) and applies it to the Au-C60-Au junctions of Kim et al. The one-level Landauer fit to the Seebeck data puts the conductance peak at the wrong gate voltage; the EP+U calculation gets both G and S in the right ballpark. That simultaneous description is a genuine constraint, and the predicted electronic thermal conductance gives an independent check. The writing is clear and the parameter regime is explicit.\n\nThe soft spots are real but not fatal. The main one is missing the U-only control. The abstract and conclusions say that only the combined effect of electron-vibration and electron-electron interactions reproduces the data, but the paper shows only the EP-only and EP+U curves. A U-only calculation (even with the shifted E0 and a plausible Gamma) would close the logical loop. Without it, a U-only model might also place the conductance peak near VG ~ 5 V, and the combined model would be sufficient but not necessary. The other issue is parameter hand-tuning: E0 is shifted from 0.057 to 0.065 eV, and EP and U are chosen within physically motivated ranges. With three adjustable parameters and two curves, the agreement is encouraging, but no goodness-of-fit or sensitivity analysis is given. A referee should ask for that.\n\nThe Eq. (24) two-peak spectral ansatz is a more minor concern than the reader's report suggests. U = 0.3 eV is an order of magnitude above Gamma = 0.032 eV, so the two Coulomb peaks are well resolved; the condition stated in the paper (Gamma and EP both much less than U) holds. The equilibrium P(x) assumption is worth a robustness check, but it is not the load-bearing issue.\n\nOverall, the paper is a legitimate extension of the authors' earlier work, with a clear physical point and a useful falsifiable prediction. It deserves refereeing. I would ask for the U-only control and a few sensitivity tests before accepting, but the core approach looks sound.","headline":"A solid application of the authors' adiabatic method to C60 junctions, but the 'only combined' claim needs a U-only control before it is proven.","tokens_in":15503,"tokens_out":3772,"would_cite":false,"duration_ms":39293,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["73.63.-b","72.20.Pa","73.23.Hk"],"model":"deepseek-v4-flash","headline":"A one-level model with both electron-vibration and electron-electron interactions reproduces the measured conductance and thermopower of fullerene junctions, where the coherent one-level model fails.","keywords":["thermoelectricity","molecular junctions","fullerene C60","Seebeck coefficient","Coulomb blockade","electron-vibration coupling","adiabatic approximation","self-consistent transport"],"falsifier":"Measure the full gate-voltage trace from about -40 V to +60 V at T = 100 K in an Au-C60-Au junction and look for the predicted secondary conductance peak and oscillatory Seebeck signal in the negative-gate (high-energy) region; their absence, or a peak position inconsistent with U = 0.3 eV, would falsify the two-peak weighted spectral ansatz.","tokens_in":14437,"feed_emoji":"⚡","tokens_out":8298,"duration_ms":76771,"temperature":0.7,"pith_summary":"The paper claims that the measured gate-voltage dependence of both the electrical conductance and the Seebeck coefficient in Au-C60-Au molecular junctions can be explained only when local electron-vibration and electron-electron interactions act together. Using a self-consistent adiabatic treatment of the molecule's center-of-mass vibration and a Coulomb blockade description of the local charging energy, the authors reproduce the experimental conductance peak near V_G ~ 5 V and the thermopower curve at T = 100 K. This matters because the standard coherent one-level model, fit to the thermopower data, places the conductance peak near 9 V, in clear disagreement with experiment. The result identifies the combined many-body interactions, not either alone, as the mechanism controlling charge and thermoelectric transport through large molecules.","feed_headline":"Two local interactions together explain fullerene transport","feed_subtitle":"Coulomb repulsion plus vibration coupling reproduces both measured signals; the simple one-level model cannot.","key_machinery":"The central object is the displacement-averaged electronic spectral function A(E), built from Eq. (24): for a fixed oscillator displacement x, it is a weighted sum of two Lorentzian peaks at energies epsilon + lambda x and epsilon + lambda x + U, with weights 1 - rho(x) and rho(x), where rho(x) is the self-consistently determined level occupancy per spin. This two-peak ansatz incorporates Coulomb blockade into the adiabatic approximation, in which the slow center-of-mass mode is treated as a classical field and transport coefficients are obtained by averaging over its equilibrium position distribution P(x). The mechanism carries the whole argument: the electron-vibration coupling shifts spectral weight and the Hubbard repulsion transfers weight to the second peak, jointly reshaping G and S.","core_discovery":"Within a single-level Anderson-Holstein-type model, the authors show that the conductance gap, peak position, and thermopower magnitude observed in gated fullerene junctions require the simultaneous presence of the electron-vibration coupling EP = 0.018 eV and a Hubbard repulsion U = 0.3 eV. The electron-vibration coupling alone shifts and narrows the conductance peak, and the Hubbard term suppresses the conductance amplitude and generates a secondary Coulomb-blockade feature. Together they bring the calculated conductance peak and Seebeck zero into agreement with the experimental data at 100 K for a level position E0 - mu = 0.065 eV, where the coherent model with E0 - mu = 0.057 eV fails. The paper also predicts an electronic thermal conductance whose gate-voltage profile closely follows the charge conductance.","pith_inferences":["The same two-Lorentzian weighted ansatz could be applied to other large molecules with a soft center-of-mass mode, predicting that a Coulomb blockade satellite peak should appear whenever the charging energy is comparable to the vibrational shift and the mode remains adiabatic.","Because the averaging uses the equilibrium oscillator distribution, the fit at 100 K implicitly assumes current-induced heating and forces are negligible; at higher bias or lower temperature this assumption should break down, offering a testable deviation.","The parameter set implies an effective level shift from vibrational coupling of about EP, so gate-voltage calibration in coherent fits may systematically underestimate the LUMO-to-chemical-potential distance.","Extending the comparison to the full gate-voltage range would discriminate the two-peak ansatz from alternative multi-level interference models, since the predicted secondary peak is a distinctive Coulomb-blockade signature."],"forward_implications":["For fullerene junctions at 100 K, the coherent one-level model cannot fit both conductance and thermopower; the combined interactions are the minimal required ingredient.","The conductance peak position is controlled by EP and the peak amplitude is reduced by U, so fits to G alone without Coulomb repulsion will mislocate the level.","The Seebeck coefficient stays relatively robust near resonance, meaning thermopower data alone cannot discriminate many-body mechanisms; conductance data are the discriminating probe.","The model predicts a secondary conductance peak and oscillatory thermopower at gate voltages far from resonance, signatures of Coulomb blockade that could be searched for experimentally.","The predicted electronic thermal conductance is a fraction of the thermal conductance quantum and tracks the charge conductance, setting a scale for future single-molecule thermal measurements."],"supporting_citations":[{"why":"Supplies the experimental gate-voltage traces of G and S for Au-C60-Au junctions at 100 K and the coherent one-level fit parameters.","marker":"17"},{"why":"Provides the self-consistent adiabatic formalism for charge and heat transport in soft molecular junctions, including the oscillator distribution averaging.","marker":"39"},{"why":"Extends the adiabatic approach to include strong Coulomb repulsion through the two-peak spectral function ansatz used in Eq. (24).","marker":"40"},{"why":"Gives the experimental basis for the C60 center-of-mass mode (about 5 meV) and the charging energy scale (about 0.27 eV).","marker":"23"},{"why":"Supplies the fullerene junction context and the estimate that the polaron energy EP is on the order of 0.030 eV.","marker":"10"},{"why":"Supports treating C60 transport through a single dominant level by showing the LUMO levels are split by a few tenths of an eV.","marker":"44"}],"fun_headline_variants":["Fullerene thermopower needs both vibrations and repulsion","Combined interactions solve fullerene junction puzzle","Pair of local effects drive fullerene thermoelectrics","Single-level model works only with electron and phonon terms","Vibration plus Coulomb gives right thermoelectric signals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the strongly interacting spectral function can be represented, at every oscillator displacement, as a weighted sum of two independent Lorentzian peaks whose occupation weight is self-consistently determined, a two-peak ansatz that also assumes the oscillator distribution stays in equilibrium; the fitted parameters, with Gamma = 0.032 eV and EP = 0.018 eV, sit at the edge of the regime where that ansatz is expected to hold.","fun_headline_variants_meta":{"raw":{"variants":["Fullerene thermopower needs both vibrations and repulsion","Combined interactions solve fullerene junction puzzle","Pair of local effects drive fullerene thermoelectrics","Single-level model works only with electron and phonon terms","Vibration plus Coulomb gives right thermoelectric signals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000197,"raw_usage":{"total_tokens":1316,"prompt_tokens":846,"completion_tokens":470,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":462,"completion_tokens_details":{"reasoning_tokens":394}},"tokens_in":462,"tokens_out":470,"duration_ms":4597,"temperature":1.0,"reasoning_tokens":394,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:38:00.247531+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the full gate-voltage trace from about -40 V to +60 V at T = 100 K in an Au-C60-Au junction and look for the predicted secondary conductance peak and oscillatory Seebeck signal in the negative-gate (high-energy) region; their absence, or a peak position inconsistent with U = 0.3 eV, would falsify the two-peak weighted spectral ansatz.","supporting_citations":[{"cited_title":"Electrostatic control of thermoelectricity in molecular junctions","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental gate-voltage traces of G and S for Au-C60-Au junctions at 100 K and the coherent one-level fit parameters."},{"cited_title":"Charge and heat transport in soft nanosystems in the presence of time-dependent perturbations","cited_arxiv_id":null,"evidence_quote":"Provides the self-consistent adiabatic formalism for charge and heat transport in soft molecular junctions, including the oscillator distribution averaging."},{"cited_title":"A.; Ninno, D.; Cataudella, V","cited_arxiv_id":null,"evidence_quote":"Extends the adiabatic approach to include strong Coulomb repulsion through the two-peak spectral function ansatz used in Eq. (24)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the experimental basis for the C60 center-of-mass mode (about 5 meV) and the charging energy scale (about 0.27 eV)."},{"cited_title":"Thermoelectric efficiency of molecular junctions","cited_arxiv_id":null,"evidence_quote":"Supplies the fullerene junction context and the estimate that the polaron energy EP is on the order of 0.030 eV."},{"cited_title":"H.; Louie, S","cited_arxiv_id":null,"evidence_quote":"Supports treating C60 transport through a single dominant level by showing the LUMO levels are split by a few tenths of an eV."}],"review_version":1}