{"id":"c462de08-e1e8-4c88-b133-426633ef8301","arxiv_id":"1909.01196","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A generalized Landauer formula including strong vibronic decoherence explains the high, temperature-independent, distance-independent conductance observed in protein junctions.","lead":"This paper derives a new formula for electrical conductance through proteins attached to metal electrodes, combining coherent tunneling with strong vibrational decoherence. It aims to explain why some proteins conduct electricity over long distances with little temperature or distance dependence, contrary to standard theory.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Strong-decoherence formula (66) rests on an unverified Γ≪|R| hierarchy; without an estimate of vibrational relaxation rates in proteins, the claimed distance and temperature independence are not established.","rationale":"The reader's weakest_assumption identifies exactly the same concern: the derivation of Eq. (66) assumes that the inverse Redfield operator is dominated by the perturbatively shifted zero eigenvalue, requiring all other relaxation eigenvalues to be large compared to the lead couplings. The paper postulates Γ≪|R| but supplies no estimate of the vibrational relaxation rates for proteins. I agree that this is the most load-bearing point. My read sharpens it: the parameter values used to obtain nS-scale conductance (Γ ~ 0.1 eV) may already put the system outside the strong-decoherence regime if realistic protein decoherence rates are in the meV range, and no internal check or external estimate is provided. The concern is not that the derivation is internally inconsistent within the stated limit, but that the applicability to the experiments is unsupported. The paper does have independent support: the derivation is transparent, it recovers the coherent Landauer-Büttiker and Marcus limits, and the sum rules in Appendices A and D are consistent. The conductance magnitude in Section 5.5 arises from a parameter-free estimate using plausible Γ values, which is a positive feature. However, the missing spectral-gap estimate is the crux: without it, the explanation of the stylized facts rests on an unverified asymptotic condition. The reader's CONDITIONAL verdict already reflects this uncertainty, so no change in verdict is needed; the concrete test would provide the missing evidence.","tokens_in":17969,"tokens_out":12986,"duration_ms":132774,"concrete_test":"Compute the smallest nonzero eigenvalue magnitude of the Redfield operator L0 for a representative protein (e.g., the Streptavidin extended-Hückel orbitals from Section 6.3) using a realistic vibrational spectral density (reorganization energy λ ≈ 0.1-1 eV, cutoff ω_c ≈ 50-200 meV) and compare it with the Γ ≈ 0.1 eV values used in Section 5.5. If the smallest |λ_i| is not at least several times Γ, the approximation in Eq. (87) fails and the distance/temperature predictions of Eq. (66) are not supported in the claimed regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central formula (66) is derived in the strong-decoherence limit via Appendix B, where Eq. (87) approximates L^{-1} by the reciprocal of the perturbatively shifted zero eigenvalue of L0. This is valid only if all other eigenvalues of L0 are large compared to the lead couplings Γ_n; Section 5.2 frames this as Γ≪|R|. The paper does not estimate these eigenvalues for proteins. This is load-bearing because the stylized facts (distance independence, weak temperature dependence) are consequences of Eq. (66), not of the coherent limit Eq. (52), where distant localized orbitals would give exponentially suppressed conductance. Section 5.5 applies the formula with Γ in the 0.1 eV range, while realistic vibrational relaxation/decoherence rates in proteins, estimated from electron-phonon couplings and reorganization energies, are typically 1-100 meV. If |R| is not at least several times Γ, the perturbation expansion in Appendix B fails and Eq. (66) does not follow. Moreover, the Redfield tensor itself is a weak-coupling (Born-Markov) object; using it to describe strong decoherence is internally in tension unless the large |R| still emerges from a second-order treatment. The paper offers no spectral-gap estimate, no convergence check, and no comparison with an exact or higher-order calculation in a representative model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives a generalized Landauer-Büttiker formula for molecular conduction starting from a Markovian Redfield master equation. After obtaining a general expression for conductance in terms of the inverse of the Redfield operator (Eq. (28)), the authors consider two limits: weak contacts, where conductance is proportional to the electron-transfer rate (recovering Nitzan-type relations), and strong contacts. The central result is the strong-decoherence limit, Eq. (66), in which the conductance depends only on the equilibrium Boltzmann distribution and the electrode coupling strengths, leading to nS-scale conductance, weak temperature dependence, and distance independence when both contacts are strong. The paper applies this formula to Myoglobin and Cytochrome C temperature data, to the PNA distance-exponent relation, and to the conductance distribution of Streptavidin.","tokens_in":18265,"tokens_out":7595,"duration_ms":72191,"significance":"If correct, Eq. (66) is a significant simplification: it replaces a many-parameter microscopic description with a handful of macroscopic parameters and offers a concrete mechanism (vibrational mixing) for the experimentally observed high, temperature-insensitive, distance-insensitive conductance of proteins. The derivation from the Redfield equation is transparent, and the appendix sum rules (Eqs. (80) and (98)) are useful. The paper also provides falsifiable predictions, such as distance independence for two strong contacts, and states that simulation code is provided as supplementary material. However, the validity of the strong-decoherence reduction is not quantitatively established, and the experimental comparisons are mostly consistency checks with fitted parameters, so the current evidence for the central claim is incomplete.","major_comments":[{"comment":"The approximation (87) for L^{-1} retains only the perturbatively shifted zero eigenvalue of L0, which is valid only if all other eigenvalues of L0 (the vibrational relaxation rates) are large compared with the lead couplings Γ_n. The paper never estimates these eigenvalues for proteins; Section 5.5 uses Γ in the 0.1 eV range, so the assumed hierarchy Γ ≪ |R| is a postulate rather than an established condition. This is load-bearing because Eq. (66), and the claimed distance and temperature independence, follow only in this strong-decoherence limit. The authors should provide a quantitative estimate of |R| from electron-phonon coupling parameters or a numerical check of the spectral gap in a representative model.","section":"Appendix B, Eq. (87); Section 5.2"},{"comment":"The explanation of the PNA exponent β_G ≈ 0.66 β_ET is obtained by choosing x ≈ 0.17 after the fact; Eq. (75) with a fitted x that matches the measured exponent is a post-hoc fit, not a prediction of the theory. Moreover, the assumption that couplings scale as Γ_HOMO/LUMO ∝ e^{-β_ET x l} introduces the electron-transfer decay exponent β_ET as an input. The section should be reframed as a consistency check or as a demonstration that a plausible x can account for the data, with x treated as a free fit parameter.","section":"Section 6.2, Eq. (75)"},{"comment":"The fits in Figs. 1 and 2 use two adjustable parameters per curve (I0 and IT), and in the electrostatic Cytochrome C case the activation energy is taken from the experimental value 0.105 eV rather than from the computed orbital structure; the computed HOMO-LUMO gap/2 = 0.161 eV is not used. The agreement is therefore weaker evidence than the text suggests, and the activation-energy identification is partly circular. The authors should report goodness-of-fit, parameter uncertainties, and ideally make predictions for at least one curve with parameters fixed by the theory.","section":"Section 6.1, Tables 1 and 3"},{"comment":"There is an internal tension between using the Redfield equation, which is derived under a Born-Markov weak system-bath coupling assumption, and the strong-decoherence regime Γ ≪ |R|, which requires large R. The paper should discuss the parameter regime in which both assumptions hold, or demonstrate with a specific model that the Redfield tensor can produce relaxation rates exceeding Γ without leaving its domain of validity.","section":"Section 5.2"}],"minor_comments":[{"comment":"The phrase 'Landauer fromula' in the keywords is a typo for 'Landauer formula'.","section":"Abstract and Keywords"},{"comment":"The text contains 'LOMO' in place of 'LUMO' in the sentence about the smallest orbital energy differences.","section":"Section 5.3"},{"comment":"'Matlab and Phyton' should be 'Matlab and Python'.","section":"Section 8"},{"comment":"The authors state that start and end points of curves were visually extracted from original figures; they should indicate the estimated digitization uncertainty and provide the raw extracted data points in the supplementary material.","section":"Figures 1 and 2"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a real attempt at a new result, and Eq. (66) is worth taking seriously, but the paper overclaims on the experimental side and the key decoherence assumption is unchecked. I would send it to review, but the referee should push hard on the spectral gap problem.\n\nThe genuinely new thing is the strong-decoherence limit of a Redfield-generalized Landauer formula. For a closed system the Redfield Liouvillian has a zero eigenvalue corresponding to Boltzmann equilibrium; the paper perturbs that eigenvalue in the lead couplings and gets a conductance that depends only on Boltzmann weights and the Γ's, with vibrational details dropping out. That is a surprising and potentially important result, and it recovers the right coherent limit when |R|≪Γ. This is not circular: the formula is derived from stated physical assumptions, not from the data. They also provide code, so the curve in Fig. 3 is reproducible.\n\nThe soft spots are real. The perturbation step in Appendix B keeps only the shifted zero eigenvalue and assumes all other eigenmodes of L0 are fast. The paper never estimates those relaxation rates for a protein, and the parameter range used later (Γ≈0.1 eV) makes the required hierarchy Γ≪|R| far from obvious. There is also an internal tension: the Redfield tensor is itself a weak-coupling object, so \"strong decoherence\" has to come out of a Born-Markov term that is assumed large. Without a spectral-gap estimate or a check against a solvable model, Eq. (66) is conditional.\n\nThe experimental sections do not fix that. Temperature curves are reconstructed from two endpoints and fitted with I0/IT; the distance-exponent argument selects x=0.17 to match the measured ratio; the conductance distribution uses Gaussian couplings with means and variances set by visual optimization. Those are consistency checks, not validation. The authors mostly say that too, but the abstract and discussion lean on the experimental agreement.\n\nWho should read it: people working on decoherence in molecular junctions and protein transport. They will find the strong-decoherence formula a useful starting point even if the quantitative claims are not yet solid. I would not cite (66) as established, but I would cite the paper as the origin of the formula.\n\nRecommendation: send to peer review, with the requirement that the authors estimate or bound the spectral gap of the Redfield operator in a realistic model and present the experimental fits as exploratory.","headline":"A genuinely new strong-decoherence conductance formula that would be important if its unverified relaxation-rate assumption holds; the experimental claims outrun the evidence.","tokens_in":18813,"tokens_out":2265,"would_cite":false,"duration_ms":24496,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Protein-junction conductance is set by vibrations and Boltzmann weights, not by tunneling through the molecule.","keywords":["Landauer formula","conductance of biomolecules","metallic contacts","decoherence","Redfield equation","electron transfer","protein conductance","stylized facts"],"falsifier":"Compute the eigenvalues of the isolated-molecule Redfield operator $L_0$ from a realistic electron–phonon model of a protein junction. If any nonzero eigenvalue is of the same order as or smaller than $\\Gamma_n$, the strong-decoherence formula (66) is not justified. Experimentally, a direct contradiction would be a symmetric two-strong-contact junction whose conductance decays exponentially with electrode separation or follows an Arrhenius factor with half the HOMO–LUMO gap.","tokens_in":17784,"feed_emoji":"🧬","tokens_out":7872,"duration_ms":78248,"temperature":0.7,"pith_summary":"The paper sets out to show that the puzzling experimental behavior of protein junctions—conductance in the nano-Siemens range, almost no drop over 1–10 nm, and near constancy between 30 and 300 K—follows from a single generalized Landauer formula once strong coupling to vibrations is included. In the regime the authors call strong decoherence, the conductance is governed not by tunneling through the whole molecule but by the equilibrium Boltzmann occupation of molecular orbitals together with the electrode coupling strengths $\\Gamma_n$. The result reduces the problem from a full quantum-chemical transport calculation to a handful of macroscopic parameters, and it recovers the familiar Landauer and weak-coupling electron-transfer limits when decoherence is absent or when contacts are weak. A sympathetic reader would care because it offers a parameter-light explanation of why long protein wires can conduct as well as they do.","feed_headline":"Protein conductance obeys a vibration-driven formula","feed_subtitle":"It predicts nS currents that barely change with temperature or electrode distance when contacts bind strongly.","key_machinery":"The engine of the derivation is the Bloch–Redfield Liouvillian $L_{nmkl}$, whose inverse feeds the generalized current formula (28) through terms like $L^{-1}_{kknn}$. In the strong-decoherence limit the paper evaluates that inverse by perturbation theory: the zero eigenvalue of the isolated molecule's Redfield operator is shifted by the contacts to a small value, and the inverse is dominated by the reciprocal of that shift, $$$L^{{-1}}$_{nnmm} \\approx -\\frac{\\hbar\\,\\$varrho^{0}$_{nn}}{\\sum_p \\Gamma_p \\$varrho^{0}$_{pp}},$$ with $\\varrho^0$ the Boltzmann density matrix. This single identity is what removes the microscopic detail: it turns the many couplings of the molecule to its vibrations into the statement that the electron or hole wanders through the Boltzmann-weighted orbitals and exits through whichever electrode offers the larger escape rate.","core_discovery":"The paper's central claim is that in the strong-decoherence regime, the conductance of a protein junction is given by Eq. (66), $$G = \\frac{$e^{2}$}{h}\\left(T_h^R(E_F)P_h^L(T) + T_h^L(E_F)P_h^R(T) + T_e^R(E_F)P_e^L(T) + T_e^L(E_F)P_e^R(T)\\right),$$ where each $T$ is a tunneling term built from the orbital energies and level broadenings $\\Gamma_n$, and each $P$ is an equilibrium Boltzmann probability that a charge leaves through a given electrode. The formula is derived from the Redfield equation by inverting its Liouvillian in the limit where vibrational relaxation dominates the electrode couplings, and the derivation shows that all microscopic details of the electron–vibration interaction cancel out. The remaining input is just the molecular orbital spectrum and the two contact strengths. From this the paper obtains conductance values in the $0.01$–$10$ nS range, almost no temperature dependence below the nearest orbital gap, and no systematic dependence on the distance between two strong contacts.","pith_inferences":["Inference: if the strong-decoherence formula is right, isotopic substitution or matrix changes that alter vibrational relaxation rates without changing orbital energies should barely affect the conductance, since the vibration details cancel; this is a testable extension the paper does not perform.","Inference: the formula effectively describes the junction as a classical random walk on Boltzmann-weighted orbitals with electrode escape as the absorbing step, so the same logic could connect this result to environment-assisted transport models used in photosynthetic energy transfer.","Inference: the assumption $\\Gamma_n \\ll |R|$ could be checked directly by computing the full Redfield spectrum of a realistic electron–phonon model of a protein; real proteins may sit in the strong-decoherence regime or in a crossover where Eq. (66) needs correction terms."],"forward_implications":["With two strong contacts, conductance through a single protein should stay flat as electrode separation grows from 1 to 10 nm, because none of the terms in (66) depends systematically on distance.","With one strong and one weak contact, distance dependence can reappear with an exponent set by where the HOMO and LUMO sit relative to the strong contact, matching the power-law relation observed in peptide nucleic acids.","The electron-transfer rate and the conductance are not proportional in the strong-decoherence regime: the transfer rate keeps an Arrhenius factor and decays with length, while conductance does not.","The temperature dependence, when present, should be set by the smallest HOMO/HOMO-1 or LUMO/LUMO+1 gap, not by half the HOMO-LUMO gap, which is how the Myoglobin and Cytochrome C datasets are reproduced.","Predictions for new protein junction experiments can be made from orbital energies and contact couplings computed semiempirically, avoiding expensive full quantum transport calculations."],"supporting_citations":[{"why":"It provides the molecular-junction Landauer/Breit-Wigner derivation whose finite-temperature, decoherence-aware generalization is the paper's central object.","marker":"[20, 21]"},{"why":"It establishes the proportionality between conductance and electron-transfer rate that the paper recovers in the weak-contact limit and shows breaks down for strong contacts.","marker":"[13, 14, 15]"},{"why":"It gives the sequential-hopping electron-transfer relation recovered here when contacts are weak and the partition function is dominated by a single barrier.","marker":"[24]"},{"why":"It supplies the Redfield master equation whose inverse Liouvillian carries the strong-decoherence calculation.","marker":"[22]"},{"why":"It reports Myoglobin current-temperature data that the paper reproduces with the temperature-dependent part of the formula.","marker":"[25]"},{"why":"It reports Cytochrome C current-temperature data reproduced with activation energies identified from the orbital spectrum.","marker":"[26]"},{"why":"It reports the distribution of single-protein conductance data that the paper compares with simulations using Eq. (66).","marker":"[9]"},{"why":"It reports the power-law relation between electron-transfer rate and conductance that the paper explains via one strong and one weak contact.","marker":"[27]"}],"fun_headline_variants":["Protein conductance formula ignores distance and temperature","Simple formula predicts protein conductance in nS range","Landauer formula extended to protein transport","Protein junction conductance: new Landauer variant","Formula explains temperature and distance independence in proteins"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the vibrational relaxation rates inside the protein are all much larger than the electrode coupling strengths $\\Gamma_n$, so the inverse Redfield operator is dominated by a single perturbed zero eigenvalue; the paper does not compute or estimate those relaxation rates in actual proteins.","fun_headline_variants_meta":{"raw":{"variants":["Protein conductance formula ignores distance and temperature","Simple formula predicts protein conductance in nS range","Landauer formula extended to protein transport","Protein junction conductance: new Landauer variant","Formula explains temperature and distance independence in proteins"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000234,"raw_usage":{"total_tokens":1497,"prompt_tokens":948,"completion_tokens":549,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":564,"completion_tokens_details":{"reasoning_tokens":484}},"tokens_in":564,"tokens_out":549,"duration_ms":6475,"temperature":1.0,"reasoning_tokens":484,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:25:23.054623+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the eigenvalues of the isolated-molecule Redfield operator $L_0$ from a realistic electron–phonon model of a protein junction. If any nonzero eigenvalue is of the same order as or smaller than $\\Gamma_n$, the strong-decoherence formula (66) is not justified. Experimentally, a direct contradiction would be a symmetric two-strong-contact junction whose conductance decays exponentially with electrode separation or follows an Arrhenius factor with half the HOMO–LUMO gap.","supporting_citations":[{"cited_title":"The relationship between electron transfer rate and molecular conduction 2","cited_arxiv_id":null,"evidence_quote":"It gives the sequential-hopping electron-transfer relation recovered here when contacts are weak and the partition function is dominated by a single barrier."},{"cited_title":"The theory of open quantum systems","cited_arxiv_id":null,"evidence_quote":"It supplies the Redfield master equation whose inverse Liouvillian carries the strong-decoherence calculation."},{"cited_title":"Protein electronic conductors: hemin–substrate bonding dictates transport mechanism and efficiency across myoglobin","cited_arxiv_id":null,"evidence_quote":"It reports Myoglobin current-temperature data that the paper reproduces with the temperature-dependent part of the formula."},{"cited_title":"Solid-state electron transport via cytochrome c depends on electronic coupling to electrodes and across the protein","cited_arxiv_id":null,"evidence_quote":"It reports Cytochrome C current-temperature data reproduced with activation energies identified from the orbital spectrum."},{"cited_title":"Electronic decay length in a protein molecule","cited_arxiv_id":null,"evidence_quote":"It reports the distribution of single-protein conductance data that the paper compares with simulations using Eq. (66)."},{"cited_title":"The single-molecule conductance and electrochemical electron-transfer rate are related by a power law","cited_arxiv_id":null,"evidence_quote":"It reports the power-law relation between electron-transfer rate and conductance that the paper explains via one strong and one weak contact."}],"review_version":1}