{"id":"db3a1c34-b86c-4d94-b148-9a2c646a5122","arxiv_id":"2509.08605","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Fitting a modified ATLCW quiet-Sun model to ALMA, Metsähovi, and Nobeyama observations yields prominence densities 60-163 times higher and temperatures 155-163 times lower than the quiet Sun, with non-hydrostatic (magnetostatic) support favored.","lead":"By modifying a standard quiet-Sun atmospheric model to match millimeter-wavelength radio images, this paper infers that solar prominences are cooler and denser than their surroundings and appear in absorption. It also argues that the data favor magnetic rather than gas-pressure support of prominences, a result relevant to future ALMA studies of the solar chromosphere.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed 'highly statistically significant' preference for the non-hydrostatic model is contradicted by Table 2: Δχ² is only 2.3–2.7 for one extra parameter (p≈0.10–0.13), so the magnetic-support conclusion is not statistically established.","rationale":"The reader selected dataset heterogeneity as the weakest assumption. That concern is real: five measurements come from different prominences, instruments, epochs, and beam sizes, and beam dilution is acknowledged but not corrected. However, the statistical significance claim in Section 4.2 is a more direct and internal problem: even if the dataset were homogeneous, the reported χ² minima do not support the statement that the non-hydrostatic model is 'highly statistically significant' preferred. Since the hydrostatic model is nested within the non-hydrostatic model (one extra free parameter), the improvement of Δχ²≈2.3–2.7 is marginal (p≈0.10–0.13), far below conventional thresholds. This directly affects the central interpretive claim that prominence stability is most likely maintained by the magnetic field rather than by hydrostatic equilibrium. The additional issue that the 'hydrostatic' case is actually an isobaric (constant-pressure) model rather than a true hydrostatic stratification reinforces the need for caution, but the concrete fix is to report a proper significance test. Because the paper still provides a valuable new ALMA measurement and a plausible qualitative picture (cooler, denser, absorbing prominences), the conditional verdict stands; the authors should retract or soften the statistical-preference claim unless a proper test supports it.","tokens_in":18467,"tokens_out":12388,"duration_ms":442479,"concrete_test":"Perform a likelihood-ratio test using the χ² values in Table 2. With 5 data points and assumed 5% errors, compute p = 1 - CDF_χ²(Δχ²; df=1). For procedure a, Δχ² = 11.69 - 14.39 = 2.70 (p≈0.10); for procedure b, Δχ² = 5.84 - 8.11 = 2.27 (p≈0.13). If p>0.05 for both procedures, replace the 'highly statistically significant' claim in Section 4.2 with a statement that the non-hydrostatic model is not significantly preferred, and downgrade the magnetostatic-support conclusion accordingly. Optionally also report AIC/BIC, which penalize the extra free parameter more strongly than the raw Δχ².","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing weakness is the claim in Section 4.2 that the non-hydrostatic fit is 'highly statistically significant' relative to the hydrostatic fit. Table 2 gives χ²_min(hydro)=14.39 (procedure a) and 8.11 (b), while χ²_min(non-hydro)=11.69 (a) and 5.84 (b). These models are nested: the hydrostatic case fixes f_T=1/f_n, leaving one free parameter, while the non-hydrostatic case allows f_n and f_T independently, adding one parameter. With five measurements and the stated 5% errors, the likelihood-ratio statistic is Δχ²=2.70 (a) and 2.27 (b) on one degree of freedom, corresponding to p≈0.10 and p≈0.13. This is not significant at the 5% level, so the phrase 'highly statistically significant' is not supported by the paper's own numbers. In fact, for procedure (a) the absolute χ² values (11.69 with 3 dof; 14.39 with 4 dof) indicate a poor fit even for the preferred model, so the data do not strongly favor either model. The conclusion that magnetic (non-hydrostatic) support is statistically preferred therefore rests on an untested and, on the face of Table 2, unsupported significance claim. In addition, the 'hydrostatic equilibrium' baseline is implemented as constant gas pressure (f_T=1/f_n, so P_PR=P_QS), not as a self-consistent hydrostatic stratification with dP/dh=-ρg; this mislabels the baseline model, though correcting the significance analysis is the key check.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper adapts the 1D semi-empirical ATLCW quiet-Sun model to compute the brightness temperature of solar prominences on the disk at millimeter wavelengths. The authors use ALMA single-dish Band 3 and Band 6 measurements, supplemented by Metsähovi 8.1 mm and Nobeyama 3.1/8.3 mm observations, and apply multiplicative factors to the electron density and temperature of the ATLCW model within a 40,000–50,000 km layer. They fit the factors by χ² minimization under two assumptions: a 'hydrostatic' case with f_T = 1/f_n and a 'non-hydrostatic' case with f_n and f_T free. The best fits give PR densities 60–68 times (hydrostatic: 159–163 times) the quiet-Sun value and temperatures roughly 155–163 times lower, with the non-hydrostatic model yielding, the authors claim, a statistically significant better fit. From this they conclude that prominences appear in absorption at mm wavelengths, that their stability is likely magnetic rather than hydrostatic, and that thermal bremsstrahlung is the dominant mm emission mechanism.","tokens_in":18742,"tokens_out":6279,"duration_ms":52598,"significance":"The radiative-transfer calculation is standard and the data are presented transparently in Table 1, with a useful sensitivity check excluding Metsähovi discussed in Section 5. If the model comparison were statistically sound, the paper would provide a simple semi-empirical description of prominence mm emission and a testable prediction of absorption contrast. However, the central statistical claim is not supported by the paper's own numbers, the 'hydrostatic' baseline is not a true hydrostatic stratification, and the five measurements are heterogeneous in structure and beam size. As presented, the work is a promising framework rather than an established result.","major_comments":[{"comment":"The statement that the difference in χ²_min between the hydrostatic and non-hydrostatic fits is 'highly statistically significant' is not supported. The models are nested with one extra free parameter; from Table 2, Δχ² = 14.39 − 11.69 = 2.70 (procedure a) and Δχ² = 8.11 − 5.84 = 2.27 (procedure b), corresponding to p ≈ 0.10 and p ≈ 0.13 on one degree of freedom. These are not significant at the 5% level, and the absolute χ² values for procedure a indicate poor fits even for the preferred model (χ²_min = 11.69 with 3 degrees of freedom). The conclusion that the data 'strongly favor' the non-hydrostatic/magnetic-support model must be revised or supported by a proper likelihood-ratio test with reported p-values.","section":"§4.2, Table 2"},{"comment":"The 'hydrostatic equilibrium' case is not modeled as hydrostatic equilibrium. Setting f_T = 1/f_n imposes constant gas pressure (n_e T_e = const.), not the hydrostatic relation dP/dh = −ρg with a stratified atmosphere. A true hydrostatic model would have scale-height pressure and density profiles varying with height, which is incompatible with applying a single uniform f_T across the 40,000–50,000 km layer. Because the comparison between the two cases is the basis for the paper's main conclusion, the baseline model must be either correctly implemented or explicitly relabeled as 'constant pressure' and the conclusions adjusted accordingly.","section":"§4.1, Table 2"},{"comment":"The joint fit treats five measurements of different prominences from different instruments and epochs as a single homogeneous dataset, despite the paper's own statement in Section 4.1 that 'ALMA, Metsähovi, and Nobeyama all observed a different PR structure.' Beam sizes range from 17″ to 144″ and beam-convolution effects are acknowledged but not corrected. Under these conditions, the fitted f_n and f_T values and the model comparison could be biased by unknown structure-to-structure variations and by resolution-dependent contrast loss. A quantitative treatment—for example, convolving the modeled profiles with the respective beams or fitting per-structure normalization parameters—is needed before the reported factors can be taken at face value.","section":"§3, §4.1"},{"comment":"The statement that the good agreement 'confirms that thermal bremsstrahlung is indeed the dominant radiation mechanism' is partly circular. Thermal bremsstrahlung is assumed in the optical-depth expression (Eq. 1), and the density/temperature factors are fitted to the same measured brightness temperatures used for the comparison. The agreement therefore validates internal consistency but does not independently confirm the emission mechanism; a direct test would require comparing against an alternative emission model or using independent constraints on n_e and T_e.","section":"§5"}],"minor_comments":[{"comment":"The ranges '60–68' and '159–163' for density factors should be labeled by stability assumption at first mention to avoid apparent inconsistency with the abstract's combined '60–163' range.","section":"Abstract and §4.2"},{"comment":"The Nobeyama 3.10 mm entry displays '= 𝑻𝐛(PR)' in an unusual way; the column should be cleaned up, and the date range '1984-07-16 (start) ... 1984-07-22 (end)' should be presented more clearly.","section":"Table 1"},{"comment":"The symbols for the three instruments are listed without a legend definition in the caption; add the symbol definitions as is done for Figure 3.","section":"Figure 2 caption"},{"comment":"The statement that only two height points fall in the 40,000–50,000 km range, followed by cubic-spline interpolation by a factor of 100, should be accompanied by a quantitative convergence test showing that the interpolated profiles are stable over the wavelength range of interest; the text asserts this verbally but gives no numerical check.","section":"§2"},{"comment":"Equation (5) could benefit from a short derivation, and the approximation B_QS² ≈ 0 should be stated as an explicit assumption before the ΔB calculation, not only in the following paragraph.","section":"Eq. (5)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is technically competent in its radiative-transfer machinery and transparent in its data presentation, but the central claim of a statistically significant preference for the non-hydrostatic model is overstated. The heterogeneity of the dataset and the mislabeled hydrostatic baseline compound the problem. I would ask the authors to either supply a correct significance analysis and a properly formulated hydrostatic model, or substantially soften the conclusions. The paper could then be a useful contribution to solar radio physics."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper has two genuine contributions: a new 2.8 mm ALMA single-dish prominence brightness temperature measurement, and a systematic two-parameter fit of the ATLCW quiet-Sun model to a composite radio dataset. The reported density and temperature factors are in line with earlier prominence work, and the authors are honest about the heterogeneity of the data, including the different instruments, beams, and structures, and the likely beam-convolution effect. The radiative transfer calculation is standard, and Table 1 gives enough information for someone to reproduce the fit. That is real value.\n\nThe soft spots are mostly in the interpretation. The paper claims the non-hydrostatic fit is 'highly statistically significant' over the hydrostatic fit, but the numbers in Table 2 do not support that. With Δχ² = 2.7 (procedure a) and 2.3 (procedure b) for one extra free parameter, the p-values are about 0.10 and 0.13. That is not significant at the 5% level. The authors even note that the non-hydrostatic model has an extra free parameter, so they should know better; this is not a minor wording issue, it is the main support for the magnetic-support conclusion. Relatedly, the 'hydrostatic equilibrium' baseline is not actually a hydrostatic atmosphere. Fixing f_T = 1/f_n enforces constant gas pressure between the PR and the quiet Sun, not dP/dh = -ρg. That is a different physical assumption, and calling it hydrostatic muddies the comparison. A true hydrostatic stratified model would be the right baseline, and the paper should either implement one or clearly state that it is testing constant-pressure versus independent density/temperature.\n\nThe other concern is the dataset. Five measurements from three instruments at different epochs, with beams from 17\" to 144\", are fit as if they sample one homogeneous prominence. The authors acknowledge the issue but do not correct or model the beam convolution, so the fitted factors may be biased toward lower contrast at long wavelengths. This does not destroy the modeling exercise, but it means the quoted density and temperature factors should carry larger effective uncertainties than the nominal parameter ranges in Table 2.\n\nI would not defend the 'confirms thermal bremsstrahlung' statement either: the model assumes bremsstrahlung and fits to the data, so a good fit is not independent confirmation.\n\nWho is this for? Solar radio astronomers working on prominence millimeter emission and semi-empirical modeling. The paper deserves a serious referee because it contains new data and a reproducible analysis, but it needs major revision on the statistics, the hydrostatic baseline, and the strength of the conclusions. I would send it to peer review, not desk reject, and would ask the authors to report a proper significance test (p-value or AIC) and resubmit.","headline":"Useful new ALMA Band 3 prominence measurement and a transparent modeling exercise, but the statistical preference for non-hydrostatic support is overstated and the hydrostatic baseline is mislabeled.","tokens_in":19373,"tokens_out":1639,"would_cite":false,"duration_ms":16400,"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":"The paper claims that a perturbed ATLCW quiet-Sun model with independently varied density and temperature factors, not hydrostatic coupling, explains ALMA and legacy radio measurements of prominence brightness temperatures, implying…","keywords":["solar prominences","brightness temperature","ALMA single-dish","ATLCW model","thermal bremsstrahlung","hydrostatic equilibrium","magnetostatic support","millimeter radio emission"],"falsifier":"Observe a single prominence with interferometric ALMA at arcsecond resolution across multiple bands and measure its electron density independently from spectral lines; if the beam-corrected brightness temperatures no longer prefer free density and temperature factors over the hydrostatic coupling, or the measured density falls outside 15–103 times the quiet Sun, the central claim fails.","tokens_in":18192,"feed_emoji":"☀️","tokens_out":7139,"duration_ms":51285,"temperature":0.7,"pith_summary":"The paper aims to show that a single perturbed quiet-Sun model, the ATLCW atmosphere, can explain millimeter-wave brightness temperatures of solar prominences on the disk, and that the perturbation must treat density and temperature as independent parameters. The best non-hydrostatic fit gives prominence densities 60–68 times the quiet-Sun value and temperatures 155–159 times lower, reproducing the observed absorption at ALMA wavelengths while the hydrostatic fit (densities 159–163 times higher, temperatures the inverse) produces unphysical short-wavelength emission. If the claim holds, prominences are cool, dense plasma supported by magnetic fields rather than gas pressure alone, and their millimeter emission is thermal bremsstrahlung. A reader should care because ALMA single-dish data are sparse but this analysis turns five measurements into quantitative constraints on prominence density, temperature, and stability.","feed_headline":"Radio fits favor magnetic support, not pressure, for prominences","feed_subtitle":"Free density and temperature factors fit the millimeter absorption better than pressure equilibrium.","key_machinery":"The load-bearing object is the modified ATLCW quiet-Sun atmosphere: a 1D semi-empirical model of temperature and density versus height that the paper multiplies by constant factors $f_n$ and $f_T$ in a single 10,000-km layer at 40,000–50,000 km. The calculation then integrates the thermal-bremsstrahlung optical-depth formula (Equation 1) through the whole model atmosphere and converts to brightness temperature via the Rayleigh-Jeans radiative-transfer equation (Equation 2). A $\\chi^2$-minimization over $f_n$ and $f_T$ against five radio measurements yields the best fits. The key identity is the hydrostatic constraint $f_T = 1/f_n$ compared with the non-hydrostatic case where $f_n$ and $f_T$ are free; the paper's evidence for magnetic support is the statistical and physical preference for the free fit.","core_discovery":"On the paper's own terms, the central claim is that the 1D semi-empirical ATLCW quiet-Sun model, with electron density and temperature multiplied by constant factors over the height range 40,000–50,000 km and integrated over 0–57,797 km, reproduces the measured brightness temperatures of prominences on the solar disk at 1.21, 2.80, 3.10, 8.10, and 8.30 mm. The non-hydrostatic fit, in which the density factor (60–68, with uncertainty 15–103) and the temperature factor (1/159 to 1/155) are independent, is statistically preferred over the hydrostatic fit, in which temperature is forced to be the inverse of density (159–163). The paper interprets this as evidence that prominence stability is maintained by magnetic fields obeying magnetostatic conditions rather than by hydrostatic equilibrium, and that thermal bremsstrahlung is the dominant emission mechanism at these wavelengths. The claim includes derived prominence densities of about $0.35$–$3.9\\times10^{10}$ cm$^{-3}$, consistent with earlier optical and EUV estimates, and prominence temperatures of about 6,290–6,957 K, on the cool side of previously reported ranges.","pith_inferences":["If the magnetic-support interpretation is right, the same modeling framework could be inverted: given measured brightness-temperature contrasts, one could map the excess magnetic field required for pressure balance, potentially yielding prominence field estimates from routine ALMA single-dish data without polarimetry.","The weak density dependence of the fitted brightness temperature noted in the paper suggests that millimeter observations primarily constrain prominence temperature, not density; combining ALMA with a density-sensitive diagnostic (e.g., EUV line ratios or H-alpha emission measure) for the same structure could break this degeneracy.","The method's reliance on perturbing a quiet-Sun model in a single layer could be extended to other magnetically supported cool structures, such as filaments in transition or coronal rain, where the same height range and free density and temperature factors may apply.","Future ALMA Band 5 and Band 7 observations, by filling the 3–10 mm gap now bridged by Metsähovi and Nobeyama, would test whether the fitted wavelength trend holds within a single instrument and without cross-calibration differences."],"forward_implications":["Prominences on the solar disk should appear darker than the quiet Sun at all ALMA wavelengths, with contrast growing toward longer wavelengths; the non-hydrostatic model keeps this true even below 2 mm, where the hydrostatic model predicts excess emission.","Electron densities in quiescent prominences at chromospheric heights are about 60–68 times the quiet Sun (up to 103 with uncertainties), and temperatures about 155–159 times lower, giving values around $10^{10}$ cm$^{-3}$ and roughly 6,300–7,000 K.","The statistical preference for the non-hydrostatic fit implies that the global minimum of the brightness-temperature fit is reached with the plasma not in hydrostatic equilibrium, favoring magnetostatic support as the stability mechanism.","Thermal bremsstrahlung, with the same Gaunt-factor treatment, is sufficient to explain millimeter and sub-millimeter prominence emission; no additional emission mechanism is required."],"supporting_citations":[{"why":"Supplies the ATLCW quiet-Sun atmosphere that the paper perturbs.","marker":"Avrett et al. (2015)"},{"why":"Provides the brightness-temperature calculation method and the chi-square fitting procedure used here, validated for quiet Sun, active regions, and coronal holes.","marker":"Matković et al. (2024)"},{"why":"Supplies the ALMA Band 6 (1.21 mm) prominence brightness temperature measurement and the H-alpha boundary extraction method.","marker":"Brajša et al. (2018b)"},{"why":"Supplies the Nobeyama 3.10 and 8.30 mm prominence measurements used in the fitting dataset.","marker":"Hiei et al. (1986)"},{"why":"Supplies the Metsähovi 8.10 mm measurement and the earlier FAL-based prominence modeling approach this work extends.","marker":"Brajša et al. (2009)"},{"why":"Supplies the Gaunt-factor interpolation used in the thermal-bremsstrahlung optical-depth calculation.","marker":"van Hoof et al. (2014)"},{"why":"Provides the comparison ranges for prominence density and temperature against which the derived values are checked.","marker":"Parenti (2014)"},{"why":"Supplies the thermal-bremsstrahlung optical-depth formula and the pressure-balance equation used to estimate the excess magnetic field.","marker":"Benz (2002)"}],"fun_headline_variants":["ALMA data favor magnetic support for solar prominences","Prominences: magnetic fields beat pressure in new ALMA model","Magnetic support wins for solar prominences, ALMA shows","ALMA model: prominences held up by magnetism, not gas pressure","Non-hydrostatic fit: solar prominences lean on magnetic fields"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The five radio measurements come from different prominences observed by different instruments at different epochs and spatial resolutions (beams from 17 to 144 arcsec), yet the fitting treats them as one dataset and applies one pair of multiplicative factors over a fixed height range; the paper itself notes that different structures were observed and that beam convolution is not corrected.","fun_headline_variants_meta":{"raw":{"variants":["ALMA data favor magnetic support for solar prominences","Prominences: magnetic fields beat pressure in new ALMA model","Magnetic support wins for solar prominences, ALMA shows","ALMA model: prominences held up by magnetism, not gas pressure","Non-hydrostatic fit: solar prominences lean on magnetic fields"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000261,"raw_usage":{"total_tokens":1651,"prompt_tokens":1063,"completion_tokens":588,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":679,"completion_tokens_details":{"reasoning_tokens":501}},"tokens_in":679,"tokens_out":588,"duration_ms":4936,"temperature":1.0,"reasoning_tokens":501,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:01:04.341865+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe a single prominence with interferometric ALMA at arcsecond resolution across multiple bands and measure its electron density independently from spectral lines; if the beam-corrected brightness temperatures no longer prefer free density and temperature factors over the hydrostatic coupling, or the measured density falls outside 15–103 times the quiet Sun, the central claim fails.","supporting_citations":[{"cited_title":"(1986), Dark filamentsobservedat8.3mmand3.1mmwavelengths","cited_arxiv_id":null,"evidence_quote":"Supplies the Nobeyama 3.10 and 8.30 mm prominence measurements used in the fitting dataset."},{"cited_title":"(2014),MNRAS, 444, 420-428","cited_arxiv_id":null,"evidence_quote":"Supplies the Gaunt-factor interpolation used in the thermal-bremsstrahlung optical-depth calculation."},{"cited_title":"2002, AstrophysicsandSpaceScienceLibrary.Vol.279, Plasma Astrophysics (2nd ed., Dordrecht: Kluwer)","cited_arxiv_id":null,"evidence_quote":"Supplies the thermal-bremsstrahlung optical-depth formula and the pressure-balance equation used to estimate the excess magnetic field."}],"review_version":2}