{"id":"810125fd-2f18-4fc1-a7c4-500145de31cb","arxiv_id":"2505.23636","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":10,"one_line_summary":"For a model donor-acceptor rectifier, Fisher information shows the acceptor energy is the most estimable parameter with an optimal time, the donor energy is weakly estimable, and the vibrational frequency is best estimated at steady state and low temperature.","lead":"A theoretical study computes how precisely three parameters of a donor-acceptor molecular rectifier can be estimated from the dynamics of its vibrational mode. It finds the acceptor energy is the most measurable parameter with an optimal time, while the vibration frequency is best measured in the steady state at low temperature.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's Fisher-information claims depend on an underived master equation (Eqs. (8)-(9)) whose reduced system is ambiguous: Eq. (7) treats the vibrational TLS as the system, yet p1,p2 are called donor/acceptor occupations, and H_int does not match H_I in Eq. (6).","rationale":"I consider this the most load-bearing concern because the central claim—that ε_a can be estimated better at high bias—has no stable meaning unless the probability distribution in Eq. (12) is the reduced density matrix of the donor-acceptor system (or some explicitly identified observable). The paper's own Eq. (7) identifies the reduced system as the vibrational TLS while the text identifies p1,p2 as donor/acceptor occupations. The interaction operator mismatch (σ_z vs σ_x) and the cross-lead Fermi-function products further indicate that the master equation may not follow from the stated Hamiltonian. The numerical FI curves, even if computed correctly from Eqs. (10)-(11), are therefore not evidence that the donor-acceptor energies can be estimated as claimed. I agree with the reader's weakest_assumption; the σ_z/σ_x discrepancy sharpens the same underlying ambiguity. The proposed derivation test would settle the interpretation. If the test confirms the vibrational-TLS reading, the donor/acceptor claims collapse; if it confirms the donor/acceptor reading, the model needs a corrected derivation before the numerics are trustworthy. In either case the current manuscript does not justify its physical conclusions, so I would keep the reader's REJECT verdict.","tokens_in":11190,"tokens_out":10963,"duration_ms":95648,"concrete_test":"From the total Hamiltonian (1)-(6), derive the Born-Markov-Redfield master equation explicitly, treating (a) the vibrational TLS alone as the reduced system (tracing out both electronic sites and leads) and (b) the combined donor-acceptor-vibration four-state system as the reduced system (tracing out only the leads). Compare the resulting closed population equations with Eqs. (8)-(9) and the rate definitions α±. If neither derivation reproduces Eqs. (8)-(9), the claimed Fisher-information curves are not for the claimed donor-acceptor populations. If a derivation does reproduce them, identify unambiguously which physical states p1 and p2 denote. As a secondary check, recompute the I(ε_a) curve of Fig. 1(c) with the corrected Fermi function f_R = (exp((ε_a − µ_R)/k_B T_R)+1)^{-1} instead of the printed µ_L, and record whether the bias-dependence claim survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The entire Fisher-information analysis is built on the two-population master equation (8)-(9), which is neither derived nor consistent with the model Hamiltonian. In Section 2.1 the reduced density matrix is introduced as ρ_TLS for a two-level vibrational Hamiltonian, and Eq. (7) writes the interaction as H_int = σ_z F_c(t), i.e., the vibrational TLS is the open system and the electronic leads act as a bath. If this is the correct reading, p1(t) and p2(t) are vibrational ground and excited state populations, not donor and acceptor occupations; the subsequent I(ε_a) and I(ε_d) curves would then be Fisher informations of the vibrational state with respect to these parameters, not the electronic site populations the text repeatedly claims. If instead p1 and p2 are donor/acceptor electronic populations, Eq. (7) does not describe that system: the physical electron-vibration coupling in Eq. (6) is κ(c†_d c_a + c†_a c_d) σ_x, not σ_z F_c(t), and a first-order Born-Markov treatment would not produce rates such as α+_da = f_L[1−f+_R] that multiply left- and right-lead Fermi functions; such cross-lead products suggest a higher-order cotunneling process or a different reduced system. The derivation is relegated to Refs. [26,27,30] without reproducing the steps. Additional signs of unreliability include the definition of f_R in Eq. (9) erroneously using µ_L instead of µ_R, and Section 4 dismissing Franck-Condon factors with the sentence 'since the space is spanned by a single photon.' Because the physical interpretation of every Fisher-information curve depends on the identity of the states in Eq. (12), this unresolved ambiguity is the most load-bearing weakness.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies Fisher information for parameter estimation in a nonequilibrium donor-acceptor rectifier coupled to an anharmonic vibrational mode. The authors adopt a two-level Pauli master equation for donor/acceptor populations (Eqs. 8–9) from Refs. [26,27,30], compute time-dependent classical Fisher information for the donor energy ε_d, acceptor energy ε_a, and vibrational frequency ω0, and report an optimal measurement time for ε_a and ε_d but steady-state saturation for ω0. They then introduce a five-level model that includes vibrational states |0⟩ and |1⟩ and repeat the steady-state Fisher information calculation for ω0. The central claims are that ε_a is the most precisely estimable parameter, especially at high bias, and that ω0 is best estimated at low temperatures and long times.","tokens_in":11627,"tokens_out":12692,"duration_ms":110636,"significance":"If the results are correct, the paper provides a concrete example of parameter-specific Fisher information profiles in a molecular transport model, extending prior work on optimal time in open quantum systems to a vibrationally coupled rectifier. The distinction between parameters best estimated at finite time and those requiring steady state could inform metrological protocols for molecular junctions. However, the significance is currently limited because the master equation is not derived in the manuscript and the text is ambiguous about which physical degrees of freedom form the reduced system; without resolving these points, the Fisher information curves cannot be uniquely interpreted. The paper does present the classical Fisher information formula for diagonal states clearly and performs a systematic parameter scan, which are positive features.","major_comments":[{"comment":"The master equation used for all Fisher information calculations is not derived; the authors simply refer to Refs. [26,27,30]. More seriously, the text introducing Eq. (7) is internally inconsistent: the reduced density matrix ρ_TLS is described as a two-level vibrational system (H_vib = ω0/2 σ_z, H_int = σ_z F_c(t)), but p1(t) and p2(t) are immediately called donor and acceptor occupation probabilities. The electron-vibration coupling in Eq. (6) is κ(c†_d c_a + c†_a c_d)σ_x, which is not of the form σ_z F_c(t). If ρ_TLS is meant to describe the donor/acceptor charge state, Eq. (7) should involve the tunneling Hamiltonian H_c of Eq. (4) rather than σ_z F_c(t). This ambiguity makes the physical interpretation of the Fisher information curves in Figs. 1–3 unclear. The authors must either derive the master equation from the model Hamiltonian or provide a precise mapping to the cited works, including the explicit form of the interaction Hamiltonian used.","section":"Section 2.1, Eqs. (7)–(9)"},{"comment":"In the definitions after Eq. (9), f_R is written as (exp((ε_a − μ_L)/(k_B T_R)) + 1)^{-1}, with the left chemical potential μ_L appearing in the right-lead Fermi function. This is thermodynamically inconsistent and changes the rates α. The authors should correct this to μ_R and explicitly state whether the numerical results were computed with the corrected expression. In addition, the rates such as α+_da = f_L[1 − f+_R] multiply Fermi functions from different leads without an accompanying derivation; since this product form is characteristic of a sequential or cotunneling process, the authors should explain how it arises from the Hamiltonian (1)–(6).","section":"Eq. (9), Fermi function definitions"},{"comment":"The five-level model neglects Franck-Condon factors with the justification 'since the space is spanned by a single photon.' This is incorrect: the vibrational space consists of two phonon states, and the factors F_n = |⟨n|e^{λ/ω0(b†−b)}|0⟩|^2 are precisely the quantities that encode the electron-vibration coupling strength. Neglecting them without a small-λ argument removes the vibrational coupling from the rates, so the five-level calculation no longer represents the anharmonic rectifier introduced in Section 2. The additional use of δ_n0 in the acceptor-lead rates (Eq. (15)) eliminates vibrationally assisted acceptor tunneling altogether. Please include the Franck-Condon factors or provide a physically valid justification for their neglect.","section":"Section 4, Eq. (15) and text after Eq. (16)"},{"comment":"The paper's central claim that 'ε_a can be estimated better if the bias is high' and that I(ε_a) ≫ I(ε_d) is based on comparing curves across Figs. 1(c,d) and 2(a–d) that are plotted in arbitrary units. The authors should provide normalized plots or numerical values to support the magnitude comparison. This is secondary to the previous comments, but it is necessary to make the main conclusion quantitative.","section":"Section 3, Figs. 1–2 and Section 5"}],"minor_comments":[{"comment":"The abstract states that the Fisher information is calculated 'by deriving a quantum master equation,' but the master equation is not derived in the paper; it is taken from Refs. [26,27,30]. Please align the abstract with the main text.","section":"Abstract and Section 2.1"},{"comment":"The paper consistently misspells 'Fisher' as 'Fischer'; please correct throughout.","section":"Throughout"},{"comment":"In the caption of Fig. 1b, 'p1(1)' should read 'p1(t)'.","section":"Fig. 1b caption"},{"comment":"The sentence 'Since, the model is well accepted and studies, for a complete derivation ... we simply refer to the following works' is ungrammatical and should be rewritten.","section":"Section 2.1"},{"comment":"The phrase 'This can confirmed from the peak positions' should be 'This can be confirmed from the peak positions'.","section":"Section 3"},{"comment":"The notation 'Frank-Condon' should be 'Franck-Condon', and 'single photon' should be 'single phonon' (though the sentence remains incorrect as noted in the major comments).","section":"Section 4"},{"comment":"The figures are not always referenced in the order they appear; for example, Fig. 2 is discussed before Fig. 3 in the text, and some parameter values appear only in captions. A table of parameters for each figure would improve readability.","section":"Figures and captions"}],"recommendation":"major_revision","confidential_remarks":"The paper has a potentially useful idea but is not publishable in its current form because the master equation is underived and the text is internally inconsistent about the reduced system. A thorough revision that derives or precisely references the master equation, fixes the typo in f_R, and properly handles Franck-Condon factors is required. I would be willing to review a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper reports a clear and potentially interesting numerical observation—I(ω0) saturates rather than peaking before steady state, while I(ε_a) and I(ε_d) peak at an optimal time. I think that contrast is real, given the equations in the paper. But the equations themselves are the problem. The master equation is never derived, and the paper is internally ambiguous about whether p1, p2 are vibrational populations or donor/acceptor occupations. Eq. (7) treats the vibrational TLS as the open system with H_int = σ_z F_c(t), while the text calls p1, p2 donor and acceptor probabilities. The rates in Eq. (9) multiply left- and right-lead Fermi functions, which does not look like a standard Born–Markov–Redfield result for that interaction, and the derivation is just cited to Refs. [26–28]. The f_R function in Eq. (9) uses μ_L instead of μ_R, which is probably a typo, but it's not the only sign of carelessness: Section 4 justifies dropping Franck–Condon factors by saying 'the space is spanned by a single photon,' which is nonsense. The donor/acceptor comparison is also confounded by using asymmetric lead temperatures (T_L = 2, T_R = 1) in the key figures.\n\nWhat is genuinely new: the saturation of I(ω0) is not in earlier work, and the systematic scan over bias and temperature is a reasonable first exploration. The anthracene-PMDA parametrization grounds the numbers. I would not dismiss the numerical work as faked; the FI curves probably follow from the stated equations. But because the physical identity of the states is ambiguous, every curve loses its interpretation, and the central claims about estimating donor/acceptor energies are unsupported.\n\nMy recommendation: this is not ready for peer review. A serious editor should desk reject unless the authors fix the model description: derive or properly invoke the master equation for the actual reduced system, clarify what p1 and p2 are, correct the f_R typo, and replace the 'single photon' statement with a real justification. If that happens, the saturation observation might be worth a second look. For now, I would not cite it.","headline":"Ambiguity about whether the reduced system is the vibration or the electronics sinks the Fisher-information claims, despite a genuinely new saturation observation.","tokens_in":12203,"tokens_out":4107,"would_cite":false,"duration_ms":36101,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper finds that in a donor–acceptor rectifier model, the acceptor energy is the most precisely estimable parameter, best measured at high bias before the system reaches steady state, while the vibrational frequency is best estimated…","keywords":["quantum Fisher information","donor-acceptor rectifier","vibrationally assisted transport","parameter estimation","nonequilibrium steady state","anharmonic vibrational mode","molecular junction metrology"],"falsifier":"Measure the donor and acceptor populations in an anthracene–PMDA-style junction as a function of time at several biases and temperatures and compare the empirical Fisher information (from repeated population estimates) with the predicted I(ε_a) and I(ω_0) curves. A direct check: at high bias, I(ε_a) should peak at an earlier time than at low bias, and I(ω_0) should keep rising to steady state rather than peaking—if the peak times do not shift with bias or if I(ω_0) shows a maximum before steady state, the central claims fail.","tokens_in":10982,"feed_emoji":"⚛","tokens_out":7264,"duration_ms":63889,"temperature":0.7,"pith_summary":"This paper asks which parameters of a nonequilibrium donor–acceptor molecular junction can be estimated most precisely from measurements of its two-site occupation probabilities. Using the classical Fisher information of the time-dependent populations, it reports three parameter-specific behaviors: the acceptor energy shows a sharp maximum at an optimal interrogation time that shifts earlier as the bias grows; the donor energy is weakly estimable and barely affected by vibrations; and the anharmonic vibrational frequency shows no optimal time—its Fisher information instead saturates at steady state and grows at low temperature and low frequency. The authors conclude that acceptor energy is the best-estimated parameter, especially under high bias and strong vibronic coupling, which matters because it suggests which experimental settings (bias, timing, temperature) would most effectively reveal each molecular-junction parameter.","feed_headline":"High bias makes acceptor energy the most measurable junction parameter","feed_subtitle":"Acceptor energy peaks before steady state; vibrational frequency is best read in steady state at low temperature.","key_machinery":"The central object is the classical Fisher information I(θ)=Σ_i (1/p_i)(∂p_i/∂θ)^2 built from the time-dependent donor and acceptor populations p_1(t), p_2(t) of a two-level (spin-fermion) master equation, with rates determined by Fermi functions of the leads and by vibrational excitation or de-excitation. Because the reduced density matrix is diagonal in the site basis, the quantum Fisher information coincides with the classical Fisher information, so population measurements are optimal for parameter estimation. The paper also uses a five-state Liouville-space extension (empty state, donor or acceptor with ground or excited vibration) to check the steady-state behavior of I(ω_0). This machinery converts parameter sensitivity into a time-dependent scalar whose peak-or-plateau shape decides when and under what bias or temperature a parameter is best measured.","core_discovery":"On the paper's own terms: for a junction in which electron transfer between donor and acceptor is mediated entirely by an anharmonic vibrational mode, the classical Fisher information computed from the time-dependent populations of the reduced master equation yields three parameter-specific estimation regimes. I(ε_a) peaks at a finite optimal time before steady state, grows with vibrational frequency and bias, and its peak shifts to earlier times as bias increases; I(ε_d) is consistently smaller and insensitive to the vibrational mode; I(ω_0) does not peak—it saturates monotonically toward the steady-state value, and its magnitude increases as temperature drops and as the mode frequency decreases. A five-level extension of the model confirms that steady-state I(ω_0) rises at low temperature and decays nonlinearly to zero as ω_0 grows. The paper's claim is that acceptor energy is the most precisely estimable system parameter, especially under high bias and strong vibronic coupling.","pith_inferences":["An implicit consequence of the parameter-specific time profiles is that a single measurement-time choice cannot optimise all three parameters: protocols tuned for acceptor-energy readout (early times, high bias) will be poor for vibrational-frequency readout, which requires steady state.","The resonance structure visible in the I(ω_0) contours—maxima near ε_a + ω_0 = μ_R—suggests a testable strategy: tuning the right-lead chemical potential to this resonance should enhance vibrational-frequency estimation, a prediction the paper notes in the contours but does not elevate to a standalone claim.","The same classical Fisher formalism could be extended to estimate the electron–vibration coupling strength κ or lead temperatures, though the diagonal-density-matrix assumption would need re-checking if coherences were included, since then classical and quantum Fisher information would no longer coincide."],"forward_implications":["If the model is right, time-resolved transport or spectroscopic measurements on a donor–acceptor junction under high bias and with high-frequency vibrations will determine the acceptor energy with the highest precision, but the measurement must be taken before the system relaxes to steady state.","Lowering the temperature improves the precision of vibrational-frequency estimation, with smaller vibrational energies giving larger Fisher information, so low-temperature steady-state measurements are the right protocol for reading out the mode frequency.","Donor-energy estimation is intrinsically harder: I(ε_d) is orders of magnitude smaller than I(ε_a) and does not respond to the vibrational mode, so donor energies would require longer integration times or complementary techniques.","The optimal time for measuring acceptor and donor energies shifts earlier as the bias window is widened, meaning high-bias protocols gain precision only if measurements are fast.","In the five-level model, the steady-state information about the vibrational frequency decays to zero for very large ω_0 because the mode cannot participate in electron transfer, so estimation fails for modes too energetic to be accessed.","The saturation of I(ω_0) rather than a peak means that, unlike site energies, the mode frequency should be read after the system equilibrates, reversing the usual 'measure before steady state' rule for this parameter."],"supporting_citations":[{"why":"Supplies the earlier observation that system-parameter Fisher information is maximal before steady state, which the paper extends to site energies and contrasts with the vibrational mode's steady-state saturation.","marker":"[16]"},{"why":"Provides the spin-fermion mapping and the two-level master equation (Eqs. 8-9) whose populations the entire Fisher-information analysis is built on.","marker":"[26]"},{"why":"Justifies including vibrational anharmonicity in molecular conduction and underlies the master-equation rates used for the donor-acceptor-vibration model.","marker":"[27]"},{"why":"Gives the full counting statistics and rate structure for vibrationally assisted conduction, supporting the reduced dynamics and rate definitions.","marker":"[30]"},{"why":"Supplies the anthracene-PMDA donor/acceptor energies and the three vibrational frequencies used as the numerical prototype for the Fisher-information curves.","marker":"[31]"},{"why":"Motivates treating the localized vibrational mode as an anharmonic few-level system rather than a harmonic oscillator.","marker":"[12]"}],"fun_headline_variants":["Acceptor-energy precision peaks before steady state","Vibrational mode best read in steady state at low temperature","High bias sharpens acceptor-energy measurement window","Optimal time exists for donor and acceptor energy estimation","Acceptor energy most precise; donor insensitive"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results depend on the master equation in Eqs. (8)–(9) (cited to earlier works) being the correct reduced dynamics for this junction; if the two tracked populations are not the donor and acceptor occupations, or if the Born–Markov approximation misses essential coherences, the Fisher information curves no longer describe what the paper says they do.","fun_headline_variants_meta":{"raw":{"variants":["Acceptor-energy precision peaks before steady state","Vibrational mode best read in steady state at low temperature","High bias sharpens acceptor-energy measurement window","Optimal time exists for donor and acceptor energy estimation","Acceptor energy most precise; donor insensitive"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000919,"raw_usage":{"total_tokens":3911,"prompt_tokens":881,"completion_tokens":3030,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":497,"completion_tokens_details":{"reasoning_tokens":2958}},"tokens_in":497,"tokens_out":3030,"duration_ms":22495,"temperature":1.0,"reasoning_tokens":2958,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:42:57.710066+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the donor and acceptor populations in an anthracene–PMDA-style junction as a function of time at several biases and temperatures and compare the empirical Fisher information (from repeated population estimates) with the predicted I(ε_a) and I(ω_0) curves. A direct check: at high bias, I(ε_a) should peak at an earlier time than at low bias, and I(ω_0) should keep rising to steady state rather than peaking—if the peak times do not shift with bias or if I(ω_0) shows a maximum before steady state, the central claims fail.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the earlier observation that system-parameter Fisher information is maximal before steady state, which the paper extends to site energies and contrasts with the vibrational mode's steady-state saturation."},{"cited_title":"Simine, D","cited_arxiv_id":null,"evidence_quote":"Provides the spin-fermion mapping and the two-level master equation (Eqs. 8-9) whose populations the entire Fisher-information analysis is built on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Justifies including vibrational anharmonicity in molecular conduction and underlies the master-equation rates used for the donor-acceptor-vibration model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the full counting statistics and rate structure for vibrationally assisted conduction, supporting the reduced dynamics and rate definitions."},{"cited_title":"Fonari, N","cited_arxiv_id":null,"evidence_quote":"Supplies the anthracene-PMDA donor/acceptor energies and the three vibrational frequencies used as the numerical prototype for the Fisher-information curves."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Motivates treating the localized vibrational mode as an anharmonic few-level system rather than a harmonic oscillator."}],"review_version":1}