{"id":"b1909d11-17c1-41e6-b1cb-0953f33402a6","arxiv_id":"1908.09028","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A WKB analysis gives the change in establishment probability and mean establishment time for a stochastic population under temporary and periodic environmental perturbations.","lead":"This paper computes, using a WKB semi-classical approximation, how a temporary or periodic environmental change alters the probability and mean time for a small population to establish. The analytical formulas could help predict when species, cells, or other stochastic populations switch states in fluctuating environments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (21) depends on an asserted optimal path; the paper never rules out lower-action paths that complete switching during the perturbation window.","rationale":"The reader correctly flags the temporary-perturbation optimal-path ansatz as the weakest point. I refine it: the piecewise structure is largely forced by Hamilton's equations, so the real missing step is the global minimization over the switching/crossing time, including trajectories that establish during the perturbation. The paper asserts rather than proves this, and because the result is an exponential probability, the omission is load-bearing. The concrete numerical/analytical test would settle it. The adiabatic and weak-periodic results are better supported (the adiabatic prefactor is parameter-free and agrees with simulations), so I do not recommend rejection; a conditional verdict with a request for this variational check is appropriate. Agreement is partial because the reader's stated concern about the concatenation is somewhat too broad; the concatenation itself is derivable, while the missing global-optimality check is the precise issue.","tokens_in":13625,"tokens_out":48129,"duration_ms":469669,"concrete_test":"For representative parameters (e.g., f0=0.40, F below the deterministic-establishment threshold, and T spanning the range of Fig. 3), compute the boundary-path action S_boundary = ∫_{q1}^{q_p1} p0 dq + ∫_{q_p1}^{q2} pp(q;E') dq - E' T, where E' solves T = ∫_{q_p1(E')}^{q2} dq / [λp e^{pp} - μp e^{-pp}], and compare it with Eq. (21). Also, if desired, solve the discretized time-dependent variational problem by shooting or optimal-control methods. If S_boundary is never smaller than Eq. (21) and all numerical minimizations match Eq. (21), the asserted optimality is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The temporary-perturbation result S = S0 - E_p T - ∫_{q_p1}^{q_p2} [p0(q)-pp(q;E_p)] dq rests on the claim in Sec. IV A that the optimal path is the concatenation of the unperturbed H0=0 trajectory before/after the perturbation and a single Hp=E_p trajectory during it, with E_p fixed by Eq. (20). The concatenation itself follows from continuity of q and p and conservation of each autonomous Hamiltonian, but this does not establish global optimality. In particular, a path that crosses the establishment threshold q2 during the favorable period (so the final time is free and no post-perturbation segment is needed) is not considered. The paper states 'It turns out that...' and gives no transversality or variational argument showing such paths have larger action. Since the establishment probability is exponentially sensitive to S (Eq. 15), any lower-action alternative would make the predicted EP increase exponentially wrong. This is not merely a prefactor issue; it concerns the exponent itself. The gap is most pressing at larger F or T, where the action reduction is large and the WKB exponent approaches O(1), the regime where crossing during the perturbation could plausibly compete with the constructed path.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies noise-induced establishment in a stochastic Verhulst-Allee model with immigration and a strong Allee threshold. In the WKB limit N≫1, it computes the leading-order action for switching from the metastable pre-established state to the established state. For a temporary environmental change, the authors propose a piecewise optimal path constructed from the unperturbed zero-energy trajectory before and after the perturbation and a constant-energy trajectory H_p=E_p during the perturbation; the energy is fixed by the duration condition Eq. (20), and the action is S=S0-E_pT-∫_{q_p1}^{q_p2}(p0-pp)dq (Eq. 21). Near bifurcation this reduces to Eq. (25), which depends only on the product FT. For a periodic perturbation they derive the linear-in-amplitude correction Eq. (32), an adiabatic expression with prefactor Eq. (36), and a square-wave analogue Eq. (39). The analytical exponents are compared with Gillespie simulations for time-dependent rates, using arbitrary prefactors for plotting.","tokens_in":13896,"tokens_out":28141,"duration_ms":286388,"significance":"If the temporary-perturbation ansatz is correct, the paper is a useful contribution: it gives explicit closed-form expressions for the exponential change of establishment probability and mean establishment time under two types of environmental variability, identifies scaling laws (FT dependence near bifurcation, amplitude-dominated adiabatic limit, exponential suppression with frequency), and extends earlier extinction work to switching. The derivation is self-contained from the master equation, the analytical exponents contain no fitted parameters, and the simulations independently implement a time-dependent Gillespie algorithm. The main caveat is that the central temporary-perturbation result rests on an unproved optimal-path ansatz; until that is justified, the significance of Eq. (21) is conditional.","major_comments":[{"comment":"The optimality of the piecewise path is asserted, not proved. The text says \"It turns out that...\" and then constructs the path by matching the H0=0 trajectory before and after the perturbation with a single Hp=Ep trajectory; continuity and conservation produce a candidate, not a minimizer of the action in Eq. (14). In particular, the paper never compares against paths that reach q2 during the favorable window, so that the post-perturbation segment is absent and the final time is free, nor against paths carrying nonzero H0 energy before the perturbation. Because Eq. (15) makes the EP exponentially sensitive to S, any lower-action competitor changes the predicted exponent, not just a prefactor. This is not a pedantic issue: at larger F or T, where the action reduction is large and S can approach O(1), the excluded alternatives are precisely the ones that might compete. Please add a variational or transversality argument, or a numerical minimization of Eq. (14) over path families, or state a regime in which such alternatives are provably subdominant.","section":"Sec. IV A, Eqs. (18)-(21)"},{"comment":"The manuscript does not state the domain of (F,T) for which the duration equation has a solution and for which q_p1 and q_p2 lie in the physical interval (q1,q2). The near-bifurcation limit gives the condition FT<2, but the full Eq. (20) is not analyzed. Since the central formula (21) is meaningful only when such an intersection exists, the domain of validity of the temporary-perturbation result is currently undefined; the numerical failures at large F,T are attributed to WKB breakdown without separating the two possible causes, namely absence of an admissible solution versus action becoming O(1). Please provide this analysis or at least an explicit statement of the assumed parameter range.","section":"Sec. IV A, Eqs. (19)-(20)"}],"minor_comments":[{"comment":"As written, σ(t)=1-δ tanh(δt/2) describes a trajectory that starts at q2=1+δ and ends at q1=1-δ, opposite to the stated boundary conditions in the same paragraph. The intended optimal path from q1 to q2 requires σ(t)=1+δ tanh(δt/2). The phase in Eq. (31) suggests the authors used the plus sign in the actual integration; please correct Eq. (8) and the time parameterization.","section":"Eq. (8)"},{"comment":"The symbol T is used both for the actual perturbation duration and for the rescaled duration T/δ. Please use a distinct symbol, e.g. \\mathcal T, in Eqs. (24)-(25) and in the caption of Fig. 4, because the condition FT<2 is easily misread as a condition on the physical duration.","section":"Sec. IV A.1, Eqs. (24)-(25)"},{"comment":"In the paragraph introducing the adiabatic approximation, \"w≪δ\" should read \"ω≪δ\".","section":"Sec. IV B.2"},{"comment":"The word \"crushed lines\" should be \"dashed lines\" or another standard term.","section":"Fig. 2 caption"},{"comment":"The multiplicative temporary perturbation is dismissed with \"similar results (not shown)\", although the paper claims to treat both additive and multiplicative variability. Please include the analog of Eqs. (18)-(21) for the multiplicative case, or state explicitly that the calculation is identical in form with the corresponding λp and μp.","section":"Sec. IV A"},{"comment":"For the temporary perturbation the paper computes only the change in EP; the MTE is asserted in footnote [46] to be almost unaffected. Please make this explicit in Sec. IV A when the temporary protocol is introduced, since the abstract mentions both EP and MTE.","section":"Abstract and Sec. IV A"},{"comment":"The comparisons use linear vertical scales together with arbitrary multiplicative prefactors. Since the analytical claims are for exponents only, plotting N^{-1}\\log ΔP or \\log τ would make the agreement and its limits more transparent, especially where the prefactor is not predicted.","section":"Figs. 3-5"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely acceptable after the authors provide a variational justification or numerical verification for the piecewise optimal path in Sec. IV A, and after correcting the sign inconsistency in Eq. (8). The missing proof is the only substantive obstacle; the remaining issues are presentation and parameter-range questions. I would not reject the paper because the rest of the derivations and the simulation validation are consistent."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know before reading: the paper is a legitimate extension of the extinction WKB machinery to establishment (upward switching) under temporary and periodic environmental changes; and the headline temporary-perturbation formula Eq. (21) depends on an optimal-path ansatz that is asserted, not derived. If you work on stochastic population dynamics, the adiabatic and periodic results are worth your time; the temporary case is plausible but not pinned down.\n\nWhat's actually new: to my knowledge this is the first systematic treatment of establishment probability and mean time under time-dependent environments in this model. The periodic linear correction Eq. (32) is a clean application of existing theory, and the adiabatic result Eq. (36) with its prefactor is a solid, testable prediction. The near-bifurcation reduction Eq. (25) is a nice simplification showing a scaling collapse on FT. The simulations use a modified Gillespie with time-dependent rates and the code description is adequate. All of this is credible.\n\nThe soft spot is exactly where the stress-test note lands. In Sec IV A the authors say 'It turns out' the optimal path is the concatenation of the H0=0 trajectory, a constant-energy Hp=Ep trajectory, and the H0=0 trajectory again, with Ep set by the duration T. They never rule out a path that completes the switch during the perturbation, so the post-perturbation segment is unnecessary. That is not a prefactor worry; it changes the exponent. The paper doesn't give a variational argument or a check. The numerical agreement in Fig. 3 confirms the exponential scaling, but the theory is multiplied by an arbitrary prefactor, so it can't distinguish between Eq. (21) and another exponent in the tested range. The authors themselves note the theory fails when the action becomes O(1), which is exactly when a crossing-during-perturbation path might compete. So the central temporary result is under-supported.\n\nThe periodic and adiabatic parts don't have this problem; they are standard and the simulations support them without hand-tuned prefactors in the adiabatic case (Fig. 7). The paper cites its own prior work heavily, but that's appropriate because the method is an extension rather than a borrowed result.\n\nBottom line: this deserves a serious referee. I'd send it to review with a request to either justify the ansatz (even a simple upper bound on alternative actions) or test it numerically by simulating the actual path probabilities. The paper is honest about its limitations and the physics community in population dynamics would cite the periodic and adiabatic results even if the temporary case turns out to need revision.","headline":"New WKB results for switching under time-varying environments, but the temporary-perturbation exponent rests on an unproved path ansatz that needs testing or proof.","tokens_in":14330,"tokens_out":2400,"would_cite":true,"duration_ms":24462,"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":"A changing environment exponentially alters the odds of population establishment.","keywords":["population switching","establishment probability","mean establishment time","WKB approximation","Allee effect","time-varying environment","stochastic population dynamics","semi-classical action"],"falsifier":"Numerically minimize the action functional of Eq. (14) over all smooth heteroclinic paths for a temporary perturbation, without assuming the piecewise constant-energy form, and compare the minimizing action with Eq. (21); any discrepancy beyond the WKB accuracy would falsify the concatenation assumption and hence the predicted establishment probability.","tokens_in":13424,"feed_emoji":"🌱","tokens_out":7475,"duration_ms":72056,"temperature":0.7,"pith_summary":"This paper analyzes a stochastic population in a pre-established state, subject to a strong Allee effect and constant immigration, and asks how a deterministically time-varying environment changes the rare event of establishment. It claims that under a temporary favorable change, the action governing the exponential establishment probability is reduced by an explicit piecewise-path formula, $S = S_0 - E_p T - \\int_{q_p^1}^{q_p^2} [p_0(q) - p_p(q;E_p)]\\,dq$, and that under a weak periodic perturbation the mean establishment time is exponentially shortened by the frequency-dependent correction of Eq. (32). These results give ecologists and cell biologists a way to compute, within exponential accuracy, how pulse-like or seasonal environmental changes influence switching rates. The paper verifies the formulas against a modified Gillespie algorithm with explicitly time-dependent rates.","feed_headline":"Time-varying environments reshape population establishment odds","feed_subtitle":"Pulse-like and periodic changes rewrite the exponential rate of rare establishment events, confirmed by simulation.","key_machinery":"The machinery is the WKB Hamilton–Jacobi formulation of the master equation, with Hamiltonian $H(q,p,t) = [\\lambda(q,t)-\\mu(q,t)e^{-p}](e^p-1)$, whose zero-energy heteroclinic trajectory between the fixed points $(q_1,0)$ and $(q_2,0)$ carries the optimal switching path, and the action $S = \\int p\\,dq - H\\,dt$ sets the exponential rate. For the temporary perturbation, the new object is the constant-energy perturbed trajectory $p_p(q;E_p)$, whose energy is fixed by requiring the perturbed segment to last exactly the perturbation duration $T$, via Eq. (20). For the periodic case, the key step is expanding the action to first order in the amplitude and minimizing over the initial phase $t_0$, which yields the closed-form correction of Eq. (32).","core_discovery":"The central discovery is that the semi-classical WKB approach, previously used for extinction under time-modulated environments, extends to population switching, with the action computed along a heteroclinic path of an explicitly time-dependent Hamiltonian. For a temporary environmental change, the optimal path is the concatenation of the unperturbed zero-energy trajectory before and after the perturbation and a constant-energy trajectory of energy $E_p$ during it; the action is $S = S_0 - E_p T - \\int_{q_p^1}^{q_p^2}[p_0(q) - p_p(q; E_p)]\\,dq$, with $E_p$ fixed by the duration constraint of Eq. (20). For a weak sinusoidal perturbation, the leading correction to the action, after minimizing over the phase of the perturbation, is $\\Delta S = -\\varepsilon \\pi \\sqrt{1-\\delta^2}\\,\\mathrm{csch}(\\pi\\omega/\\delta)\\,\\sinh[(\\omega/\\delta)\\arcsin(\\delta)]$, so the establishment time decreases exponentially with the perturbation amplitude in a frequency-dependent way. The paper verifies the exponential scaling of both results numerically.","pith_inferences":["I infer that the piecewise constant-energy ansatz generalizes to a sequence of temporary switches: the action would accumulate terms $-E_i T_i$ plus pairwise area corrections, providing a route to model randomly switching environments by averaging over switch protocols.","Equation (32)'s factor $\\mathrm{csch}(\\pi\\omega/\\delta)$ suggests a cutoff frequency $\\omega \\sim \\delta$ beyond which periodic forcing becomes exponentially ineffective; a microbial evolution experiment with controlled periodic stressors could test whether the measured switching rate saturates to the unperturbed rate at high frequency.","The minimization over the initial phase $t_0$ in the periodic case is analogous to a stochastic-resonance condition: the optimal escape attempt locks to the phase of maximal perturbation, which may be observable in time-resolved switching statistics accumulated over many cycles.","The near-bifurcation result that only the product $FT$ matters hints at a broader principle for weak-noise escape: the leading exponential effect of a temporary change may be governed by the time-integrated perturbation, a conjecture worth testing in other escape models."],"forward_implications":["A temporary favorable change of magnitude $F$ and duration $T$ exponentially raises the establishment probability, with the exponent given by Eq. (21); longer or stronger perturbations monotonically reduce the action.","Close to the bifurcation limit, the action depends only on the product $FT$ (as long as $FT < 2$), so pulses with the same duration–magnitude product have the same effect on the establishment probability.","A weak periodic perturbation exponentially reduces the mean establishment time; the correction is linear in the amplitude and is largest when the perturbation is slow compared to the relaxation rate.","In the adiabatic limit, the shape of the periodic perturbation, sinusoidal versus square wave, does not change the leading-order exponent; only the amplitude matters.","The same theoretical framework applies to phenotypic switching in gene regulatory networks under extrinsic fluctuations, where switching times depend sensitively on the magnitude, frequency, and duration of environmental variations."],"supporting_citations":[{"why":"Provide the unperturbed mean time to establishment and establishment probability, including the pre-exponential factor, used as the baseline for the perturbed results.","marker":"[7, 9]"},{"why":"Supplies the time-dependent WKB framework and adiabatic averaging for extinction under a time-modulated environment, which the paper extends to switching.","marker":"[31]"},{"why":"Gives the temporary-perturbation WKB action formula that the paper generalizes to arbitrary perturbation magnitude and duration.","marker":"[32]"},{"why":"Provides the WKB Hamilton–Jacobi formalism of metastable decay, including the heteroclinic path and the action integral.","marker":"[37]"},{"why":"Establish the method of minimizing the action over the initial phase $t_0$ for periodic perturbations, used to obtain Eq. (32).","marker":"[22, 23, 47]"},{"why":"Gives the exact sampling criterion for the Gillespie algorithm with time-dependent rates, on which the modified simulation is based.","marker":"[40]"},{"why":"Provides the standard Gillespie algorithm that the modified time-dependent version builds on.","marker":"[39]"}],"fun_headline_variants":["Environmental pulses and cycles shift establishment probability","Rare establishment times respond exponentially to environmental swings","Periodic environments exponentially speed establishment","Time-varying environments bend population switching odds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The temporary-perturbation formula assumes that the optimal switching path is exactly the unperturbed zero-energy trajectory before and after the change joined to a constant-energy trajectory during the change, with the energy chosen only so that the path lasts the perturbation duration; the paper does not prove that no other path gives a smaller action.","fun_headline_variants_meta":{"raw":{"variants":["Environmental pulses and cycles shift establishment probability","Rare establishment times respond exponentially to environmental swings","Periodic environments exponentially speed establishment","Time-varying environments bend population switching odds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001064,"raw_usage":{"total_tokens":4446,"prompt_tokens":917,"completion_tokens":3529,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":533,"completion_tokens_details":{"reasoning_tokens":3477}},"tokens_in":533,"tokens_out":3529,"duration_ms":24809,"temperature":1.0,"reasoning_tokens":3477,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:23:48.670050+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically minimize the action functional of Eq. (14) over all smooth heteroclinic paths for a temporary perturbation, without assuming the piecewise constant-energy form, and compare the minimizing action with Eq. (21); any discrepancy beyond the WKB accuracy would falsify the concatenation assumption and hence the predicted establishment probability.","supporting_citations":[{"cited_title":"Lande, The American Naturalist 142, 911 (1993)","cited_arxiv_id":null,"evidence_quote":"Supplies the time-dependent WKB framework and adiabatic averaging for extinction under a time-modulated environment, which the paper extends to switching."},{"cited_title":"Kamenev, B","cited_arxiv_id":null,"evidence_quote":"Gives the temporary-perturbation WKB action formula that the paper generalizes to arbitrary perturbation magnitude and duration."},{"cited_title":"Billings and E","cited_arxiv_id":null,"evidence_quote":"Provides the WKB Hamilton–Jacobi formalism of metastable decay, including the heteroclinic path and the action integral."},{"cited_title":"Be’er and M","cited_arxiv_id":null,"evidence_quote":"Provides the standard Gillespie algorithm that the modified time-dependent version builds on."}],"review_version":1}