{"id":"ce401d54-3ad3-417e-b1b4-de8217ab1b5a","arxiv_id":"2608.06306","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Fixing the stellar inclination angle in spectropolarimetric inversions can produce visually perfect Stokes fits with substantially wrong magnetic field parameters, as shown by large synthetic tests with a neural-network surrogate.","lead":"This paper tests how accurately stellar magnetic fields can be recovered from spectropolarimetric data when the star's rotation angle is assumed instead of measured, using a fast neural-network model of the Stokes profiles. It finds that fixing the rotation angle to a wrong value can still produce excellent-looking fits while severely biasing the inferred magnetic field, a degeneracy that affects existing magnetic-mapping studies.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unverified hinge: 'satisfactory fits' under wrong inclination are scored against MAPNet's own error-normalized metric (App. B), never re-checked with cossam; if MAPNet's ~7% Q/U errors absorb inclination differences, the claimed degeneracy could be partly a surrogate artifact.","rationale":"The paper's contribution is a numerical demonstration that fixing the stellar inclination angle in spectropolarimetric inversions can still produce excellent Stokes fits while biasing the recovered magnetic parameters. For that demonstration to support the stated conclusion, the word 'satisfactorily' must mean satisfactory with respect to the true forward model, cossam, since the observed profiles are generated by cossam. Instead, the fitness function EN-WMAPE normalizes each Stokes WMAPE by the ANN's expected synthesis error, estimated from MAPNet's own validation data. This is a reasonable device for balancing Stokes components during inversion, but it converts 'agreement at the level of MAPNet's typical error' into an acceptable fit. Appendix A shows that low-amplitude Q and U profiles, the very profiles most affected by inclination-dependent projection, carry the largest ANN errors, around 7% WMAPE. In the fixed-inclination experiments, the optimizer is therefore allowed to stop at models whose polarized profiles differ from the observed cossam profiles by roughly the ANN's error floor, and no independent check establishes that those models would also be accepted by cossam. The paper does validate MAPNet's overall synthesis accuracy (Fig. 1), and the inversion tests are internally consistent, but that does not close the gap: a surrogate can be globally accurate yet locally smooth over the specific inclination-induced differences that the degeneracy claim depends on. The missing cossam-level re-synthesis is a concrete, feasible check, and it is the natural condition under which the paper's central claim should be accepted. This is the same concern the reader identified, and I agree with the CONDITIONAL verdict: the qualitative conclusion is plausible and important, but it is not yet fully demonstrated at the true forward-model level. No code, weights, or training data are released, which further limits independent verification, although that is secondary to the missing cossam check.","tokens_in":13551,"tokens_out":3893,"duration_ms":46064,"concrete_test":"Take a random subset (≥50, including the Fig. 5 and Fig. D1 examples) of the fixed-inclination best-fit parameter vectors from Tables 4/5. For each, synthesize all seven phases with cossam and compare to the original cossam 'observed' profiles using raw WMAPE and reduced chi2 (no EN-WMAPE normalization). Also run ~10 fixed-inclination inversions using cossam directly as the forward model for a subsample and check whether the same wrong-inclination family of solutions appears. If cossam-level errors for wrong-inclination solutions are similar to free-inclination solutions (Table 3; e.g., Q/U WMAPE below a few percent), the physical degeneracy holds; if they exceed MAPNet's median errors or fail a chi2 threshold, the ANN surrogate manufactured the degeneracy.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that, for noise-free cossam-generated Stokes profiles, fixing a wrong inclination still yields satisfactory fits, so the inversion is degenerate. The fitness function used to define 'satisfactory' is EN-WMAPE (Eq. B1), which divides each Stokes WMAPE by the ANN's empirically fitted expected error \\WMAPE_s(A). The coefficients come from MAPNet's own test data, and Appendix A shows Q/U WMAPE reaches roughly 7% for low amplitudes. Thus a candidate profile that differs from the observed cossam profile by about the ANN's typical synthesis error is treated as a good fit. Since the wrong-inclination experiments are exactly the regime where polarized profiles can be weak, the normalization can make physically distinct profiles look equivalent. The paper never re-synthesizes the wrong-inclination best-fit parameter vectors with cossam and compares them with the original cossam profiles using a model-independent metric. Consequently, a key part of the demonstrated degeneracy—that the same profiles can be fitted by wrong-inclination models at the true forward-model level—remains unverified. If MAPNet's interpolation smooths over inclination-dependent profile differences, the degeneracy could be substantially a surrogate artifact rather than a property of the physical model. This is the load-bearing concern because it targets the definition of 'satisfactorily fitted' in the paper's strongest claim and is not settled by the ANN validation plots in Fig. 1, which measure MAPNet against cossam but not the fixed-inclination solutions against cossam.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces MAPNet, an artificial neural network surrogate for the cossam polarized radiative transfer code, and combines it with particle swarm optimization to invert synthetic Stokes profiles for a set of 12 magneto-atmospheric parameters. Using single-line profiles, the inversion recovers magnetic parameters well but struggles with Teff and log g; using multi-line (SVD) profiles, all parameters are recovered with R^2 > 0.9 when the inclination angle is free. The central claim is that fixing the rotational inclination angle to an incorrect value still allows the observed multi-line profiles to be satisfactorily fitted, while the recovered magnetic moment and dipole position are seriously biased, demonstrating a degeneracy in the inversion problem. The paper also discusses the preference for SVD over LSD line addition when atmospheric parameters are free, and it reports a large computational campaign of over 6,600 inversion cases.","tokens_in":13817,"tokens_out":7990,"duration_ms":80606,"significance":"If the degeneracy claim holds at the true forward-model level, this is an important and cautionary result for stellar spectropolarimetric inversions, since many mapping codes fix the inclination angle a priori. The experimental design is a strength: the observed profiles are generated with the external cossam code, the ground-truth parameters are known, and the free-inclination multi-line inversion results (Table 3) are convincing evidence of the method's basic validity. The paper also demonstrates a practical technique (SVD line addition with free atmospheric parameters) and a computationally efficient GPU-based inversion strategy. However, the central degeneracy claim is evaluated with a metric that is normalized by the surrogate's own error statistics, and the wrong-inclination solutions are never checked against the true radiative transfer code; this leaves the main conclusion less firmly established than the abstract suggests.","major_comments":[{"comment":"The claim that the observed profiles 'can be satisfactorily fitted regardless of the assumed stellar inclination angle' is evaluated with the EN-WMAPE fitness metric, whose amplitude-dependent normalization coefficients are fitted to MAPNet's test data. The wrong-inclination best-fit solutions are never re-synthesized with cossam and compared with the original observed profiles using a model-independent metric. Given that MAPNet's median WMAPE reaches about 7% for low-amplitude Q and U profiles (Appendix A), the surrogate could smooth over inclination-dependent profile differences, making the degeneracy partly an artifact of the surrogate rather than a property of the physical forward model. I request that the authors take a sample of the fixed-inclination solutions (including the two shown in Fig. 5) and recompute their Stokes profiles with cossam, then compare them to the observed cossam profiles using a metric not normalized by the ANN's error statistics. This verification is necessary to support the paper's central claim.","section":"Section 4, Fig. 5, Appendix B, Eq. (B1)"},{"comment":"The statement that the wrong-inclination inversions produce 'satisfactory' fits is supported only by two illustrative examples (Fig. 5 and Fig. D1). The manuscript does not report the distribution of the fitness values (EN-WMAPE) for the fixed-inclination runs, nor does it compare them with the distribution for the free-inclination runs. If the fixed-inclination fits are systematically worse than the free-inclination fits, the word 'satisfactorily' would require qualification, and the severity of the degeneracy would need to be assessed differently. The authors should provide a statistical summary, for instance the median and a percentile range of the EN-WMAPE over the full sample for each configuration, and state what fraction of the fixed-inclination inversions achieve a fit quality comparable to the free-inclination median.","section":"Section 4, Tables 4 and 5"}],"minor_comments":[{"comment":"The abstract states that fixing 'the atmospheric parameters and the stellar inclination angle' is shown to induce very high errors, but the experiments in Section 4 only fix the inclination angle; the atmospheric parameters remain free. Please revise the abstract and the opening of Section 4 to match the actual scope of the experiment.","section":"Abstract"},{"comment":"The sentence 'only the metallicity and the vsini parameters (in the case of four Stokes profiles) preserves R2 > 0.9' is contradicted by Table 4, where vsini has R2 = 0.801. Please correct the text or the table.","section":"Section 4, paragraph after Tables 4 and 5"},{"comment":"The denominator of WMAPE, the sum of absolute observed values, can be zero for Stokes Q and U in continuum regions, making the metric undefined. Please clarify the wavelength range over which the sum is taken (e.g., only over the line region) or introduce a small regularization floor in the denominator.","section":"Appendix B, Eq. (B2)"},{"comment":"The phrase 'we repeated the precendet test' contains a typo; it should be 'precedent test'.","section":"Section 4, first paragraph"},{"comment":"The paper mentions 1,298 'observed' stars for the single-line test and 1,008 for the multi-line test, but the total number of inversion cases is given as 6,628 in the conclusions. Please provide a breakdown of how this total is composed (including the fixed-inclination runs and possibly other configurations).","section":"Section 3, inversion samples"},{"comment":"The caption for Fig. 2 says 'percentile 75 arranged in order of their fitness'; it would be clearer if the selection criterion for the displayed case were spelled out (e.g., the case with the 75th-percentile fitness value). The same applies to Fig. C1.","section":"Figure 2 and Appendix C"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about the surrogate-based fitness metric is well-founded and should be addressed in the revision. The requested cossam re-synthesis of a subsample of wrong-inclination solutions is a concrete, feasible check that would either support or substantially weaken the paper's central claim. The abstract's mention of fixed atmospheric parameters overstates the experimental coverage and should be corrected regardless. Overall, the paper presents a valuable methodology and a potentially important result, but the central degeneracy claim needs additional verification before it can be fully accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper delivers a real result with one hinge that needs tightening. It shows, over more than a thousand synthetic stars, that fixing the inclination angle to a wrong value can bias the recovered dipole offset and moment severely while the Stokes fits still look good. That quantitative demonstration is new and worth taking seriously.\n\nThe good stuff first. The study design is honest: ground-truth parameters are known, profiles come from the external cossam code, and the inversion uses a fast ANN surrogate (MAPNet) validated against cossam. The multi-line all-free inversions recovering all parameters with R2 > 0.9 is a genuinely useful sanity check, and the large-scale search is a legitimate use of the surrogate. The paper also scopes itself clearly to a de-centered dipole and noise-free data, and flags the higher-order generalization as pending.\n\nThe soft spot is the one the stress-test flags, and it is real. The fitness metric, EN-WMAPE, normalizes each Stokes WMAPE by the ANN's empirically fitted expected error. A candidate profile that disagrees with the observed cossam profile by roughly the ANN's typical synthesis error is therefore treated as a good fit. MAPNet's Q/U errors reach about 7% for low amplitudes, which is exactly the regime that matters when wrong inclination makes the polarized profiles weak. The paper never re-synthesizes the wrong-inclination best-fit solutions with cossam and compares those with the original profiles using a model-independent metric. So the claim that the same profiles can be 'satisfactorily fitted' under a wrong inclination is not yet pinned down at the true forward-model level. If MAPNet smooths over inclination-dependent differences, part of the degeneracy could be a surrogate artifact rather than a property of the physical model.\n\nThat said, the tables tell a consistent story even without the examples. The R2 values for the dipole position collapse when inclination is fixed, and the magnetic moment errors grow by an order of magnitude. Those aggregate results do not depend on the visual quality of any particular fit. So I don't think the qualitative conclusion collapses; it just needs to be stated more cautiously, or better, backed by a cossam cross-check.\n\nWho gets value from this: anyone using Zeeman-Doppler imaging codes that fix inclination (INVERS, iMap, PIMMS, TIMES). The paper is a caution flag worth referee time. The main fix is straightforward: take a sample of the wrong-inclination best-fit parameter vectors, synthesize them with cossam, and show they still fit the observed profiles. If that works, the paper becomes much stronger. If it doesn't, they need to soften the 'satisfactorily fitted' language.\n\nRecommendation: send it to review, but require the cossam cross-check as a major revision.","headline":"Fixed-inclination inversions can badly bias dipole parameters even with excellent-looking fits, but the fit-quality metric is calibrated on the ANN's own errors and never cross-checked with the forward code.","tokens_in":14429,"tokens_out":2880,"would_cite":true,"duration_ms":30283,"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":"With the inclination angle fixed to a wrong value, Stokes profiles from a star can still be fitted satisfactorily while the recovered magnetic moment and dipole position are wrong by hundreds of gauss and up to a tenth of a stellar radius.","keywords":["stellar magnetic fields","spectropolarimetry","Stokes profiles","magnetic field inversion","degeneracy","neural network surrogate","inclination angle","particle swarm optimization"],"falsifier":"Take the best-fit solutions from the fixed-inclination inversions (for example the cases shown in Fig. 5) and synthesize their Stokes profiles with the original polarized radiative transfer code used to create the training data. Compute the same EN-WMAPE or a standard chi-square against the observed profiles; if the wrong-inclination solutions fail this ground-truth comparison while the true-inclination solutions pass, the degeneracy is at least partly an artifact of the neural-network surrogate.","tokens_in":13335,"feed_emoji":"🧲","tokens_out":8668,"duration_ms":73986,"temperature":0.7,"pith_summary":"The paper sets out to show how accurately stellar magnetic fields can be recovered from spectropolarimetric data when atmospheric parameters are allowed to vary, and what happens when the common practice of fixing the stellar inclination angle is followed. Using a fast neural-network surrogate for polarized radiative transfer, the authors invert noise-free synthetic Stokes profiles from a single spectral line and from multi-line averages. They find that multi-line data let all twelve magnetic and atmospheric parameters be recovered with high accuracy when inclination is free. But when inclination is fixed to a wrong value, the Stokes profiles can still be fitted almost perfectly while the inferred magnetic moment and dipole position are off by hundreds of gauss and up to a tenth of a stellar radius. The paper's central claim is that this degeneracy makes visual agreement between observed and fitted profiles an unreliable validation of a magnetic map.","feed_headline":"Fixing the inclination angle skews stellar magnetic maps","feed_subtitle":"Stokes profiles match at any assumed inclination while the recovered dipole strength and position go badly wrong.","key_machinery":"The argument is carried by MAPNet, a fully connected neural network (seven hidden layers of 4,096 neurons) that synthesizes the four Stokes profiles from twelve magneto-atmospheric parameters, trained on 1.5 million profiles computed with a polarized radiative transfer code. Because the network evaluates profiles in milliseconds, a particle-swarm optimizer with 2,048 particles can explore 102,400 candidate parameter combinations per inversion run in about 9 seconds on a GPU. For multi-line data the paper uses singular-value-decomposition (SVD) mean profiles instead of least-squares deconvolution, since SVD does not require prior knowledge of line depths and thus leaves temperature, gravity, and metallicity free. The fitness function EN-WMAPE normalizes each profile's weighted mean absolute percentage error by the network's empirically fitted amplitude-dependent synthesis error, so that weak Q and U profiles do not dominate the optimization.","core_discovery":"The paper demonstrates that, for a star observed at seven rotation phases with noise-free multi-line Stokes profiles (I, V, and optionally Q, U), a de-centered dipole field can be fitted satisfactorily for any assumed inclination angle of the rotation axis. The cost is that the magnetic solution itself is no longer trustworthy: in one example, assuming an inclination wrong by 93 degrees overestimates the dipole magnetic moment by 2,200 gauss (270%) and moves the dipole position outward by 0.1 stellar radii; in another, a 17.5-degree error underestimates the moment by 1,246 gauss (443%). With all four Stokes parameters, a 26-degree inclination error overestimates the moment by 1,168 gauss (40%). Quantitatively, fixing the inclination destroys the recovery of the dipole position ($R^{2}$ < 0.12) and raises the magnetic-moment RMSE from below 50 gauss to above 500 gauss, while Teff, log g, and vsini degrade moderately.","pith_inferences":["A direct check on the physical reality of the degeneracy would be to re-synthesize the best-fit wrong-inclination solutions with the full radiative transfer code; if those solutions no longer match the observed profiles, part of the effect is a surrogate artifact rather than a property of the inverse problem.","Because the EN-WMAPE fitness is calibrated to the network's own amplitude-dependent errors, the inversion's definition of a 'good fit' is relative to the surrogate; the reported magnitudes of the degeneracy should be read with that calibration in mind.","If the degeneracy holds under ground-truth synthesis, the safe practice becomes treating inclination as a free but prior-constrained parameter, and re-interpreting published maps from fixed-inclination inversions as envelopes rather than point estimates.","The anticipated higher-order multipole degeneracy could be tested directly by repeating the fixed-inclination experiment with a dipole-plus-quadrupole field and comparing how many (l,m) coefficient combinations fit the same phase-resolved Stokes series."],"forward_implications":["Published stellar magnetic maps that fixed the inclination angle a priori carry unquantified errors in the dipole moment and especially in the dipole position; with a wrong inclination the position cannot be recovered at all (R^2 < 0.12).","Real observations include noise, and the paper notes that noise would likely enlarge the degeneracy, so the problem is expected to be worse, not better, in practice.","Using the full Stokes vector (I,Q,U,V) instead of only (I,V) reduces the magnetic-moment error (40% at 26 degrees versus 270% at 93 degrees) but does not eliminate it.","Metallicity and vsini remain reliably recovered even with a fixed wrong inclination, while Teff and log g degrade only moderately; the degeneracy is concentrated in the magnetic geometry.","Since the effect appears already for the simplest de-centered dipole, the authors anticipate similar or stronger degeneracies among higher-order spherical harmonic coefficients, a claim left for future work."],"supporting_citations":[{"why":"Paper I; supplies the ANN synthesis architecture and training methodology that MAPNet extends to 12 parameters.","marker":"Raygoza-Romero et al. (2025)"},{"why":"Provides the cossam radiative transfer code that generates the ground-truth training and test profiles.","marker":"Stift et al. (2012)"},{"why":"The particle swarm optimization algorithm used for the inversion search.","marker":"Kennedy & Eberhart (1995)"},{"why":"The SVD line-addition method that builds multi-line profiles without fixing atmospheric parameters.","marker":"Martínez González et al. (2008)"},{"why":"The widely used formalism that motivates fixing the inclination angle a priori, which the paper challenges.","marker":"Donati & Brown (1997)"}],"fun_headline_variants":["Inclination angle choice ruins stellar magnetic inversion","Don't fix star tilt: magnetic field recovery fails","Assuming wrong tilt breaks magnetic mapping","Stellar magnetograms unreliable if tilt is presumed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the neural network's synthesis errors are small enough, and its amplitude-dependent error normalization accurate enough, that the fitness landscape the optimizer sees resembles the true radiative-transfer landscape; if the network smooths away inclination-dependent differences in the weak Q and U profiles, the claimed degeneracy could be partly an artifact of the surrogate.","fun_headline_variants_meta":{"raw":{"variants":["Inclination angle choice ruins stellar magnetic inversion","Don't fix star tilt: magnetic field recovery fails","Assuming wrong tilt breaks magnetic mapping","Stellar magnetograms unreliable if tilt is presumed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000786,"raw_usage":{"total_tokens":3492,"prompt_tokens":996,"completion_tokens":2496,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":612,"completion_tokens_details":{"reasoning_tokens":2439}},"tokens_in":612,"tokens_out":2496,"duration_ms":19443,"temperature":1.0,"reasoning_tokens":2439,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:49:58.427475+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the best-fit solutions from the fixed-inclination inversions (for example the cases shown in Fig. 5) and synthesize their Stokes profiles with the original polarized radiative transfer code used to create the training data. Compute the same EN-WMAPE or a standard chi-square against the observed profiles; if the wrong-inclination solutions fail this ground-truth comparison while the true-inclination solutions pass, the degeneracy is at least partly an artifact of the neural-network surrogate.","supporting_citations":[],"review_version":1}