{"id":"07e11512-1b55-442f-a39a-fc62aa8bf6f7","arxiv_id":"2511.08689","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A trapped-ion protocol that balances controlled heating and cooling realizes thermal baths with independently tunable temperature and dissipation rate, demonstrated in spin-boson charge- and exciton-transfer simulations.","lead":"This paper demonstrates a trapped-ion method to create thermal reservoirs whose temperature and decay rate can be tuned independently, by balancing laser cooling with electric-field noise heating. The authors benchmark it on charge-transfer and exciton-transfer simulations, observing temperature-dependent transfer rates and thermally activated pathways.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Direct mode-resolved thermometry for non-COM modes is missing; two-mode local-temperature claim rests on uncalibrated mode selectivity of broadcast RF noise.","rationale":"The reader's weakest assumption is the same one I identify: the independence and mode-selectivity of the two engineered dissipators away from the single-COM demonstration. I do not see an internal inconsistency in the Lindblad derivation or in the data for the COM mode; the reported steady states and dynamics are consistent with Eq. (2), and the high-n fitting caveat is explicit. However, the paper's own supplement concedes that non-COM heating is empirical and not characterized in detail. Since the headline applications are two-mode local-temperature simulations, a direct mode-resolved thermometry check is necessary to fully support the claim. If that check passes, ACCEPT is warranted; until then, CONDITIONAL is the appropriate verdict.","tokens_in":21014,"tokens_out":12824,"duration_ms":144737,"concrete_test":"On the same two-ion chain, apply the exact preparation used for Fig. 4 (sympathetic cooling plus noise tone at ω2 with no LVC drive), then measure the phonon-number distribution of each axial mode independently via blue-sideband Rabi flopping (Eq. 3). Check (i) the ω2 steady state is thermal with nbar equal to γ_h/γ_c, and (ii) the ω1 steady state remains at its base temperature (e.g., Δnbar<0.05). If either fails, the local-temperature interpretation of Fig. 4 must be re-evaluated by repeating the experiment with direct mode-resolved calibration.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central relation Eq. (1) requires the cooling and heating channels to act as independent Lindblad dissipators with rates γ_c and γ_h. This is verified for a single COM mode in Fig. 2, but the spin-boson demonstrations use an out-of-phase mode at 3.776 MHz. The Supplemental section 'Controlled heating of non-COM motional mode' states that near-field antenna noise couples predominantly to COM and that heating non-COM modes requires roughly 20× larger amplitude and an empirically measured antenna profile. No direct blue-sideband thermometry of the non-COM modes is reported under the simultaneous cooling+heating used in Figs. 3-4. If the ω2 heating tone also drives the ω1 mode, or if the sympathetic cooling and broadcast noise do not combine as independent dissipators on the shared mode, the quoted local temperatures (nbar1,nbar2)=(0.10,0.80) and dissipation rates would be nominal, and the thermally activated interference pathway in Fig. 4 would not be a verified local-temperature effect.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes and experimentally demonstrates a scheme to engineer thermal reservoirs for trapped-ion motional modes by combining continuous resolved-sideband cooling at rate γ_c with broadcast RF electric-field noise at rate γ_h. For a single mode, the steady-state mean occupation is n_ss=γ_h/γ_c and the relaxation rate is γ_c, enabling independent control of temperature and dissipation. On a COM mode of a single 171Yb+ ion, the authors verify exponential phonon dynamics consistent with Eq. (2) for multiple initial temperatures and heating rates, and find steady-state phonon distributions consistent with thermal states. They then apply the method to a dual-species chain to simulate finite-temperature charge transfer and two-mode vibrationally assisted exciton transfer, observing a broadening of the charge-transfer rate spectrum at higher temperature and a thermally activated mixed-mode resonance in the two-mode system.","tokens_in":21277,"tokens_out":9272,"duration_ms":87819,"significance":"If the results hold, this provides a practical and versatile tool for trapped-ion quantum simulation of open systems at finite temperature, extending earlier cooling-only or noise-only reservoir-engineering methods. The COM-mode benchmark is convincing: independent γ_h and γ_c control, dynamics following Eq. (2), and thermal steady-state distributions are demonstrated with quantitative agreement. The applications to charge-transfer and two-mode exciton-transfer models illustrate the method's relevance to chemical dynamics, and the temperature-dependent transfer-rate spectra in Fig. 3 are consistent with numerical simulations. The paper also includes a self-contained analytic derivation of the thermal steady state and an alternative all-laser protocol, though the latter is not experimentally demonstrated. The main weakness is the lack of direct thermometry for the non-COM modes used in the simulations, which tempers the strength of the local-temperature claims.","major_comments":[{"comment":"The central reservoir demonstration (Fig. 2) is performed on a COM mode of a single ion, with direct phonon-number thermometry. However, the finite-temperature LVC simulations in Figs. 3-4 are performed on the out-of-phase mode at 3.776 MHz (and the ω1,ω2 modes in Fig. 4). The Supplemental states that broadcast RF noise couples predominantly to the COM mode and that heating non-COM modes requires ~20× amplitude and an empirically measured antenna profile, but no direct blue-sideband thermometry is reported for these modes under simultaneous cooling+heating. The values (nbar1,nbar2)=(0.10,0.80) in Fig. 4 and nbar=0.15,0.80 in Fig. 3 thus appear to be nominal settings rather than independently calibrated temperatures. Since the theory curves are numerical simulations of the same master-equation model, agreement with data tests self-consistency of the model, not the local-temperature claim.","section":"Supplemental 'Controlled heating of non-COM motional mode'; Figs. 3-4 and Eq. (6)"}],"minor_comments":[{"comment":"The dissipator is introduced as D_c[ρ] but then used as D_a and D_{a†}; please clarify the notation by defining D_c for a general operator c and then applying it to a and a†.","section":"Eq. (1)"},{"comment":"The statement that a single-phonon excitation rate is 'proportional to p_ni(ni+1)' is ambiguous: the rates for absorption and emission involve p_ni(n_i+1) and p_{ni+1}(n_i+1) respectively. Please specify the two cases.","section":"End Matter, 'Finite-temperature vibrationally assisted exciton transfer'"},{"comment":"Unlike Fig. 3A, no error bars are mentioned for the experimental data points in Fig. 4A. If error bars were omitted for clarity, state this and point to the Supplemental; if not available, the significance of the thermally activated resonance near 0.30ω1 is hard to assess.","section":"Fig. 4 caption"},{"comment":"The dephasing rates γ_z and γ_m are extracted by comparing simulations to the same experimental data used in Figs. 3-4. Please state whether the engineered nbar and γ values were held fixed during this extraction, and provide a brief sensitivity analysis or uncertainty estimate for the plotted theoretical curves.","section":"Supplemental, 'Numerical calculations for finite-temperature excitation transfer'"}],"recommendation":"major_revision","confidential_remarks":"The COM-mode demonstration is solid and the method is clearly useful. The main gap is the absence of direct non-COM thermometry; if the authors can supply those calibration data or clearly mark the nbar values as fitted, I would be happy to support acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: the central claim holds up, and the paper deserves a serious referee. The group combines continuous sympathetic/resolved-sideband cooling with broadcast RF noise to give a motional mode a thermal steady state with n_ss = γ_h/γ_c, tunable independently of the dissipation rate γ_c. The Lindblad equation itself is textbook, but the experimental realization is careful: phonon-number dynamics from two initial temperatures follow Eq. (2), steady-state populations match thermal states, and n_ss scales with γ_h/γ_c across two cooling rates.\n\nThe applications are the real substance. In the strongly adiabatic charge-transfer model, raising n̄ from 0.15 to 0.80 broadens the transfer-rate spectrum, and the two-mode experiment — a mixed-mode resonance near 0.30ω1 that appears only when the ω2 mode is hot — is a genuinely nice observation. The interference-pathway explanation in the End Matter is clear.\n\nThe soft spot is the one the stress-test flagged, and it is real but not fatal. All direct thermometry in Fig. 2 is on a COM mode of a single 171Yb+ ion. The spin-boson experiments run on a non-COM mode at 3.776 MHz in a two-species chain, where the Supplemental concedes RF noise couples about 20× weaker and requires an empirically measured antenna profile. No direct blue-sideband characterization of that mode under simultaneous cooling+heating is shown, so the n̄ = 0.15 and 0.80 values in Fig. 3, and especially the local temperatures (0.10, 0.80) in Fig. 4, are nominal values from calibrated rates rather than directly measured ones. The theory-data agreement mitigates this: the mixed-mode resonance requires ω2 excited-state population specifically, so its appearance in the hot case is decent evidence the ω2 mode really is hotter. But a referee should ask for direct mode-resolved thermometry of the non-COM modes, or at least a sensitivity analysis over (n̄1, n̄2). This is a disclosed characterization gap, not a demonstrated error.\n\nMinor items: the theory curves use two small fitted dephasing rates (disclosed, small, fine), no code or data are released, and the paper honestly reports that free-parameter phonon fits break down above n̄ ≈ 3. The citation pattern is fair — Refs. [40, 41] are their own prior work, and temperature control is a genuine extension of those.\n\nWho this is for: anyone doing trapped-ion quantum simulation of open systems or reservoir engineering. It should go to peer review; the experiments look carefully done and the requested fix is modest.","headline":"Solid experimental step: independent temperature and dissipation control for trapped-ion reservoirs, with the main gap being missing direct thermometry on the non-COM modes used in the applications.","tokens_in":21716,"tokens_out":5948,"would_cite":true,"duration_ms":61818,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"By simultaneously heating and cooling a trapped ion's vibration, this paper creates a thermal reservoir whose temperature and dissipation rate are set independently, and uses it to reveal how temperature reshapes charge- and exciton-transfe","keywords":["trapped ions","thermal reservoir engineering","spin-boson simulation","charge transfer","exciton transfer","motional modes","Lindblad master equation","tunable temperature"],"falsifier":"Measure the steady-state phonon distribution of the out-of-phase mode for n_ss ≈ 5–10 using an independent thermometer, and check the detailed-balance ratio p_n/p_{n−1} = n_ss/(n_ss+1) for all n, while also monitoring neighboring modes; a deviation from the geometric distribution, or heating of spectator modes, would invalidate the single-mode thermal-bath model.","tokens_in":20956,"feed_emoji":"🌡️","tokens_out":5007,"duration_ms":50851,"temperature":0.7,"pith_summary":"This paper demonstrates a way to engineer a thermal reservoir for the vibrational motion of trapped ions, with the bath's temperature and its dissipation rate set by two independent knobs. By simultaneously laser-cooling a motional mode at rate γ_c while broadcasting electric-field noise that heats it at rate γ_h, the ion's steady state becomes a thermal state whose mean phonon number is n_ss = γ_h/γ_c; the cooling rate alone sets how fast the system equilibrates. The authors verify the predicted exponential approach to equilibrium and the thermal (geometric) steady-state population distribution. They then apply the reservoir to simulate finite-temperature charge transfer and two-mode vibrationally assisted exciton transfer, observing that higher temperature broadens the transfer-rate spectrum and that local temperature can activate otherwise-suppressed interference pathways.","feed_headline":"Ion-trap experiment dials in a thermal bath's temperature","feed_subtitle":"Balancing laser cooling and broadcast noise gives independent control of temperature and dissipation, enabling finite-temperature quantum ch","key_machinery":"The central object is the Lindblad master equation for a harmonic oscillator subject to simultaneous cooling and heating, with jump operators a and a†. Cooling contributes a rate γ_c and heating contributes γ_h = γ_c n_ss; detailed balance gives p_n/p_{n−1} = n_ss/(n_ss + 1), the defining ratio of a thermal (geometric) state. In the experiment, the cooling channel is resolved-sideband Raman laser cooling on a 171Yb+ ion (or sympathetic cooling via a 172Yb+ ancilla for the out-of-phase mode), and the heating channel is a radio-frequency noise signal broadcast from an antenna; a blue-sideband probe extracts the full phonon distribution.","core_discovery":"The paper's central claim is that a trapped-ion motional mode coupled to two controlled dissipative channels—resolved-sideband laser cooling and broadcast electric-field noise—relaxes to a true thermal state, not merely a heated or squeezed one, with mean occupation n_ss = γ_h/γ_c. Because γ_c and γ_h can be tuned independently, the temperature and the equilibration (dissipation) rate are independently controllable. The authors prove this by deriving the Lindblad master equation and its detailed-balance steady state, and they confirm it experimentally: the phonon number follows ⟨n(t)⟩ = n_ss + (n_0 − n_ss)e^{−γ_c t} from two different initial temperatures, and the steady-state populations ma","pith_inferences":["The relation n_ss = γ_h/γ_c ties temperature to the ratio of two rates; a natural next step is to sweep γ_c and γ_h together along curves of constant n_ss to measure how charge-transfer rates depend on dissipation alone while keeping temperature fixed, separating friction effects from thermal effects.","Because the heating channel is a broadcast classical noise field, the method may be transferable to other platforms with laser-coolable harmonic oscillators, such as trapped neutral atoms or optomechanical systems, provided the noise can be shaped to the target mode.","The reported difficulty of measuring steady states above n_ss ≈ 3 suggests a practical test: with a more sensitive phonon thermometer (e.g., using multiple sideband probes or ancilla-based measurements), one could check whether the thermal distribution holds at higher temperatures where free-parameter fits failed.","The thermally activated mixed-mode resonance points toward a design principle: local temperature differences between vibrational modes can route excitation along specific interference-enabled channels, which could be exploited to control energy flow in multi-mode simulators."],"forward_implications":["Finite-temperature open-system dynamics (thermal-state preparation, heat flow, thermal entanglement) can now be studied in trapped-ion simulators with independent control of temperature and dissipation.","For the charge-transfer model, higher bath temperature broadens the transfer-rate spectrum: rates drop at small donor-acceptor gaps and increase at large gaps, a temperature dependence beyond the earlier dissipation-only picture.","In the two-mode exciton-transfer model, raising the local temperature of one vibrational mode activates a mixed-mode coherent pathway that is otherwise suppressed, demonstrating thermally activated interference.","With multiple addressed modes, each mode can be given its own temperature, allowing local-temperature gradients and multi-bath spin-boson models to be realized.","Resonantly driving higher-order sideband processes with noise and cooling tones could extend the scheme to nonlinear and structured (non-Gaussian) reservoirs."],"fun_headline_variants":["Trapped-ion bath temperature dialed independently of dissipation","Thermal baths with tunable temperature for trapped-ion simulators","Ion trap realizes thermal reservoir with adjustable temperature","Independent thermal control in trapped-ion spin-boson simulator","Ion-trap experiment offers independent temperature and dissipation control"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the broadcast electric-field noise acts as an independent heating Lindblad channel on the targeted motional mode alone, adding linearly to the laser cooling; for the non-center-of-mass modes actually used in the simulations, this required empirical calibration of the near-field antenna profile and roughly a factor of 20 in extra noise amplitude, and if bystander modes heat up or the two drives do not act independently, the engineered temperatu","fun_headline_variants_meta":{"raw":{"variants":["Trapped-ion bath temperature dialed independently of dissipation","Thermal baths with tunable temperature for trapped-ion simulators","Ion trap realizes thermal reservoir with adjustable temperature","Independent thermal control in trapped-ion spin-boson simulator","Ion-trap experiment offers independent temperature and dissipation control"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00024,"raw_usage":{"total_tokens":1320,"prompt_tokens":673,"completion_tokens":647,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":417,"completion_tokens_details":{"reasoning_tokens":568}},"tokens_in":417,"tokens_out":647,"duration_ms":7440,"temperature":1.0,"reasoning_tokens":568,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T22:45:59.749108+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the steady-state phonon distribution of the out-of-phase mode for n_ss ≈ 5–10 using an independent thermometer, and check the detailed-balance ratio p_n/p_{n−1} = n_ss/(n_ss+1) for all n, while also monitoring neighboring modes; a deviation from the geometric distribution, or heating of spectator modes, would invalidate the single-mode thermal-bath model.","supporting_citations":[],"review_version":1}