{"id":"3d381d55-bb17-4c85-b431-19a0f17b2059","arxiv_id":"2411.10215","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Thermal nonequilibrium limit cycles, previously seen in closed coronal loops, also occur in a transonic solar wind when foot-point heating exceeds about 3.5 times the global heating.","lead":"Using a 1D simulation of the solar wind from the chromosphere out to 30 solar radii, this paper shows that when foot-point heating is strong enough relative to global heating, no steady state exists: the plasma runs through repeated cycles of condensation, fallback, and evaporation every 3 to 24 hours. The result sets a quantitative bound on steady-state foot-point heating and may explain the variability of dense solar wind streams from coronal hole boundaries and plumes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The existence of TNE in a transonic wind is credible, but the quantitative threshold (alpha_c=3.5) and the 24-to-3 h period curve are set by the prescribed foot-point heating scale length lb, so the paper's 'limit on steady-state heating' is not yet a robust solar constraint.","rationale":"The reader's weakest_assumption correctly identifies the prescribed two-component heating profile, specifically the fixed small scale height lb, as the most sensitive input. My stress-test pass agrees: this is the single most load-bearing concern because the paper's novelty includes a quantitative bound on steady-state foot-point heating and a distinct cycle-period signature, both of which derive from the chosen lb. The existence proof itself does not rest on this choice—if lb is varied while keeping the strong stratification, the same qualitative behavior should emerge. I found no internal inconsistency that invalidates the central claim; the acknowledged inconsistency of the Field (1965) analysis in Section 5.4 is a diagnostic caveat, not a load-bearing step. The paper also gives due credit to prior observational and modeling context and explicitly concedes the parameter dependence in Section 8. The absence of code/data and of shown convergence tests weakens reproducibility but does not change the scientific assessment. Therefore the reader's CONDITIONAL verdict remains appropriate, and the concrete parameter sweep proposed above would settle whether the quantitative claims are robust enough to stand as solar constraints.","tokens_in":25080,"tokens_out":9995,"duration_ms":102272,"concrete_test":"Re-run the marginal and cycle cases with lb = 0.005 R_sun and lb = 0.02 R_sun (keeping lg = R_sun and Hg fixed), and measure (i) the critical alpha_c at which steady states cease, (ii) the location sn of the thermal sink, and (iii) the TNE cycle period at alpha = 4 and alpha = 6. If alpha_c shifts by more than ~20% from 3.5 or if the period-vs-alpha curve in Fig. 9 changes qualitatively when plotted against (alpha - alpha_c), then the quantitative claims in the abstract and Section 6 are specific to the chosen prescription; if the threshold and periods obey the expected scaling and remain within ~20%, the 'limit on foot-point heating' becomes a more robust, testable prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central existence claim is supported by the simulations: with the heating profile of Eq. (6), the model exhibits steady states up to alpha=3.45 and TNE limit cycles for alpha>=3.5 in a wind that becomes supersonic beyond the sonic point. The load-bearing weakness is the quantitative generalization. The threshold alpha_c is not an independent plasma property; it is set by the fixed small-scale heating length lb=0.01 R_sun through the location of the thermal sink. Near the sink, the local heating balance is H(sn) ~ -R, and because H(sn) = Hb exp(-(sn-s0)/lb) + Hg, the critical Hb scales roughly as exp((sn-s0)/lb). The cycle periods in Fig. 9 also depend on lb through the sink depth and the relaxation timescale tau_rho = hn/<u>_TR. The paper concedes in Section 8 that lb is a model parameter and that a physics-based heating model should determine the scale naturally. Therefore, while the demonstration of TNE in a transonic wind is robust to moderate changes in lb (as long as strong stratification remains), the specific quantitative bound (alpha<=3.45) and the period signature (3-24 h) are conditional on an unconstrained parameter and cannot yet be treated as universal solar limits.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents 1D field-aligned hydrodynamic simulations of a transonic solar wind along an open flux tube, using the HYDRAD code extended to 30 R_sun. The heating is prescribed as the sum of a small-scale exponential foot-point component (scale height ℓ_b=0.01R_sun) and a global exponential component (scale height R_sun), with the ratio α=H_bℓ_b/(H_gℓ_g) scanned from 0 to 10. For fixed H_g=2×10^-6 erg cm^-3 s^-1, the authors find steady-state wind solutions up to α=3.45 and, for α≥3.5, self-sustained TNE cycles in which coronal plasma condenses near a thermal sink at s≈25-50 Mm, precipitates to the chromosphere, ablates back into the corona, and relaxes before the next cycle. The simulations are supplemented by a Lagrangian 'cooling factor' Γ criterion (Eq. 14), a Field (1965) thermal-instability growth rate (Eq. 23), and an analysis of heliospheric signatures in white-light density fluctuations and O/Fe charge states. The paper concludes that open-field solar wind regions can host TNE when foot-point heating is strongly stratified, with implications for coronal-hole plumes and slow wind.","tokens_in":25328,"tokens_out":15499,"duration_ms":148608,"significance":"The central existence result is new and, if correct, extends TNE studies from closed loops to open, transonic geometries. The cycles are directly visible in the time-distance diagnostics of Figures 3-8, and the emergent nature of the loss of equilibrium (α is scanned, not imposed) is a strength. The analytic Γ and Field-growth-rate diagnostics are derived from the model equations and checked against the simulations, making the paper internally consistent. The claimed threshold α_c=3.45 and the period curve in Fig. 9 are potentially important quantitative predictions. However, their solar applicability is currently limited by the prescribed heating stratification scale ℓ_b, which the authors themselves identify in Section 8 as a model parameter rather than an emergent quantity. The paper also reports, but does not display, robustness checks against a more sophisticated emissivity and a two-fluid treatment. With these caveats, the paper makes a credible case for the existence of TNE in the solar wind, while the specific numerical limits should be treated as conditional on the adopted heating model.","major_comments":[{"comment":"The quantitative claims that steady-state solutions exist only for α≤3.45 and that TNE periods vary from ~24 h at threshold to ~3 h at α=10 are established only for ℓ_b=0.01R_sun. Because the thermal sink is located at s_n-s_0≈25 Mm, the local heating balance there contains H_b exp(-(s_n-s_0)/ℓ_b); changing ℓ_b by a factor of two changes the exponent by a factor of two and therefore shifts the critical H_b by a large factor. Section 8 correctly concedes that ℓ_b is a parameter and should emerge from a physics-based heating model. As a result, the abstract's phrase 'limits on the amount of foot-point heating that can be withstood under steady-state conditions' overstates the robustness of the result. I recommend adding a parameter study varying ℓ_b (e.g., 0.005 and 0.02 R_sun) and reporting the resulting α_c and period curves; if this is beyond the scope of the paper, the quantitative claims should be explicitly framed as conditional on the adopted stratification.","section":"Section 6 / Figure 9; Section 8"},{"comment":"The robustness checks reported in Section 8 are asserted but not shown. The statements that more sophisticated emissivity profiles produce 'little qualitative difference' and that a two-fluid treatment leaves the necessary conditions and timescales unchanged are load-bearing for the generality of the threshold and period results, yet no figures, tables, or numerical values are provided. Moreover, the same paragraph notes that the two-fluid treatment introduces coherent oscillations that 'precede and eventually trigger the thermal runaway,' which suggests a qualitative change that needs to be reconciled with the claim of unchanged conditions. Please present these tests quantitatively (e.g., a table of α_c and periods for each case) or clearly label them as preliminary tests that do not yet support the reported numbers.","section":"Section 8"}],"minor_comments":[{"comment":"At the thermal sink location s_n≈s_0+25 Mm, the foot-point heating term H_b exp(-(s_n-s_0)/ℓ_b) is approximately 2×10^-5 erg cm^-3 s^-1 for α=3.5, about an order of magnitude larger than H_g, so the statement H(s∼s_n)≈H_g is inaccurate; using the correct value changes the derived growth rate and condensation time by roughly 10-15%.","section":"Section 5.4"},{"comment":"The text gives inconsistent values for the density jump across the transition region: Section 3 says the density decreases by roughly an order of magnitude, while Section 5.2 uses n_c/n_b≈1/100.","section":"Section 3 and Section 5.2"},{"comment":"The sentence 'the cycle period of the marginal case (α=3.45) is about 24 h' should read α=3.5, since α=3.45 is the last steady-state solution and has no cycle.","section":"Section 6"},{"comment":"The caption refers to a 72 h period, but the time axis in the figure spans -20 to 40 h; please reconcile.","section":"Figure 3 caption"},{"comment":"The manuscript contains several typographical errors and an inconsistent date ('Accepted August 23, 2021'); examples include 'Obseratory', 'separatix', 'mas', 'snaphsot', and 'wherin'. A careful proofread is needed.","section":"General"},{"comment":"The notation |[Pu]|_TR used in Eq. (17) is not defined; please define the jump operator.","section":"Eq. (17)"}],"recommendation":"major_revision","confidential_remarks":"This is a solid numerical study and likely of interest to ApJ readers. The main issue is the gap between the strong quantitative claims (α_c=3.45, period curve) and the acknowledged dependence on the unconstrained heating scale ℓ_b. I would be willing to accept after a parameter study or after the claims are appropriately conditionalized."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this paper actually does show TNE limit cycles surviving in a transonic solar wind with the sonic point in the domain, on explicitly open field lines. That removes the ambiguity in Schlenker et al. (2021) and is the first result of its kind. The cycle is directly visible in the time-distance plots, the four-phase taxonomy (precipitation, ablation, relaxation, condensation) is physically plausible and matches the simulated evolution, and the period scaling from ~24 h near threshold to ~3 h at high alpha is a concrete, falsifiable signature. The ionization-state calculation (O versus Fe behaving differently) is a nice bonus that could matter for heliospheric diagnostics. Credit where due: the central existence claim is well-supported within the model, and the paper is honest about its own limitations.\n\nThe soft spots are real but mostly known to the authors. The threshold alpha ~ 3.5 and the period curve are set by the prescribed foot-point heating scale length lb = 0.01 Rsun through the location of the thermal sink. The stress-test note is right that this is not yet a robust solar constraint; the paper says so itself in Section 8. That weakness is structural, not a fixable error, and it mainly limits the quantitative generalization rather than the qualitative existence claim.\n\nTwo smaller issues. The robustness checks (more elaborate emissivity, two-fluid treatment) are asserted but not shown; I'd want those in an appendix or at least a quantitative summary. And the Section 5.4 Field (1965) instability analysis is acknowledged by the authors to rest on an equilibrium assumption that TNE by definition lacks; they use it as a diagnostic, which is fine, but they should present it that way up front rather than as a derivation of the condensation timescale. Also, no code or data is provided, which for a purely numerical paper makes independent verification impossible. That is a concrete reproducibility problem, not just a style complaint.\n\nVerdict: the paper deserves a serious referee. It is a solid, clearly-written modeling study with a genuinely new result and useful observable predictions. I would send it to review, with the expectation that the authors be asked to (a) release the HYDRAD setup or at least the input files and complete run parameters, (b) show the robustness checks rather than asserting them, and (c) reframe the Field analysis strictly as a diagnostic. The quantitative solar limit should be presented as conditional on the heating parameterization, which they already concede. I would cite this if I worked on coronal rain or solar wind sources.","headline":"First clear demonstration of TNE limit cycles in a self-consistent transonic solar wind, but the quantitative threshold and period curve are tied to the assumed foot-point heating scale height.","tokens_in":25947,"tokens_out":950,"would_cite":true,"duration_ms":11562,"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":"Strong foot-point heating can push a transonic solar wind out of steady state, triggering cycles of condensation and evaporation.","keywords":["thermal nonequilibrium","solar wind","coronal rain","coronal heating","foot-point heating","transition region","thermal instability","open magnetic field lines"],"falsifier":"Repeat the same HYDRAD simulation at $\\alpha = 3.5$ with the same two-scale heating but with the foot-point scale height doubled to $\\ell_b = 0.02\\,R_\\odot$; if a steady transonic solution is recovered, the quantitative threshold is a consequence of the chosen 7 Mm stratification rather than a general property of open-field heating. Alternatively, if a steady-state simulation at $\\alpha=3.5$ with the original profile exists using a different but still plausible emissivity that has no $T^{-1}$ branch below 0.47 MK, the instability trigger is an artifact of the radiative loss function.","tokens_in":24778,"feed_emoji":"☀️","tokens_out":7752,"duration_ms":71433,"temperature":0.7,"pith_summary":"This paper asks whether thermal nonequilibrium (TNE), the runaway cooling that makes hot coronal plasma condense and rain back to the Sun, can occur not just in closed magnetic loops but in the open magnetic field lines that feed the solar wind. It reports that it can: in a one-dimensional model of an open flux tube with a transonic outflow, steady solutions exist only up to a foot-point heating ratio of about 3.45, and slightly stronger heating triggers repeated cycles of condensation, precipitation, ablation, and re-evaporation. The cycle period is about 24 hours near the threshold and falls to roughly 3 hours for heating ratios around 10. If correct, steady-state solar wind models would have to respect a quantitative heating limit, and periodic density and ionization fluctuations in the heliosphere could be traced to this lower-corona instability.","feed_headline":"Solar wind can fall into thermal nonequilibrium cycles","feed_subtitle":"Simulations show steady outflow breaks down beyond a heating ratio of 3.45, driving repeated condensation and evaporation.","key_machinery":"The load-bearing setup is a prescribed two-scale heating profile, $H(s)=A(s)^{-1}[H_b e^{-(s-s_0)/\\ell_b}+H_g e^{-(s-s_0)/\\ell_g}]$, with $\\ell_b = 0.01\\,R_\\odot$ and $\\ell_g = R_\\odot$, and the control parameter $\\alpha = H_b\\ell_b/(H_g\\ell_g)$ that measures the ratio of foot-point to global energy deposition. The argument is carried by the 'thermal sink': a region near $s_0+h_n$ ($h_n$ about 25 Mm, the base-corona density scale height) where enhanced foot-point heating produces a local temperature minimum, an upward conductive heat flux, and radiative cooling that grows with density. The tipping point is the emissivity break at $T' = 10^{5.67}$ K, below which the radiative loss function $\\Lambda(T)$ becomes $\\propto T^{-1}$, so that cooling accelerates as the plasma cools; the paper uses the Field condensation-mode growth rate to estimate the runaway timescale and defines a 'cooling factor' $\\Gamma = \\langle \\tau_d^{-1}\\rangle/\\langle \\tau_u^{-1}\\rangle$ that must stay below 1 for a steady state to exist.","core_discovery":"On the paper's own terms, the discovery is that thermal nonequilibrium is not restricted to closed coronal loops: it can develop in a steadily outflowing, transonic solar wind when the heating is split into a small-scale component concentrated near the footpoints (scale height $0.01\\,R_\\odot$, about 7 Mm) and a global component (scale height $R_\\odot$) that accelerates the wind. Holding the global heating at $H_g = 2\\times10^{-6}\\,\\mathrm{erg\\,cm^{-3}\\,s^{-1}}$, the simulations find steady states for heating ratio $\\alpha = H_b \\ell_b/(H_g \\ell_g)$ up to 3.45, corresponding to $H_b = 6.9\\times10^{-4}\\,\\mathrm{erg\\,cm^{-3}\\,s^{-1}}$. A 1.5% increase to $\\alpha = 3.5$ ($H_b = 7.0\\times10^{-4}$) removes the steady state: the extra heating raises the density in the lower corona, which slows the wind and weakens the enthalpy flux, so that a radiation-dominated 'thermal sink' near the density scale height (about 25 Mm) can no longer be balanced by conduction. The plasma there cools below the emissivity break temperature $T' = 10^{5.67}\\,\\mathrm{K}$, where radiative loss increases as temperature falls, and a thermal instability produces a falling condensate; the subsequent evaporation, relaxation, and re-condensation repeat as a limit cycle. The period falls from about 24 h near threshold to about 3 h at $\\alpha = 10$, and the paper argues that the condensates' formation height directly reflects the scale height of the foot-point heating.","pith_inferences":["Editorial extension: the threshold $\\alpha\\approx3.45$ is likely to move if the small-scale heating layer is thickened; repeating the simulation with $\\ell_b$ varied by a factor of two and checking whether steady transonic solutions reappear at $\\alpha=3.5$ would map the sensitivity of this limit to heating stratification.","Editorial extension: because the falling condensate is denser and slower than the ambient wind, it should act as a transient reflector of Alfvén waves, briefly turning the sub-condensate volume into a resonant cavity; this could imprint a period of order the cycle time on wave-driven heating and on first-ionization-potential fractionation, an idea the authors suggest but do not model.","Editorial extension: if open-field TNE actually powers slow-solar-wind variability, then density and charge-state periodicity at 3-24 h periods should appear in in-situ solar wind time series, a prediction that can be tested against existing spacecraft data.","Editorial extension: the cycle-period versus heating-ratio relation offers a remote diagnostic, because measuring a quasi-periodic coronal-rain or white-light fluctuation period would then constrain the foot-point heating rate without needing to resolve the heating scale directly."],"forward_implications":["Steady-state solar wind models must respect an upper bound on foot-point heating: for the global heating used here, heating ratios above about 3.45 cannot be made steady on open field lines.","TNE cycle periods carry a signature of the heating rate: about 24 h at the threshold, dropping to roughly 5 h at $\\alpha=4$ and asymptoting near 3 h for $\\alpha\\ge10$.","The height at which condensates form along open field lines is set by the scale height of the foot-point heating, so coronal-rain observations could be inverted to infer how heating is stratified in the lower corona.","Each cycle launches outward-propagating density and acoustic disturbances that reach white-light coronagraph fields of view with periods of about 3-4 hours, a potential source of periodic structures in streamers and pseudostreamers.","Time-dependent nonequilibrium ionization calculations show oxygen charge states return near baseline while iron remains suppressed with cyclic fluctuations, so combined oxygen/iron charge states can diagnose both the heating and the non-steadiness of the lower corona."],"supporting_citations":[{"why":"Supplies the HYDRAD code used to integrate the field-aligned Navier-Stokes equations with adaptive refinement.","marker":"Bradshaw & Cargill 2013"},{"why":"Supplies the piecewise power-law coronal emissivity whose $T^{-1}$ branch below $10^{5.67}$ K sets the thermal-runaway threshold.","marker":"Klimchuk et al. 2008"},{"why":"Provides the closed-loop TNE scaling that the paper contrasts, showing why the open-field outflow does not inhibit TNE.","marker":"Klimchuk & Luna 2019"},{"why":"Established the baseline solar wind model, equations, and geometry that this study modifies to a single-fluid treatment.","marker":"Scott et al. 2022"},{"why":"Gives the thermal-instability condensation-mode growth rate used to estimate the condensation phase timescale.","marker":"Field 1965"},{"why":"Found that increased foot-point heating slows, cools, and densifies the wind, the behavior the paper reproduces and extends.","marker":"Grappin et al. 2011"},{"why":"Earlier 2.5D simulation with a self-consistent transonic wind in which condensates formed near a helmet streamer flank, motivating the open-field question.","marker":"Schlenker et al. 2021"},{"why":"Wave-turbulence-driven heating model whose naturally stratified heating is cited as the likely physical realization and which itself produced TNE cycles in loops.","marker":"Downs et al. 2016"},{"why":"Provides the stream-limited heat flux formula used to saturate conduction in the low-density corona.","marker":"Cowie & McKee 1977"}],"fun_headline_variants":["Solar wind can enter thermal nonequilibrium cycles","Wind heating ratio 3.45 triggers thermal cycles","Solar wind loses steady state via thermal cycles","Condensation-evaporation cycles found in transonic wind","Thermal nonequilibrium cycles emerge in solar wind"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the assumption that coronal heating near the footpoints is stratified on a small scale of about 7 Mm on top of a global scale of a solar radius; the paper itself notes that in a physics-based heating model this scale should emerge rather than be prescribed, and if the Sun's heating is stratified differently the quantitative threshold and the very occurrence of TNE in open-field wind could change.","fun_headline_variants_meta":{"raw":{"variants":["Solar wind can enter thermal nonequilibrium cycles","Wind heating ratio 3.45 triggers thermal cycles","Solar wind loses steady state via thermal cycles","Condensation-evaporation cycles found in transonic wind","Thermal nonequilibrium cycles emerge in solar wind"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000263,"raw_usage":{"total_tokens":1696,"prompt_tokens":1135,"completion_tokens":561,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":751,"completion_tokens_details":{"reasoning_tokens":488}},"tokens_in":751,"tokens_out":561,"duration_ms":5965,"temperature":1.0,"reasoning_tokens":488,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:51:33.922337+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the same HYDRAD simulation at $\\alpha = 3.5$ with the same two-scale heating but with the foot-point scale height doubled to $\\ell_b = 0.02\\,R_\\odot$; if a steady transonic solution is recovered, the quantitative threshold is a consequence of the chosen 7 Mm stratification rather than a general property of open-field heating. Alternatively, if a steady-state simulation at $\\alpha=3.5$ with the original profile exists using a different but still plausible emissivity that has no $T^{-1}$ branch below 0.47 MK, the instability trigger is an artifact of the radiative loss function.","supporting_citations":[{"cited_title":"J., & Cargill, P","cited_arxiv_id":null,"evidence_quote":"Supplies the HYDRAD code used to integrate the field-aligned Navier-Stokes equations with adaptive refinement."},{"cited_title":"B., Bradshaw, S","cited_arxiv_id":null,"evidence_quote":"Established the baseline solar wind model, equations, and geometry that this study modifies to a single-fluid treatment."},{"cited_title":"M., & Pantellini, F","cited_arxiv_id":null,"evidence_quote":"Found that increased foot-point heating slows, cools, and densifies the wind, the behavior the paper reproduces and extends."},{"cited_title":"A., & Velli, M","cited_arxiv_id":null,"evidence_quote":"Wave-turbulence-driven heating model whose naturally stratified heating is cited as the likely physical realization and which itself produced TNE cycles in loops."},{"cited_title":"L., & McKee, C","cited_arxiv_id":null,"evidence_quote":"Provides the stream-limited heat flux formula used to saturate conduction in the low-density corona."}],"review_version":1}