{"id":"f048127f-ce20-482d-8e25-5680dfe822cf","arxiv_id":"1908.08373","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Bacterial photosynthesis can stay highly efficient without quantum coherence, and random chromophore arrangements work as well as natural rings once nearest-neighbor spacing is right.","lead":"A detailed simulation of purple bacteria light-hvesting shows that high efficiency does not require long-lived quantum coherence, and that random arrangements of chromophores can be just as efficient as the structured rings found in nature when the molecules are close enough. The results suggest simpler, cheaper designs for artificial light-harvesting devices.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 0.1 ps Hamiltonian-resampling protocol for slow vibrations is ad hoc and may be too fast; the mixed-state and structure-insensitivity conclusions depend on it without a sensitivity check.","rationale":"The reader's weakest-assumption analysis identifies precisely the same load-bearing point: the 0.1 ps random-Hamiltonian switching protocol is a hand-chosen approximation with no sensitivity analysis, and the paper's conclusions about high mixing, lack of long-range coherence, and structural insensitivity all depend on it. My reading of the full text reinforces this concern rather than replacing it with a different one. The internal inconsistency between the stated 'slow vibrations' and the 0.1 ps redraw interval (which is far shorter than the J_RM low-frequency periods and the sink timescale) makes the concern concrete and testable. I do not think this rises to rejection: the authors acknowledge the difficulty of modeling the vibrational environment, present two contrasting bath models, and note qualitative agreement with preliminary HEOM calculations. The qualitative picture of a highly mixed state under natural incoherent illumination is consistent with prior literature. However, the specific mechanistic claim that slow vibrations keep exciton density moving, and the quantitative efficiency values, rest on an unvalidated timescale. A sensitivity scan or a comparison with a smoother correlated-disorder process should be required before the central quantitative claims can be accepted as robust. Thus the appropriate verdict remains CONDITIONAL, matching the reader's assessment; no verdict change is warranted beyond that.","tokens_in":10976,"tokens_out":5351,"duration_ms":58489,"concrete_test":"Repeat the full LHI-LHII antenna and random-chromophore-box simulations with switching interval τ_sw ∈ {0.01, 0.05, 0.1, 0.5, 1, 5, 20 ps}, or replace sudden redraws with a smooth Gaussian stochastic process whose correlation time matches the low-frequency peak of J_RM (~20–60 ps), keeping all other rates fixed. If the time-averaged efficiency changes by more than 10% or trρ²/IPR changes qualitatively, the mixed-state and structure-insensitivity conclusions are artifacts of the 0.1 ps redraw rate. As a complementary check, compare the resulting dynamics against exact HEOM on the small benchmark system referenced in ref 42.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is the treatment of slow vibrational modes by redrawing the entire random Hamiltonian every 0.1 ps (Section 'Coupling to the environment', second paragraph). The paper's central mechanistic claims—that slow vibrations keep exciton density moving, that the state remains permanently highly mixed with no long-range coherence, and that Anderson localization at low dephasing is mitigated—are all generated by this time-dependent disorder. The timescale is not independently justified: the Renger-Marcus spectral density J_RM used for the thermal bath (Eq. 7) peaks at characteristic periods 2π/ω1 ≈ 60 ps and 2π/ω2 ≈ 17 ps, an order of magnitude longer than the 0.1 ps switching interval, and the sink timescale is ~8 ps. Sudden, independent Hamiltonian redraws every 0.1 ps therefore resemble fast white-noise fluctuations rather than the slowly varying static disorder the text describes. If physically motivated slow fluctuations occur on tens-of-picosecond timescales, the system would spend most of its time near a single thermal equilibrium, and the quantitative values of purity, IPR, and efficiency could shift substantially. No sensitivity scan over the switching interval or over a correlation-time parameter is presented, leaving the central claims exposed to this model choice.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a master-equation model of the bacterial photosynthetic antenna consisting of one LHII complex, one LHI complex, and the reaction centre (65 chromophores), including explicit weak incoherent sunlight absorption, radiative and non-radiative losses, a reaction-centre sink, and two alternative models of the vibrational bath (a local dephasing model and a global thermalizing Redfield model). Slow vibrational modes are incorporated by redrawing a random Hamiltonian every 0.1 ps, with energetic and dipole disorder taken from earlier TDDFT/MD work. The authors report a time-averaged efficiency of 96% with the global thermalizing model and 62% with the local dephasing model, and they find that the one-exciton state remains highly mixed and largely incoherent. In a second set of simulations, randomly arranged chromophores in a box around a reaction centre also give high efficiency at sufficiently high concentration, which the authors interpret as evidence that detailed structural arrangement is not important for efficiency.","tokens_in":11309,"tokens_out":4587,"duration_ms":49415,"significance":"If the central conclusions are correct, the paper makes a substantive contribution: it argues that high photosynthetic efficiency does not require long-range quantum coherence and that random dense packing of chromophores may be sufficient, with implications for artificial light-harvesting design. The model is unusually comprehensive in its simultaneous treatment of sunlight, loss channels, the sink, and two opposing bath descriptions, and the inclusion of TDDFT/MD disorder and a large combined LHII–LHI–RC system goes beyond typical single-complex studies. The random-box simulations provide a falsifiable, concentration-dependent prediction. However, the quantitative efficiency and the mechanistic claims about slow vibrations rest on a heuristic time-dependent-disorder protocol that is not validated or subjected to sensitivity analysis, and the abstract's efficiency claim is stronger than the two-bath-model results support. The paper is therefore promising but requires substantial revision.","major_comments":[{"comment":"The 0.1 ps Hamiltonian-resampling protocol for slow vibrational modes is not independently validated, and no sensitivity scan over the switching interval is presented; the characteristic periods of the Renger–Marcus spectral density in Eq. (7) are 2π/ω1 ≈ 60 ps and 2π/ω2 ≈ 17 ps, an order of magnitude longer than the 0.1 ps interval, so the protocol may be closer to rapid white-noise-like fluctuations than to the slowly varying static disorder described in the text. The conclusions that the state remains permanently highly mixed, that the IPR is small, and that disorder-driven localization is mitigated all depend on this choice, and a range of switching times (e.g., 0.01–10 ps) should be tested before these claims can be considered robust.","section":"Coupling to the environment, second paragraph; Figs. 1-2; Table II"},{"comment":"The abstract's statement that 'our model describes the experimentally observed high efficiency of light harvesting' is contradicted by the local dephasing model, which yields a time-averaged efficiency of 62% rather than the high experimental range; the two bath models therefore give qualitatively different efficiency predictions, and the abstract should report or at least qualify this discrepancy. Furthermore, because Γsink = 0.125 ps−1 is 125 times larger than Γnr = 0.001 ps−1, the efficiency measure in Eq. (8) is heavily influenced by the fixed rate ratio, and the paper does not isolate how much of the predicted high efficiency is a consequence of these chosen kinetic constants rather than of the antenna structure or dynamics; a sensitivity analysis over Γsink and Γnr is needed.","section":"Abstract; Results, Fig. 3; Eq. (8)"},{"comment":"The global thermalizing model relies on the B777-fitted spectral density JRM being transferred to the LHII/LHI/RC system, and the only validation cited is the reproduction of an LHII relaxation time; no validation is provided for LHI or the combined antenna, and the text itself admits that the global model overestimates the RC population. Since the headline 96% efficiency comes from this model, the paper should present this value with the acknowledged bias and should quantify the uncertainty arising from the spectral-density transfer.","section":"Eqs. (5)-(7); SI Fig. S2, SI Fig. S4"},{"comment":"The conclusion that chromophore orientation is not significant is inferred from the small variance of efficiency across 32 random configurations at high concentration, but this design varies positions and orientations simultaneously; an explicit control in which positions are fixed and only orientations are randomized is required to support the claim that orientation is irrelevant, and the paper should report effect sizes or confidence intervals rather than only standard-deviation bars.","section":"An artificial light-harvesting system; Fig. 4"}],"minor_comments":[{"comment":"There are several typographical errors, including 'futher' in the Discussion, 'though to imply' in the artificial light-harvesting section, and 'maintian' in the Discussion; these should be corrected.","section":"Throughout"},{"comment":"The text states that the global Lindblad approach gives 'good qualitative agreement' with HEOM calculations, but no data from the HEOM comparison are shown; a figure or a more detailed statement of the comparison should be provided in the SI.","section":"Coupling to the environment"},{"comment":"The inverse participation ratio is defined for the full n-site system, but Table II reports IPR values for subsystems of different sizes; the normalization basis for the subsystem IPR should be stated explicitly to allow interpretation of the numbers in Table II.","section":"Results, Table II; Eq. (9)"},{"comment":"The time-averaged efficiency is stated to be calculated 'from 5 ps', but the total simulation duration for the full antenna is not given; this duration should be reported for reproducibility and to allow comparison with the 10–100 ps averaging used for the random-box simulations.","section":"Results, Fig. 3"},{"comment":"The distributions for the energetic and dipole disorder are cited to reference [12], but the functional forms and parameter values are not reproduced in the main text or the SI; providing those details would make the model reproducible.","section":"Theory, Hamiltonian disorder"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope and addresses an important question, but the central mechanistic and efficiency claims rest on the ad hoc 0.1 ps Hamiltonian-resampling protocol and on rate constants that dominate the efficiency measure. I would encourage the editor to request a sensitivity analysis of the switching time and a more careful qualification of the abstract's efficiency claim. The paper may become publishable after these revisions, but in its current form the load-bearing approximations are not sufficiently validated."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is worth a serious referee. The whole-antenna model (LHI+LHII+RC) with explicit incoherent sunlight, corrected point-dipole couplings, and both local and global bath treatments is an ambitious piece of modelling. The random-box study is the most useful new bit: it shows that spacing sets the efficiency and orientation barely matters, which is a concrete design rule for artificial antennas. The authors are also honest that the no-coherence conclusion is not theirs alone, and they openly state that both bath models have regime problems. That is good scholarly behavior.\n\nThe soft spot is load-bearing: the time-dependent disorder protocol. The Hamiltonian is redrawn every 0.1 ps, but the spectral density they use has its slow modes at 60 ps and 17 ps periods. Redrawing on 0.1 ps is therefore not slow vibrational disorder; it resembles fast white noise. The text's justification for 0.1 ps via thermal occupation of higher harmonic levels does not hold up either—the corresponding frequency is above the thermal threshold, not below it. The mixed-state and structure-insensitivity conclusions are generated by this protocol, and there is no sensitivity scan over the switching interval. This is a major concern, not a quibble. The stress test's calculation is right, and until we see how slower resampling (say 1–10 ps) changes purity, IPR, and efficiency, the central mechanism claim is exposed. It is fixable, but it needs to be fixed.\n\nAlso, the abstract overstates. \"Our model describes the experimentally observed high efficiency\"—the local model gives 62%, not the 84–91% range. The efficiency is largely predetermined by the measured sink rate being 125 times the loss rate, so the model is not deriving high efficiency from first principles. That is fine, but the language should be disciplined.\n\nThe random-box result uses only the local dephasing model, which is appropriate for low concentrations, and the concentration trends look plausible. The caveats about concentration quenching and synthetic cost are sensible.\n\nThis paper deserves peer review. I would send it out, with referees asked to push on the time-dependent disorder sensitivity. As is, it is conditional: the qualitative message (coherence not required, spacing matters) likely survives, but the quantitative efficiency claims and the mechanism story need more work.","headline":"A genuinely useful whole-antenna model and a clever random-box design study, but the time-dependent disorder protocol is under-justified and the efficiency claims are overstated.","tokens_in":11781,"tokens_out":3742,"would_cite":true,"duration_ms":38793,"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":"This paper argues that bacterial light-harvesting achieves high efficiency without long-range quantum coherence, and that random chromophore packing suffices if nearest-neighbour spacing matches natural values.","keywords":["photosynthetic light harvesting","exciton transport","quantum coherence","Lindblad master equation","time-dependent disorder","purple bacteria","artificial light-harvesting devices","efficiency"],"falsifier":"Vary the Hamiltonian-reset interval in the model (for example from 0.05 ps to 1 ps) and recompute the time-averaged efficiency and purity; if the efficiency changes by more than a few percent or significant coherence appears at longer intervals, the sudden-switch representation of slow vibrations is falsified as the explanation.","tokens_in":10807,"feed_emoji":"☀️","tokens_out":12396,"duration_ms":107617,"temperature":0.7,"pith_summary":"This paper asks what actually makes a photosynthetic antenna efficient: whether the elaborate ring-like arrangement of pigment molecules is essential, or whether dense random packing is enough. The authors build a master-equation model of a complete purple-bacterium antenna (one LHII ring, one LHI ring, and the reaction centre, 65 chromophores in all) driven by weak incoherent sunlight, with radiative and non-radiative losses, a reaction-centre sink, and a vibrational bath. Their central finding is that the one-exciton state stays highly mixed at all times, with essentially no long-range quantum coherence, yet the model reproduces the experimentally observed high efficiency (96% with a thermalizing bath, 62% with local dephasing, bracketing the measured 84–91%). They then replace the natural ring structure with random boxes of chromophores and find high efficiency whenever typical nearest-neighbour distances match natural values, concluding that detailed chromophore orientation is not what governs antenna efficiency.","feed_headline":"Bacterial photosynthesis stays efficient without quantum coherence","feed_subtitle":"Randomly packed pigments match Nature's rings when spaced alike, so precise chromophore layout is not the key.","key_machinery":"The mechanism that carries the argument is a Lindblad master equation for the one-exciton density matrix, with four dissipative pieces representing absorption and emission of sunlight, loss to the reaction-centre sink, non-radiative decay, and the vibrational bath: $\\frac{d\\rho(t)}{dt} = -\\frac{i}{\\hbar}[\\hat H(t),\\rho(t)] + (\\mathcal{L}_{\\mathrm{rad}}+\\mathcal{L}_{\\mathrm{bath}}+\\mathcal{L}_{\\mathrm{nr}}+\\mathcal{L}_{\\mathrm{sink}})\\rho(t)$. The distinctive feature is that the Hamiltonian $\\hat H(t)$ is re-randomized every 0.1 ps: site energies and dipoles are redrawn from distributions taken from density-functional/molecular-dynamics snapshots, mimicking the symmetry-breaking effect of slow vibrational modes. This time-dependent disorder, combined with a thermalizing bath term that relaxes the system to the instantaneous Boltzmann state, is what keeps the exciton state highly mixed and continuously moving; it replaces the functional role that long-lived coherence was once thought to play.","core_discovery":"The central claim is that high light-harvesting efficiency in purple bacteria does not require long-range quantum coherence, and that the detailed spatial arrangement of chromophores is not a significant determinant of that efficiency. In the model the excitonic subsystem settles into a state that is close to a mixture of the instantaneous thermal states of a randomly fluctuating Hamiltonian; the average state has an inverse participation ratio near the incoherent limit, and population is spread evenly across equivalent sites within each ring. The same efficiency emerges when chromophores are packed randomly, in three dimensions or in a slab, around a reaction centre, provided their typical nearest-neighbour separation is comparable to the natural one. What matters, the authors argue, is the density of chromophores and the balance between dephasing, decay, trapping, and thermalization, not a specific architecture.","pith_inferences":["If the random-box result extends to other antenna geometries, a cheap design rule for synthetic systems is to enforce only a minimum nearest-neighbour distance while letting chromophores orient freely, which is far more scalable than current coherence-preserving designs.","A direct validation would be to compare the sudden-switch disorder model with an explicit continuous low-frequency mode calculation on a small ring; disagreement would indicate that the 0.1 ps reset interval is doing too much work.","The same modelling approach could be applied to other pigment-protein networks, such as cryptophyte antennas or cyanobacterial phycobilisomes, to test whether density rather than architecture predicts efficiency across organisms."],"forward_implications":["Artificial light-harvesting devices do not need precise control of chromophore orientation; matching near-neighbour spacing to natural values should be enough for high efficiency.","Quantum coherence is not a prerequisite for efficient energy transport under natural illumination; the combination of thermalizing relaxation and slow vibrational disorder suffices.","The key design variables for an antenna become chromophore concentration, dephasing rate, and the balance of decay and trapping, rather than ring symmetry or molecular orientation.","The apparent importance of circular ring structures in nature may be a consequence of packing: rings achieve high local density while avoiding the concentration quenching that a fully disordered solution would suffer.","Different choices for the system-bath model (global thermalizing vs. local dephasing) change the predicted efficiency by tens of percent, so quantitative predictions remain model-dependent."],"supporting_citations":[{"why":"Supplies the site-energy and dipole disorder distributions from atomistic snapshots used in the time-dependent disorder model.","marker":"[12]"},{"why":"Provides the precedent for representing slow vibrational modes by sudden random Hamiltonian changes.","marker":"[37]"},{"why":"Gives the global thermalizing bath model used as one variant of the vibrational coupling.","marker":"[38]"},{"why":"Provides the fitted spectral density used to set thermalization rates and reproduce measured LHII relaxation times.","marker":"[44]"},{"why":"Defines the inverse participation ratio used to quantify the absence of coherence in the one-exciton state.","marker":"[48]"},{"why":"Documents the weaker exciton-transfer coupling between the B875 ring and the reaction centre, explaining why the global bath model overestimates efficiency.","marker":"[49]"},{"why":"Shows atomistically that excitons keep moving around the LHII ring, supporting the proposed transport role of slow vibrations.","marker":"[53]"}],"fun_headline_variants":["Random pigment packing matches Nature's light harvest","No quantum coherence needed for bacterial photosynthesis","Efficiency holds without precise chromophore layout","Light harvesting: density beats structure in bacteria","Random rings still catch sunlight efficiently"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusions rest on the assumption that slow molecular vibrations can be represented by generating a fresh random Hamiltonian every 0.1 ps; if real slow modes act differently on that timescale, the computed mixed state and efficiencies could change.","fun_headline_variants_meta":{"raw":{"variants":["Random pigment packing matches Nature's light harvest","No quantum coherence needed for bacterial photosynthesis","Efficiency holds without precise chromophore layout","Light harvesting: density beats structure in bacteria","Random rings still catch sunlight efficiently"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000117,"raw_usage":{"total_tokens":1063,"prompt_tokens":912,"completion_tokens":151,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":528,"completion_tokens_details":{"reasoning_tokens":88}},"tokens_in":528,"tokens_out":151,"duration_ms":2498,"temperature":1.0,"reasoning_tokens":88,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:21:51.910177+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Vary the Hamiltonian-reset interval in the model (for example from 0.05 ps to 1 ps) and recompute the time-averaged efficiency and purity; if the efficiency changes by more than a few percent or significant coherence appears at longer intervals, the sudden-switch representation of slow vibrations is falsified as the explanation.","supporting_citations":[{"cited_title":"Stross, M","cited_arxiv_id":null,"evidence_quote":"Supplies the site-energy and dipole disorder distributions from atomistic snapshots used in the time-dependent disorder model."},{"cited_title":"Ostilli and C","cited_arxiv_id":null,"evidence_quote":"Gives the global thermalizing bath model used as one variant of the vibrational coupling."},{"cited_title":"Renger and R","cited_arxiv_id":null,"evidence_quote":"Provides the fitted spectral density used to set thermalization rates and reproduce measured LHII relaxation times."},{"cited_title":"Meier, V","cited_arxiv_id":null,"evidence_quote":"Defines the inverse participation ratio used to quantify the absence of coherence in the one-exciton state."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the weaker exciton-transfer coupling between the B875 ring and the reaction centre, explaining why the global bath model overestimates efficiency."},{"cited_title":"Sisto, C","cited_arxiv_id":null,"evidence_quote":"Shows atomistically that excitons keep moving around the LHII ring, supporting the proposed transport role of slow vibrations."}],"review_version":1}