{"id":"3536763e-5134-4323-8157-ebb66a2946c8","arxiv_id":"2510.01823","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Treating the oscillator potential as a state removes terminal-jump boundary conditions from minimum-work protocols and reduces bulk optimal dynamics to a universal centre-manifold turnpike.","lead":"Proposes a reformulation of optimal control for fast thermodynamic transitions, treating the trap stiffness as a state variable so work-minimizing protocols no longer need discontinuous endpoint jumps. Gives a universal \"turnpike\" description of optimal protocols and separates minimum-dissipation equilibration from minimum-work driving to non-equilibrium states.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The no-terminal-jump conclusion is a property of the state augmentation in Eq. (2), not a theorem about stiffness-as-control formulations; it is conditional on the assumed first-order actuator model.","rationale":"The reader's weakest assumption identifies the state augmentation of Eq. (2) as the load-bearing premise. I agree: the paper's internal mathematics is coherent—Eq. (28) follows from the Bolza terminal cost with the stated sign convention, the Hamiltonian system has matching boundary conditions, and the numerical comparisons support the centre-manifold reduction. The concern is about external validity: the abstract claims a general disappearance of terminal jumps, but this is a consequence of the chosen actuator model, not a theorem about the original stiffness-as-control problem. The paper provides physical motivation (finite bandwidth, signal distortion) but does not validate the first-order model against any specific experimental actuation. This does not require changing the reader's conditional verdict, because the reader already flagged the same assumption; it strengthens the need for qualification but does not overturn the central contribution. I therefore recommend no change to the verdict.","tokens_in":30533,"tokens_out":13762,"duration_ms":120945,"concrete_test":"Re-solve the C.II minimum-work problem with a second-order actuation model, e.g. dot k = a, dot a = lambda with |lambda| <= Lambda (and analogously for u), while keeping the same work functional and harmonic penalty. Derive the boundary/transversality conditions from the resulting Hamiltonian system. If an explicit terminal condition on k_tf or on the control appears, or if boundary layers do not vanish as Lambda -> infinity, the no-jump conclusion is specific to the first-order model (2). If the structural result is unchanged, the concern is answered.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—'the problem ... that optimal protocols need terminal jumps ... completely disappears'—rests on the modeling decision in Section 2.1 to treat k_t and u_t as state variables with first-order dynamics (2), dot k = lambda, dot u = gamma, and essentially bounded controls. This converts the terminal stiffness from a control endpoint into a state endpoint: in case C.I it is prescribed by the equilibrium conditions (7); in case C.II it is free and fixed only by the transversality condition (28). No terminal jump can appear because k_t is absolutely continuous by construction. But in the original problem discussed in the literature, the stiffness (or drift) is the control; its terminal value is a protocol endpoint, not a state variable, and the overdetermination the paper describes is real. The paper's own numerical comparison S.I vs S.II (Figs. 5 and 6) shows that when stiffness is the control, the direct method produces endpoint jumps, and the state formulation agrees only after smoothing. Thus the headline is not a general resolution but a property of an augmented model whose physical validity rests on the assumption that the mechanical potential responds as a first-order integrator with bounded rate. That assumption is plausible but not validated; if the true actuation dynamics is second-order or has a finite step response, endpoint layers or jumps may reappear. The abstract should qualify the claim or demonstrate robustness to the actuator model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies optimal control of swift equilibration in underdamped Langevin dynamics, focusing on a one-dimensional nanomechanical oscillator. The central reformulation is to treat the parameters of the mechanical potential, k_t and u_t, as state variables obeying first-order control equations (2), with the actual controls being their time derivatives. This places the minimum-work and minimum-dissipation problems in the canonical Bolza form of optimal control, so that terminal costs depend only on state variables. The authors show that in this formulation the terminal stiffness is not prescribed for minimum-work transitions but is selected by the transversality condition y^(4)_tf = -x^(1)_tf/2 (Eq. 28), so no terminal jumps are needed. They derive Pontryagin first-order conditions, analyze the hard-, logarithmic-, and harmonic-penalty cases, and use centre-manifold and multiscale perturbation theory to obtain a universal slow-fast normal form (Eqs. 49-50) describing the turnpike behaviour of optimal protocols. Numerical comparisons between direct and indirect methods support the analytical results.","tokens_in":30960,"tokens_out":15848,"duration_ms":121851,"significance":"If the reformulation is accepted, the paper makes a useful conceptual clarification: the terminal-jump problem discussed in earlier optimal-control treatments of swift equilibration is an artefact of modelling the stiffness itself as the control. Treating the potential parameters as state variables with bounded rates of change removes the overdetermination and yields continuous protocols. The centre-manifold/turnpike analysis provides a compact universal description of the bulk dynamics, and the explicit distinction between minimum-dissipation (equilibrium-to-equilibrium) and minimum-work (equilibrium-to-non-equilibrium) transitions is valuable. The paper also provides reproducible numerical code and compares direct and indirect optimization methods. These are genuine strengths. The main caveat is that the headline conclusion is conditional on the state-augmentation assumption, and some parts of the asymptotic derivation are not fully transparent.","major_comments":[{"comment":"The statement that the terminal-jump problem 'completely disappears' is a property of the reformulated model with k_t and u_t as state variables obeying (2), not a theorem about the original stiffness-as-control formulation. In that original formulation (case S.II), the paper's own Figs. 5-6 show endpoint jumps. The abstract and conclusions should qualify the claim, e.g., 'in the Bolza formulation with first-order actuator dynamics' or 'for finite control bandwidth'. As written, the claim overreaches and could mislead readers about the scope of the result.","section":"Abstract; Sec. 2.1; Sec. 10.2"},{"comment":"The derivation of the universal normal form (49)-(50) rests on the solvability condition of the 8x8 system (47), specifically on the claim that ker(A^(1)^T) is spanned by {e_5,...,e_8} and that the omitted block entries 'do not play a role'. Without the explicit entries of A^(1) or a rigorous argument establishing the kernel and the projection of the non-homogeneous term, the reader cannot verify this load-bearing step. Please provide the explicit blocks in an appendix or in supplementary material, or give an independent derivation of the solvability condition.","section":"Sec. 7.4"},{"comment":"The running cost in the work functional is written as ∫(x^(3)_t + x^(6)_t/2)dt - t_f. From the Itô-lemma calculation in Section 3, the running cost should be ∫E[p_t^2]dt - t_f = ∫(x^(3)_t + (x^(6)_t)^2)dt - t_f. The term x^(6)/2 appears dimensionally inconsistent with the cumulant definitions (x^(6) is a first moment, x^(3) a second moment). Although the paper restricts to protocols with x^(6)=0, the displayed formula is a general statement and should be corrected to avoid propagating an error.","section":"Eqs. (10), (11), (13)"}],"minor_comments":[{"comment":"For minimum-work (C.II) transitions, the text explains that the terminal stiffness is free and determined by (28), but the figure captions say 'replace the boundary condition for x^(4)_tf with (28)' while referring to the equilibrium boundary conditions (56). Please state explicitly that in C.II the final state is generically not an equilibrium and that the target state is specified only by x^(1)_tf (and x^(3)_tf = 1, zero cross-correlation), not by a terminal stiffness value.","section":"Sec. 10.2; Figs. 7-8"},{"comment":"The scaling h = g^{1/4} is motivated by analogy with [45], but the paper does not explain why this particular scaling is the distinguished limit for the harmonic penalty. A short justification would help readers.","section":"Sec. 7.2"},{"comment":"The description of the numerical comparison S.I versus S.II says 'Applying a standard mean filter convolution onto the data uncovers the centre manifold'. This smoothing is a post-processing step; it would be helpful to state explicitly that the raw S.II solution is noisy and the smoothing is only for visual comparison, not part of the optimization.","section":"Sec. 10.1.1"},{"comment":"Minor typographical issues: Fig. 2 legend 'solution of first order conditions direct optimisation' is ambiguous; Sec. 8.2 labels y^(4:1)_0,t2 and x^(4:1)_0,t2 as constants but the text says 'first order corrections to stiffness and corresponding drift are constant', which is consistent; please double-check subscripts in Eqs. (45)-(46) for readability.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and presents a conceptually useful reformulation with solid numerics. The main issues are the overstatement of the no-jump conclusion in the abstract and the lack of transparency in the centre-manifold derivation. Both are fixable with revisions. The availability of code is a plus. I would be willing to review a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper is a serious candidate for a slot in your reading group. The core idea is to treat the stiffness and centre of the mechanical potential as state variables with first-order actuator dynamics (Eq. 2), so the terminal values of the potential become part of the state, not control endpoints. That single move turns the work functional into proper Bolza form, cleanly separates minimum-dissipation transitions between equilibria (C.I) from minimum-work transitions to non-equilibrium states (C.II), and removes the infamous terminal-jump problem: the final stiffness in C.II is fixed by a transversality condition (Eq. 28) rather than prescribed. The paper also gives a centre-manifold/turnpike analysis showing that optimal protocols in the bulk converge to a universal slow manifold, with exponential boundary layers at the ends. The numerics back this up: direct optimization, indirect shooting/collocation, and the centre-manifold equations agree well.\n\nWhat is actually new is the formulation, not the ingredients. Pontryagin, Fenichel, multiscale perturbation—these are standard tools. But the state-augmentation move is conceptually clean and resolves a real modeling inconsistency in the literature. The authors are also honest: they flag numerical instabilities in the minimum-work case, discuss the physical realizability of negative stiffness, and openly state that their actuator model is a first-order integrator with bounded rate.\n\nNow the soft spots, in proportion.\n\nThe biggest one is the headline. The abstract says the terminal-jump problem 'completely disappears.' That is only true inside the augmented model. If you keep the stiffness as the direct control, the overdetermination is real, and the paper's own S.I vs S.II comparison (Figs. 4-6) shows the endpoint jumps reappear. The no-jump result is therefore a property of the modeling choice, not a theorem about the original problem. The authors know this—they argue in the introduction that instantaneous potential changes are unphysical—but the abstract overstates it. A qualified claim and ideally a robustness check against a second-order actuator model would be needed.\n\nSecond, the derivation of the centre manifold in Sec. 7.4 omits the explicit blocks A(1:1,1) and A(1:1,2) and the non-homogeneous terms, with a hand-wave that they 'do not play a role.' The final result is consistent with the stationarity conditions, which is reassuring, but a referee should have the actual calculation, at least in an appendix.\n\nMinor: the GitHub repository is cited without a commit hash; reproducibility would be improved. Also, the 'universal' normal form is derived for the harmonic penalty with a specific scaling; the universality across the logarithmic and hard penalties is numerical, not analytic.\n\nBottom line: this is a solid, honest paper that deserves a serious referee. My own verdict would be conditional: fix the overstatement in the abstract, add the omitted block forms, and the scientific core is defensible. I'd take it to reading group and would cite it.","headline":"Worth refereeing: the no-terminal-jump result is a real advance but it is built into the state augmentation in Eq. (2), not a theorem about the earlier stiffness-as-control formulation.","tokens_in":31345,"tokens_out":4468,"would_cite":true,"duration_ms":42282,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["49K15","82C31","93E20"],"pacs":[],"model":"deepseek-v4-flash","headline":"A state-variable treatment of the trap stiffness eliminates the need for terminal jumps in optimal nano-oscillator protocols.","keywords":["optimal control","stochastic thermodynamics","swift equilibration","underdamped oscillator","turnpike property","centre manifold","work minimization","terminal jumps"],"falsifier":"Measure the trap stiffness during a work-optimal expansion in an underdamped optical trap with finite-bandwidth actuation: if the optimal protocol still exhibits a true discontinuity at the final time as the bandwidth is increased, the central claim fails. Alternatively, construct boundary data for which the transversality condition y^(4)_{t_f} = -x^(1)_{t_f}/2 has no admissible solution, which would limit the claimed disappearance of jumps.","tokens_in":30403,"feed_emoji":"⚙️","tokens_out":5269,"duration_ms":46950,"temperature":0.7,"pith_summary":"The paper argues that the widely discussed requirement of discontinuous 'terminal jumps' in optimal protocols for swift equilibration is an artifact of treating the mechanical-force parameters as direct controls. Once the trap stiffness and center are promoted to state variables and the actual controls are their time derivatives, the work functional takes the canonical Bolza form and minimum-work transitions impose no explicit boundary condition on the terminal stiffness; that value is selected by a transversality condition. The same reformulation cleanly separates transitions of minimal dissipation between genuine equilibria from transitions of minimal work to non-equilibrium targets. Using centre manifold theory, the paper also shows that optimal protocols generically exhibit a turnpike: in the bulk of the control horizon they stay near a universal manifold set by the running cost, with exponential boundary layers at the endpoints. Numerical solutions of the full and reduced systems support the analysis.","feed_headline":"No terminal jumps needed in optimal nano-oscillator control","feed_subtitle":"When the trap stiffness becomes a state variable, endpoint jumps become an artifact and bulk protocols become universal.","key_machinery":"The central construction is the state augmentation of the mechanical potential: stiffness k_t and offset u_t become state variables with controls λ_t=˙k_t and γ_t=˙u_t. This converts the work functional into Bolza form, with a terminal cost that is a pure state function, making the terminal stiffness an optimization variable rather than an assigned boundary condition. The analytic engine is centre manifold reduction in Fenichel coordinates, which yields the universal slow-fast system (49)–(50) whose stable and unstable modes describe boundary layers and whose stationary point is the turnpike.","core_discovery":"For a one-dimensional underdamped oscillator steered by a quadratic potential whose stiffness k_t and offset parameter u_t obey ˙k_t=λ_t, ˙u_t=γ_t, the paper shows by Pontryagin's maximum principle that the first-order conditions for minimum work reduce to a Hamiltonian system with no terminal condition on the stiffness. Stationarity of the terminal cost yields only the transversality relation y^(4)_{t_f} = -x^(1)_{t_f}/2, which implicitly selects the final stiffness. Hence the boundary-value problem is well posed without any jump in the mechanical force, and the overdetermination that produced the jump artifact disappears. A second result is that the optimality conditions admit a universal","pith_inferences":["The same state-augmentation argument should apply to other stochastic thermodynamic control problems with actuation delay or memory, such as underdamped bit erasure, where endpoint discontinuities have also been debated.","A direct experimental test could compare the measured trap-stiffness protocol against the centre-manifold prediction: for small control penalties, the bulk protocol should be independent of the penalty shape, with deviations confined to short initial and final layers.","The transversality condition implies that for work-optimal driving the final stiffness may lie far from the equilibrium stiffness, so the system should undergo a measurable uncontrolled relaxation after t_f; the paper computes the associated heat release, giving a quantitative experimental signature.","Treating the potential parameters as states effectively raises the order of the control system; similar jump-elimination may appear more generally when controls act through derivatives rather than directly on physical state variables."],"forward_implications":["Optimal swift-equilibration protocols do not need to end with a discontinuous change of the mechanical potential; smooth protocols exist with the same thermodynamic cost.","In minimum-work transitions to non-equilibrium targets, the final trap stiffness is a prediction of the theory, set by the transversality condition, not an input the experimenter must impose.","Minimum-dissipation transitions between equilibria and minimum-work transitions to non-equilibrium states are thermodynamically distinct, and the distinction is fixed by the choice of terminal cost in the Bolza formulation.","Optimal controls computed with different penalty mechanisms converge to the same universal turnpike manifold in the bulk, so stiffness-as-control results should be read as turnpike predictions, valid away from the endpoints.","The overdamped optimal control conditions are recovered in the ε→0 limit after optimization, confirming that the two limits commute for this detailed-balance dynamics."],"fun_headline_variants":["Optimal nano-oscillator control without terminal jumps","Swift equilibration: no jumps, universal protocols","Turnpike property simplifies nano-oscillator control","Reformulating optimal control kills jump artifact","Universal center manifold drives optimal nano-control"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The paper's conclusion that terminal jumps disappear rests on the modelling choice that the mechanical potential parameters are state variables whose time derivatives are the controls; if one instead treats the stiffness itself as the direct control, the overdetermination and the apparent need for jumps do not disappear.","fun_headline_variants_meta":{"raw":{"variants":["Optimal nano-oscillator control without terminal jumps","Swift equilibration: no jumps, universal protocols","Turnpike property simplifies nano-oscillator control","Reformulating optimal control kills jump artifact","Universal center manifold drives optimal nano-control"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000166,"raw_usage":{"total_tokens":1114,"prompt_tokens":790,"completion_tokens":324,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":534,"completion_tokens_details":{"reasoning_tokens":254}},"tokens_in":534,"tokens_out":324,"duration_ms":3737,"temperature":1.0,"reasoning_tokens":254,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T12:47:53.216816+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the trap stiffness during a work-optimal expansion in an underdamped optical trap with finite-bandwidth actuation: if the optimal protocol still exhibits a true discontinuity at the final time as the bandwidth is increased, the central claim fails. Alternatively, construct boundary data for which the transversality condition y^(4)_{t_f} = -x^(1)_{t_f}/2 has no admissible solution, which would limit the claimed disappearance of jumps.","supporting_citations":[],"review_version":1}