{"id":"1d5247f6-c8ce-4a07-b482-085a7fc3dcd7","arxiv_id":"2607.29455","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Spectral ratio analysis of the Fe I 4377 Å line in HARPS-N Sun-as-a-star spectra recovers facular filling factors that track SDO/HMI measurements (Pearson R = 0.587–0.927).","lead":"This paper shows that the Fe I 4377 Å line in disc-integrated sunlight can be used to estimate what fraction of the Sun's surface is covered by bright facular patches. This matters for exoplanet surveys, where faculae create stellar noise that can hide or mimic small planets.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fe I 4377 Å residual is not isolated to faculae: the model assumes a single facular component, yet the inferred filling factor tracks network-inclusive SDO curves at high activity, so the 'facular filling factor' claim is not uniquely established.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: the 4377 Å feature is assumed to be a predominantly facular tracer, but the paper does not demonstrate that spots/network do not contaminate it. My analysis confirms this is the most critical point for the central claim, because the title, abstract, and conclusions explicitly claim recovery of 'facular filling factors' and a 'direct observational pathway' for stellar activity monitoring. If the line is actually a mixed tracer, the quantitative filling-factor estimates are not physically meaningful facular coverages, and transferring the method to other stars (with different network/spot ratios) would introduce unknown biases. The paper's own Figure 9 and Section 5.1 show the inferred coverage shifting between the NT and T SDO curves, which is direct evidence of network sensitivity. Other concerns I considered — e.g., the post-hoc removal of PHOENIX 2022 days, the 100 G model exclusion, and the lack of significance testing on the Pearson R — are secondary; even if the correlations are robust, they validate a mixed-tracer indicator, not a facular-specific one. The proposed concrete test would settle whether spots/network contribute significantly to the 4377 Å response. Until then, the paper should be accepted conditionally, with the claims restricted to 'magnetic bright feature filling factor' or with the additional analysis demonstrating facular specificity. This preserves the reader's CONDITIONAL verdict, so no change is needed.","tokens_in":23057,"tokens_out":7562,"duration_ms":85065,"concrete_test":"Compute daily spot and network filling factors from SDO/HMI for the two epochs (spot: B>24 G and I<0.89*I_quiet; network: NT minus T). Fit a multiple regression of the daily 4377 Å SRA residual amplitude (or the PHOENIX/MURaM-inferred filling factor) on facular filling factor, spot filling factor, and network filling factor. If the spot or network partial coefficients are statistically significant (p<0.05 after accounting for autocorrelation), the line is not a specific facular tracer and the inferred 'facular filling factors' are biased; if insignificant, the facular interpretation is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the 4377 Å SRA residual responds predominantly to photospheric faculae, so that the single-component model (quiet Sun + faculae proxy) yields the facular filling factor. The paper does not establish this. It cites Thompson et al. (2020) for the feature's stronger correlation with facular than spot filling factor, but does not test the line's sensitivity to spots or network in its own data or models. Section 5.1 admits that the inferred filling factor tracks the non-thresholded SDO curve (which includes network) at high activity and the thresholded curve at low activity, implying the tracer's response is a mix of faculae and network (and possibly spots). Since the synthetic fitting attributes all residual amplitude to a single facular component, the derived filling factor is not physically the facular filling factor; it is an effective coverage of whatever magnetic bright features the line actually sees. For Sun-as-a-star observations this may still produce a useful activity indicator, but the quantitative claim of recovering 'facular filling factors' is unsupported, and application to other stars with different network/spot mixes would be biased.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses Spectral Ratio Analysis (SRA) on HARPS-N Sun-as-a-star spectra from two ~3-rotation windows (2016 June–August and 2022 May–July) to isolate the Fe I 4377 Å residual feature and model it with composite PHOENIX spectra (facular temperature contrast 200–400 K) and MPS-ATLAS spectra from MURaM simulations with imposed vertical magnetic fields of 100, 200, and 300 G. The inferred facular filling factors are compared with SDO/HMI filling factors, with reported Pearson R values of 0.587–0.927 depending on model, epoch, and SDO definition (thresholded vs non-thresholded). The paper concludes that the Fe I 4377 Å line can trace relative variations in facular filling factors, offering a photospheric activity diagnostic for RV exoplanet work. The workflow is clearly described, the external SDO comparison is a strength, and the limitations are candidly discussed; however, the central quantitative claim that the fitted single-component coverage is specifically the facular filling factor is weakened by admitted mixing with the magnetic network and by post-hoc data selection in the PHOENIX correlations.","tokens_in":23323,"tokens_out":3695,"duration_ms":44572,"significance":"If the central claim holds, the paper provides a direct, disc-integrated spectroscopic tracer of photospheric facular coverage that could be applied to other stars and improve activity mitigation in RV surveys and transmission-spectroscopy analyses. The manuscript has genuine strengths: it uses public HARPS-N and SDO data, documents the SRA and model-fitting pipeline in detail, validates against an independent disc-resolved dataset, and uses two complementary synthetic frameworks (1D PHOENIX and 3D MURaM). The clear statement of caveats, including the temperature–filling-factor degeneracy and the limb-angle insensitivity, is commendable. However, the quantitative identification of the inferred coverage as the facular filling factor is not uniquely established, because the model attributes all residual amplitude to a single facular-like component while the comparison in Section 5.1 shows that the inferred quantity tracks network-inclusive and area-thresholded SDO curves differently at different activity levels. This weakens the main quantitative claim, though the qualitative conclusion that the line is a useful activity proxy may still survive with reframing or additional modelling.","major_comments":[{"comment":"The claim that the fitted alpha is specifically the facular filling factor is not uniquely established. Section 5.1 states that in the high-activity period the inferred coverage mainly tracks the non-thresholded SDO curve (which includes the magnetic network), while in the low-activity period it tracks the area-thresholded curve. Because the model contains only one active component, alpha is an effective coverage of whatever bright magnetic features the line responds to, not necessarily photospheric faculae alone. The paper should either include a network/spot component (or otherwise demonstrate insensitivity to such contamination) or reframe the central claim as measuring the coverage of magnetically bright photospheric features. This is load-bearing for the abstract and conclusions, which state specifically that facular filling factors are recovered.","section":"§5.1, Fig. 9"},{"comment":"The reported improvement in the PHOENIX 2022 correlations from R=0.640 to R=0.795 after removing several days is not reproducible as written. The manuscript does not state how many points were removed, which days, or the selection criterion beyond 'larger errors bars ... pulling the reduced chi2 fit into the continuum.' Since these correlations are a key validation of the method, the removal must be pre-specified or fully transparent. Please provide the removed dates, the number of points, and an analysis showing the result is not driven by post-hoc selection. Also report p-values or confidence intervals for all correlation coefficients, given the modest number of days.","section":"§5.1, Table 1"},{"comment":"The temperature–filling-factor degeneracy is acknowledged in Section 5.2, but the reported alpha values are presented without uncertainties. The reduced-chi2 maps in Figure B1 show broad, correlated minima, so the daily alpha estimates in Figure 9 likely have substantial uncertainties that are not propagated into the SDO correlations. Please provide confidence intervals or posterior bounds for the best-fit alpha (and for the MURaM fits), and account for this uncertainty when quoting Pearson R values.","section":"§4.1, Fig. B1"},{"comment":"The selection of the 200 G MURaM model as the adopted facular proxy is not fully supported by the model comparison. The text notes that the 300 G models give the lowest reduced chi2 but that 'little difference' is seen between 200 G and 300 G, and 100 G models were removed post hoc because they underestimate the amplitude. The conclusion that 'the best estimates of faculae filling factor are given by magnetic field strengths of approximately 200G' is therefore based on an external canonical value rather than on the data. Please either present the comparison as a choice with motivation (e.g., Schrijver & Harvey 1994) rather than a fit result, or provide model-selection statistics (e.g., AIC/BIC or chi2 differences) that justify preferring 200 G over 300 G.","section":"§4.2, Fig. 6"}],"minor_comments":[{"comment":"The abstract states that SRA offers a means to 'reliably track surface magnetic activity in disc-resolved spectra.' The observations are disc-integrated Sun-as-a-star spectra; 'disc-resolved' should be 'disc-integrated' (or the sentence should be reworded) to avoid confusion with SDO-style resolved imaging.","section":"Abstract"},{"comment":"The ripple model uses a Lomb-Scargle periodogram to identify the frequency f, but the text does not specify the wavelength range over which the periodogram is computed or the units of f. Please clarify, since the wavelength-scaled term lambda_med/lambda changes the effective frequency.","section":"Eq. (1)"},{"comment":"Typo: 'further investigiation' should be 'further investigation.' Also, 'MP-ATLASMURaMmodels' in the same section appears to be a spacing error for 'MPS-ATLAS MURaM models.'","section":"§3"},{"comment":"The table is referred to in the text as 'Table 5.1'; it should be 'Table 1.' The caption should also state the number of days in each period and the uncertainty on R.","section":"§5.1, Table 1"}],"recommendation":"major_revision","confidential_remarks":"The paper is appropriate for the journal and the topic is timely. The main issue is that the quantitative claim of measuring 'facular filling factors' is broader than what the single-component model can support, given the admitted network contribution and the post-hoc removal of outlier days. I would not reject the paper: the core idea and the external SDO comparison are valuable, and the requested changes (multi-component test or reframing, transparent outlier handling, uncertainty propagation, and model-selection justification for 200 G) are feasible within the scope of a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid empirical methods paper with a clearly described pipeline and an external SDO check. What is new is not SRA itself—Thompson et al. already established the 4377 Å residual and its correlation with facular filling factors, and Costes et al. used SRA for RV indicators—but the quantitative forward-model inversion: fitting PHOENIX and MURaM/MPS-ATLAS synthetic SRA residuals to HARPS-N Sun-as-a-star spectra and recovering a filling-factor time series without fitting to SDO as the target. That is genuinely useful, and the MURaM correlations (R ≈ 0.82–0.93) are impressive for a single line. The paper also provides data links and is unusually transparent about its limitations, which earns credit.\n\nWhere I would push back: the facular-specific reading is not uniquely pinned down. The model is quiet Sun plus one facular component, and everything non-quiet is assigned to faculae. Section 5.1 admits the inferred coverage tracks the network-inclusive NT curve at high activity and the thresholded T curve at low activity. That means the fitted alpha is an effective coverage of whichever bright magnetic features the line sees, not necessarily photospheric faculae. The Sol'Ex context image and the Thompson et al. correlation are suggestive, but the paper does not test spot or network sensitivity of 4377 Å in its own models. The quantitative \"facular filling factor\" claim is therefore overreach; the relative activity-indicator claim is fine.\n\nOther soft spots: the PHOENIX 2022 improvement from R = 0.640 to 0.795 comes after removing points post hoc; that should be presented as a sensitivity test, not the headline. The low-activity template is set to 1% faculae using SDO, so absolute calibration is partly anchored to the benchmark the method is validated against. The temperature–filling-factor degeneracy (Figure B1) means absolute values are model-dependent. And the abstract says \"disc-resolved\" where it should say \"disc-integrated\"—that is a factual slip that will confuse readers.\n\nBottom line: if you take the paper as establishing a relative, rotationally-modulated photospheric activity indicator from one Fe I line, it holds up. If you take it as measuring true facular filling factors, it does not yet. It deserves a serious referee; revisions should focus on renaming or reframing the derived quantity, adding a spot/network sensitivity test, and treating the 2022 exclusions transparently.","headline":"A solid, honest SRA-to-filling-factor inversion for the Sun—the MURaM correlations are convincing for a relative activity indicator, but the facular-specific quantitative claim is not uniquely established and needs a network/spot sensitivity test.","tokens_in":24231,"tokens_out":3048,"would_cite":true,"duration_ms":35663,"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":"The Fe I 4377 Å line reveals facular coverage from disc-integrated sunlight.","keywords":["spectral ratio analysis","faculae","solar activity","Sun-as-a-star","radial velocity noise","photospheric activity indicator","Fe I 4377 Å"],"falsifier":"A decisive test is to measure the 4377 Å SRA residual during a low-activity period when the solar disc contains magnetic network but essentially no active-region faculae. If the inferred filling factor rises with the network coverage (as measured by non-thresholded magnetogram pixels) while area-thresholded facular coverage stays near zero, the line is not a clean facular tracer and the inversion would need a network component.","tokens_in":22874,"feed_emoji":"☀️","tokens_out":6495,"duration_ms":61656,"temperature":0.7,"pith_summary":"Faculae are bright magnetic patches that dominate stellar activity noise in radial-velocity searches for Earth-like exoplanets, but their low contrast makes them hard to track in disc-integrated light. This paper argues that a single photospheric iron line, Fe I 4377 Å, carries a strong and coherent facular signal when high-activity spectra are divided by a low-activity template (Spectral Ratio Analysis). By fitting the residual feature with composite synthetic spectra — either temperature-perturbed or magnetically enhanced — the paper recovers fractional facular coverage whose time evolution matches disc-resolved solar images (Pearson R ≈ 0.59–0.93) over two multi-rotation epochs. If the claim holds, this provides a photospheric activity indicator that directly probes the regions responsible for RV noise, complementing chromospheric proxies like log R'HK.","feed_headline":"One iron line tracks the Sun's facular coverage","feed_subtitle":"Fitting the 4377 Å spectral ratio recovers facular coverage matching disc-resolved images — a photospheric activity probe for exoplanet sear","key_machinery":"The central object is the Fe I 4377 Å absorption line and its SRA residual, constructed by dividing a daily high-activity disc-integrated spectrum by a low-activity template. The line's sensitivity to faculae follows from atomic physics: in the quiet photosphere iron is mostly singly ionized, so neutral Fe I is a minority species; the facular temperature enhancement ionizes more Fe I, depletes the ground state (the line's lower level, E_low = 0 eV), and weakens the line. The carrying mechanism of the argument is a forward spectral model: a linear combination of quiet-Sun and facular templates — one family using hotter atmospheres, the other using magnetized convection simulations — is conver","core_discovery":"The paper's central claim is that the Fe I 4377 Å line's Spectral Ratio Analysis residual is dominated by faculae, and that its amplitude can be inverted, using forward models of composite spectra, to estimate the fractional area of the solar disc covered by faculae. Two complementary synthetic frameworks are used: one approximates faculae as LTE atmospheres 200–400 K hotter than the quiet Sun; the other uses radiation-MHD simulations of magnetic regions with ~200 G mean vertical fields. Inferred filling factors rise from ~1.5% to ~5.5% across the two studied activity states and track disc-resolved filling factors with Pearson R = 0.587–0.927. The paper concludes that the line can trace rela","pith_inferences":["A multi-line extension of the SRA fit, using several iron lines with different sensitivities, could break the degeneracy between filling factor and temperature contrast and recover limb-angle information that the single-line fit cannot.","If the facular-only interpretation survives on the Sun, the same line could serve as a spectral template for correcting transmission spectroscopy of exoplanets, where unocculted faculae bias inferred planetary radii at short wavelengths.","On more active stars, spot contamination could bias the filling-factor inversion; calibrating the method against solar epochs with substantial spot coverage would quantify that bias.","The systematic offset between inferred and thresholded coverage at low activity suggests the line may be a tunable tracer of the network-to-faculae ratio, a parameter otherwise hard to observe in disc-integrated light."],"forward_implications":["Disc-integrated monitoring of the 4377 Å line can recover the rotational modulation and activity-cycle evolution of facular coverage without spatially resolved imaging.","The inferred filling factor provides a photospheric, line-profile-based activity indicator directly tied to the regions responsible for radial-velocity noise, potentially improving RV decorrelation in exoplanet surveys.","The same fitting framework can in principle be applied to high-resolution spectra of other Sun-like stars, giving a stellar facular filling factor without resolved disc imaging.","The recovered facular temperature contrasts of roughly 200–400 K are consistent with the range used in solar activity modelling, supporting the interpretation of the residual as a temperature-sensitive photospheric signal.","The inferred coverage tracks the non-thresholded filling factor (faculae plus network) at high activity and the area-thresholded one at low activity, indicating where the line's sensitivity begins."],"fun_headline_variants":["Iron line at 4377 Å reveals Sun's facular coverage","Single spectral line measures solar faculae like disk images","Fe I 4377 line: a new proxy for solar magnetic activity","Spectral ratio of one iron line tracks faculae over cycles"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the 4377 Å SRA residual is produced almost entirely by faculae, so that its amplitude can be converted directly into a facular filling factor; if spots or the magnetic network contribute significantly, every derived filling factor is biased even if the correlations remain high.","fun_headline_variants_meta":{"raw":{"variants":["Iron line at 4377 Å reveals Sun's facular coverage","Single spectral line measures solar faculae like disk images","Fe I 4377 line: a new proxy for solar magnetic activity","Spectral ratio of one iron line tracks faculae over cycles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000193,"raw_usage":{"total_tokens":1254,"prompt_tokens":879,"completion_tokens":375,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":623,"completion_tokens_details":{"reasoning_tokens":303}},"tokens_in":623,"tokens_out":375,"duration_ms":4941,"temperature":1.0,"reasoning_tokens":303,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T06:40:33.224760+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test is to measure the 4377 Å SRA residual during a low-activity period when the solar disc contains magnetic network but essentially no active-region faculae. If the inferred filling factor rises with the network coverage (as measured by non-thresholded magnetogram pixels) while area-thresholded facular coverage stays near zero, the line is not a clean facular tracer and the inversion would need a network component.","supporting_citations":[],"review_version":1}