{"id":"4d506dc3-a160-4b85-a799-accbaa70078d","arxiv_id":"2505.05962","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"P-CORONA computes Stokes I, Q, U, V for forbidden and permitted coronal lines in 3D MHD models with full Hanle and Zeeman treatment.","lead":"A new open-source code, P-CORONA, calculates the intensity and polarization of solar coronal spectral lines inside 3D magnetohydrodynamic models, including Hanle and Zeeman magnetic effects. It is designed to help interpret polarized observations from DKIST and Aditya-L1 and to test coronal magnetic field models.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 1's domain of validity—non-dipole transitions and arbitrary B orientation—is the load-bearing assumption; the paper's own caveats and missing quantitative benchmark leave it unverified.","rationale":"The reader's weakest_assumption was the validity of applying LL04 dipole-approximation theory to magnetic-dipole and non-dipole transitions under the strong-coupling approximation. I share that concern, which is explicitly flagged in the manuscript's own discussion of Casini et al. (2024) as future work. I would add that the same Eq. 1 raises an even more universal issue: the Zeeman profile labels M along the solar-radius quantization axis, but this is only the correct Zeeman eigenbasis when the magnetic field is parallel to that axis. Because no rotation to the magnetic-field frame is described, the central Section 4 claim that the code fully accounts for Zeeman effects on Q and U is unsupported for arbitrary field orientations, which are the norm in 3D MHD models. The missing quantitative pyCELP benchmark matters for the same central claim because the one discrepancy acknowledged is precisely in the Zeeman Q,U terms that are advertised as the new capability. These concerns do not prove the code is wrong; they identify a load-bearing condition that must be demonstrated. They also reinforce the reader's CONDITIONAL verdict, so I recommend keeping that verdict unchanged rather than moving to ACCEPT or REJECT. The proposed rotation-invariance test and independent Zeeman diagonalization would settle whether the emissivity is evaluated in the correct frame, and the E2-admixed line case would test whether the rank-1 formula is being applied outside its domain.","tokens_in":21299,"tokens_out":22928,"duration_ms":271346,"concrete_test":"Run P-CORONA on a single grid cell for Fe xiii 10747 Å with fixed |B| and LOS, with the magnetic field oriented first along the local vertical and then in a perpendicular direction obtained by rigidly rotating the whole (B, LOS, local-vertical) system; the Stokes Q and U profiles must match after rotating the reference direction for positive Q by the same angle. Repeat for a case where B is not aligned with the local vertical and for the E2-admixed Fe xiii 10798 Å line. If the rotated-case profiles do not coincide, or if they disagree with an independent diagonalization of the Zeeman Hamiltonian for a two-level atom, then Eq. 1 is not being evaluated in a magnetic-field-aligned frame and the advertised transverse-Zeeman Q,U signal is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that P-CORONA 'fully accounts' for Hanle and Zeeman effects in forbidden and permitted lines requires Eq. 1 to be the correct emissivity for every transition and field geometry used. Section 2.2 states that LL04 theory is valid 'under the dipole approximation,' yet Eq. 1 is applied to magnetic-dipole and non-dipole transitions satisfying ΔJ=0,±1, while other transitions are reduced to an unpolarized 'strong coupling' population term. The paper itself notes that the multipole extension (Casini et al. 2024) is not yet implemented, so the six showcase forbidden lines—where E2/M1 admixtures can enter the emissivity—are not covered by the cited rank-1 formalism in a justified way. A second, independent problem is that Eq. 1 labels Zeeman components by M along the solar-radius quantization axis, but Zeeman eigenstates are M-diagonal only when the quantization axis is aligned with the magnetic field. The manuscript does not describe a rotation to the B frame before evaluating Eq. 1, so at the arbitrary B orientations found in MURaM/PSI models the transverse-Zeeman Q,U signatures in Figures 13–14 could be misassigned. The §2.2 pyCELP benchmark, which would partially anchor correctness, is mentioned but never shown quantitatively, and §4 reports disagreement exactly in the new Zeeman Q,U terms.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents P-CORONA, an open-source Fortran 90 code for synthesizing Stokes I, Q, U, and V of coronal emission lines in 3D MHD models. The code solves multilevel statistical equilibrium equations for the atomic density matrix following Landi Degl'Innocenti & Landolfi (2004), including anisotropic radiative pumping, Hanle and Zeeman effects, electron and proton collisions, and Doppler dimming/brightening. The authors demonstrate the code on six forbidden lines (Fe xiii, Fe xiv, Fe xi, Si ix, Si x) in MURaM and PSI coronal models, and they highlight new predictions for the Zeeman-induced linear polarization of Fe xiii 10747 Å and Si ix 39343 Å. The code is publicly available on GitLab and Zenodo, with documentation.","tokens_in":21419,"tokens_out":9985,"duration_ms":117284,"significance":"If the implementation is correct, P-CORONA fills a real community need: a single public tool that treats both forbidden and permitted lines with a consistent Hanle-Zeeman framework in 3D coronal models, going beyond the saturated-Hanle approximation used by CLE/pyCELP for forbidden lines. The paper has notable strengths: the theoretical basis is established LL04 formalism, no ad-hoc fitted parameters enter the results, the inputs are standard (CHIANTI atomic data, MURaM/PSI models), and the authors explicitly quantify the detectability limitations of the predicted transverse-Zeeman signals. The code availability and strong-scaling tests are also positive. The main concerns are that the claimed benchmark against pyCELP is not shown quantitatively, and the domain of validity of the key emissivity expression is under-specified for non-dipole transitions and for arbitrary magnetic-field orientations.","major_comments":[{"comment":"The paper states that LL04 theory is valid 'under the dipole approximation,' yet it applies Eq. (1) to all radiative transitions satisfying ΔJ=0,±1, including 'non-dipole transitions.' For an electric quadrupole or higher multipole transition, the emitted photon carries angular momentum beyond rank 1, and the (1 1 K) angular factors and 3-j symbols in Eq. (1) are not the correct emissivity. Since the six showcase forbidden lines are not explicitly classified as pure M1 or E1 transitions, and since the Casini et al. (2024) multipole extension is stated not to be implemented, the code's computed Q, U, and V for any modeled line with non-negligible E2/M1 or higher multipole admixture would be systematically incorrect. The authors should either restrict the stated applicability to pure electric/magnetic dipole (rank-1) transitions, explicitly demonstrate that each modeled line is dipole-dominated, or implement the multipole treatment they cite as future work.","section":"Section 2.2 (Eq. 1)"},{"comment":"The text specifies that the solar radius vector is the quantization axis for total angular momentum, but Zeeman sublevel energies are M-diagonal only when the quantization axis is aligned with the magnetic field. The manuscript does not describe a rotation to the magnetic-field reference frame before evaluating Eq. (1), nor does it explain how the LL04 magnetic-kernel treatment accounts for the arbitrary field orientation. In the MURaM/PSI models the magnetic field direction changes significantly, so the transverse-Zeeman Q and U signatures shown in Figures 13–14 could be misassigned if the profile functions in Eq. (1) are evaluated in the vertical-axis frame without including the required level mixing or frame rotation. The authors need to specify the exact reference-frame procedure, or cite the precise LL04 equations that make the vertical-axis expression exact for arbitrary B.","section":"Section 2.2 (quantization axis)"},{"comment":"The claimed validation against pyCELP is presented only qualitatively: the text states 'good agreement between the results obtained using both codes, except when accounting for the impact of the Zeeman effect on the Stokes Q and U profiles for some IR lines.' This is a load-bearing point for a code paper, and the disagreement is exactly in the new Zeeman Q,U terms. A quantitative comparison, such as maximum relative differences in I, Q, U, and V as functions of wavelength and line-of-sight position for Fe xiii and Si ix, should be included so that the central correctness claim can be verified.","section":"Section 2.2 (benchmark)"}],"minor_comments":[{"comment":"Several figure panels contain corrupted or unreadable text (e.g., 'R/uni2299' in Figure 12 and repeated mangled labels in the colorbar axis text of Figures 5–11); the final version should use clean vector labels.","section":"Figures 5–12"},{"comment":"The quantity called 'total linear polarization' should be defined precisely; because Q and U vary across the line profile, the frequency-integrated P = sqrt(Q^2+U^2) is not a simple integrated linear-polarization amplitude unless an explicit convention is given.","section":"Table 2"},{"comment":"For the six forbidden lines the manuscript should state the dominant radiative multipole character (M1, E2, or mixed) for each line, since this is directly relevant to the validity discussion in Section 2.2.","section":"Section 3.1"},{"comment":"The fixed-LOS geometry is clearly acknowledged as unsuitable for near-Sun observers; the text should also state whether the code can be rerun with a rotated LOS to approximate Solar Orbiter/Parker Solar Probe viewing, since several readers may attempt such an use.","section":"Section 2.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a useful code paper built on established theory, and the concerns raised are fixable in a revision: provide the pyCELP benchmark quantitatively, clarify the multipole and quantization-axis scope of Eq. (1), and state the multipolarity of the modeled lines. I would not reject on the basis of the physics assumptions themselves, because they are explicit and testable. If the authors can supply the benchmark and the frame-transformation details, the paper would be suitable for acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague, quick read of P-CORONA. It is exactly what it claims: a new open-source tool for synthesizing Stokes I,Q,U,V of coronal lines in 3D MHD models with Hanle and Zeeman effects, going beyond the saturated-Hanle-only CLE/pyCELP. The public release on Zenodo (with input files), the GUI, and the parallelization are real infrastructure. The application to six forbidden lines in four MURaM and two PSI models is a solid demonstration, and the negative result that transverse-Zeeman Q,U are practically undetectable with current instruments is useful and honestly stated. Do not skip the credit: this fills a genuine gap for DKIST/Aditya-L1 forward modeling. Soft spots: the pyCELP benchmark is mentioned but never quantified, and given that the reported discrepancy is exactly the new Zeeman Q,U behavior, is part of the paper's central claim. The abstract says on-disk and off-limb capability, but the paper only shows plane-of-sky syntheses; either show an on-disk case or soften the wording. The dipole-approximation limitation is stated (strong coupling for non-dipole transitions; Casini et al. 2024 not implemented), but the authors do not estimate the impact for the M1/E2 lines they showcase. On the stress-test quantization-axis question: the paper does not explicitly describe the rotation to the magnetic-field frame for the Zeeman profiles; in LL04 the magnetic kernel in the SEE handles field orientation, so I would not call it an error, but a one-paragraph clarification would prevent a justified reader complaint. Finally, the detectability statements deserve an instrument sensitivity baseline, or at least a caveat that they are signal-level only. None of this is fatal. The theory is standard, the code is public, the models are external, and there is no fitting. I would send this to a serious referee and expect acceptance after moderate revision. Bring it to reading group if you care about coronal magnetometry forward modeling. My verdict: conditional accept, mainly because the benchmark needs to be visible.","headline":"A useful, honest code paper that deserves refereeing; the missing quantitative benchmark and a few scope overclaims are the real soft spots.","tokens_in":22134,"tokens_out":4667,"would_cite":true,"duration_ms":55085,"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":"P-CORONA is a new open-source code that calculates the intensity and polarization of forbidden and permitted coronal lines in three-dimensional solar models, treating the Hanle and Zeeman effects consistently.","keywords":["Solar corona","Solar magnetic fields","Solar coronal lines","Stokes parameters","Hanle effect","Zeeman effect","Forward modeling","Density matrix polarization"],"falsifier":"Recompute the Stokes Q and U profiles of Si ix 39343 Å in the same active-region MURaM snapshot with a code that handles all radiative transition types in full; if the transverse-Zeeman amplitudes differ materially from P-CORONA's, the strong-coupling approximation fails for this line.","tokens_in":20982,"feed_emoji":"☀️","tokens_out":13813,"duration_ms":127331,"temperature":0.7,"pith_summary":"P-CORONA is a new open-source code that computes the intensity and polarization of coronal emission lines in three-dimensional models of the solar corona. It solves the statistical equilibrium equations for the multipolar components of the atomic density matrix and uses the resulting emissivities to synthesize all four Stokes parameters, including the Hanle and Zeeman effects and Doppler dimming and brightening. The paper argues that by using the general emissivity expression with magnetic splitting, the code accounts for the Zeeman effect on Stokes Q and U, not only on Stokes V, as earlier forbidden-line codes did. If correct, the community gains a single tool for forward modeling both forbidden and permitted lines in the 3D corona, with synthetic maps that can be compared with DKIST and other polarimetric observations. Applied to MURaM and Predictive Science models, the code predicts that transverse-Zeeman linear polarization in the tested forbidden lines is present but very weak, generally below current detectability.","feed_headline":"P-CORONA computes full Stokes polarization in 3D coronal models","feed_subtitle":"New code predicts what DKIST should see from 3D solar corona models","key_machinery":"The engine is the density-matrix theory of spectral line polarization: the atomic state is represented by multipolar components $\\rho^K_Q(\\alpha J)$, and the statistical equilibrium equations are solved for them, with radiative rates built from radiation-field tensors that encode the anisotropy of the incident radiation, including Doppler dimming and brightening. The emissivity formula is the load-bearing identity: it couples the density matrix to the emergent Stokes parameters through 3-$j$ symbols and the magnetically split line profile, so the Zeeman effect enters Q, U, and V, not just V. Non-dipole radiative transitions are handled under the strong coupling approximation, meaning they only affect level populations, and a Hanle-saturated mode is offered for faster computations when appropriate.","core_discovery":"The paper's central claim is that P-CORONA correctly synthesizes Stokes I, Q, U, and V for coronal lines in 3D models by implementing the density-matrix theory of spectral line polarization in full, rather than restricting to the saturated Hanle regime or to the weak-field approximation for Stokes V. The code solves the statistical equilibrium equations for the multipolar components $\\rho^K_Q$ of the atomic density matrix for a multilevel atom, then computes frequency-dependent emissivities via the general expression $\\epsilon_i(\\nu,\\Omega)$ that retains the magnetic sublevel splitting in the Voigt and Faraday-Voigt profiles. This is what allows the Zeeman effect to enter all Stokes parameters, including the transverse-field contribution to Q and U; the paper demonstrates the effect on Fe xiii 10747 Å and Si ix 39343 Å in an active-region MURaM model, finding that the transverse Zeeman signatures are visible in the synthesized profiles but too weak for current instruments. The same implementation is applied to six forbidden lines in four MURaM snapshots and to Fe xiii 10747 Å in two large-scale Predictive Science models, and the code is released openly.","pith_inferences":["A natural extension of the code would be to apply it to EUV permitted lines such as O vi 1032 Å and Ne viii 770 Å, where the Hanle effect is not saturated; if the implementation handles them, P-CORONA could produce maps of field strength rather than only field orientation.","The strong coupling approximation is the main place where the computed Q, U, and V could be underestimated; implementing a full multipole treatment for a line like Si ix 39343 Å would give a direct test of how much polarization the non-dipole transitions actually carry.","The paper's conclusion that transverse Zeeman linear polarization is unobservable in these forbidden lines suggests that observers should prioritize longitudinal-field diagnostics via Stokes V or seek permitted lines with larger Zeeman sensitivity; P-CORONA could be used to survey candidate lines before committing telescope time."],"forward_implications":["P-CORONA provides synthetic Stokes I, Q, U, and V maps for forbidden and permitted coronal lines in 3D MHD models, enabling direct comparison with DKIST, UCoMP, and future coronagraphic polarimetry.","The full-emissivity treatment means that the Zeeman effect contributes to linear polarization (Q and U) as well as circular polarization (V), so the code can be used to test whether transverse-field signatures are observable in a given line and model.","The paper's calculations imply that in the lower corona, scattering linear polarization of the six forbidden lines is weak, with Fe xiv 5303 Å and Fe xiii 10747 Å the strongest; such predictions can guide line selection and integration times for observations.","For Si ix 39343 Å, transverse Zeeman signatures appear at low heights in strong-field active regions but remain near the detectability limit, giving a concrete upper bound on what current instruments can expect.","Because the code handles spherical Predictive Science models extending to several solar radii, the same synthesis can be applied to off-limb and large-scale coronal observations, not only disk-center geometries."],"supporting_citations":[{"why":"Supplies the density-matrix theory, statistical equilibrium equations, and emissivity formula that P-CORONA implements.","marker":"LL04"},{"why":"Provides the pyCELP code and benchmark results against which P-CORONA's forbidden-line Stokes profiles are compared.","marker":"Schad & Dima 2020"},{"why":"Provides the four MURaM coronal snapshots used for the illustrative synthesis.","marker":"Rempel 2017"},{"why":"Demonstrated P-CORONA's permitted-line capability and supplies the Doppler-dimming expression and PSI model description.","marker":"Supriya et al. 2021"},{"why":"Supplies CHIANTI version 10 atomic data for radiative and collisional rates.","marker":"Del Zanna et al. 2021"},{"why":"Recent DKIST observations of Fe xiii 10747 Å define the observational context for the detectability discussion.","marker":"Schad et al. 2024"},{"why":"Introduces the strong coupling approximation used for non-dipole radiative transitions.","marker":"Judge et al. 2006"},{"why":"Presents the multipole extension of the theory that the authors plan to implement in the future.","marker":"Casini et al. 2024"}],"fun_headline_variants":["P-CORONA synthesizes full Stokes IQUV in 3D coronal models","New code P-CORONA computes coronal line polarization in 3D models","P-CORONA: open-source synthesis of coronal Stokes profiles","Hanle and Zeeman effects in 3D coronal synthesis with P-CORONA","P-CORONA predicts DKIST-visible polarization from 3D coronal models"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the strong-coupling approximation: any radiative transition not satisfying the dipole selection rules is assumed to change only level populations, never polarization; if such transitions do polarize the upper levels of the forbidden lines, the computed Stokes Q, U, and V maps would be systematically too small.","fun_headline_variants_meta":{"raw":{"variants":["P-CORONA synthesizes full Stokes IQUV in 3D coronal models","New code P-CORONA computes coronal line polarization in 3D models","P-CORONA: open-source synthesis of coronal Stokes profiles","Hanle and Zeeman effects in 3D coronal synthesis with P-CORONA","P-CORONA predicts DKIST-visible polarization from 3D coronal models"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000736,"raw_usage":{"total_tokens":3348,"prompt_tokens":1062,"completion_tokens":2286,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":678,"completion_tokens_details":{"reasoning_tokens":2173}},"tokens_in":678,"tokens_out":2286,"duration_ms":18705,"temperature":1.0,"reasoning_tokens":2173,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:51:46.791198+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the Stokes Q and U profiles of Si ix 39343 Å in the same active-region MURaM snapshot with a code that handles all radiative transition types in full; if the transverse-Zeeman amplitudes differ materially from P-CORONA's, the strong-coupling approximation fails for this line.","supporting_citations":[{"cited_title":"D., Trujillo Bueno, J., de Vicente, ´A., & del Pino Alem´ an, T","cited_arxiv_id":null,"evidence_quote":"Demonstrated P-CORONA's permitted-line capability and supplies the Doppler-dimming expression and PSI model description."},{"cited_title":"G., Low, B","cited_arxiv_id":null,"evidence_quote":"Introduces the strong coupling approximation used for non-dipole radiative transitions."},{"cited_title":"A unifying polarization formalism for electric- and magnetic-multipole interactions","cited_arxiv_id":"2409.01197","evidence_quote":"Presents the multipole extension of the theory that the authors plan to implement in the future."}],"review_version":1}