{"id":"89bb8cba-2686-4f8d-8529-2a64ac61f9ca","arxiv_id":"2602.01364","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"3D simulations of WASP-121b show a photoevaporative outflow bent by orbital motion into two spiral arms whose ~40 km/s velocities explain observed high-velocity Na and Hα absorption without super-rotation jets.","lead":"This paper simulates the ultra-hot Jupiter WASP-121b in three dimensions and argues that the fast sodium and hydrogen-alpha absorption signals seen during transits are produced by two spiral arms of escaping gas rather than by fast equatorial jet streams. It matters because multi-species transit spectra could become a practical probe of a star's ultraviolet output and wind environment around evaporating giant planets.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Hα spiral-arm signature depends on a single unvalidated atomic-rate choice (Janev et al. vs Osterbrock & Ferland) that can suppress the signal by ~40×","rationale":"The paper's kinematic case for spiral arms (Eq. 8) is parameter-free and compelling, and the Na velocity pattern is plausibly robust to the ×0.5 scaling and FUV uncertainty because the recombination balance (Eq. 7) and geometry set the peak positions. The Hα signal, however, lives at the interface where H 2s is populated by recombination and depleted by collisional 2s→2p transfer; a factor-40 difference in that rate determines whether any Hα secondary feature emerges. Since the paper explicitly flags this sensitivity (§3.2.3) but provides no sensitivity run, the central claim's Hα component is not yet established. This is a correctness risk, not a circularity or consensus issue, and it is testable by recomputing spectra with the alternative rate set. The lack of a super-rotation control run is a secondary concern because the paper's positive claim (arms can explain the features) is independent of whether a jet might also; the rate sensitivity directly threatens the positive claim.","tokens_in":22208,"tokens_out":4787,"duration_ms":50477,"concrete_test":"Recompute the fiducial synthetic Hα transmission spectra (§2.3, Fig. 9) using the Osterbrock & Ferland (2006) collisional 2s–2p rate coefficients in place of Janev et al. (2003) in the thermochemical network, holding all other parameters fixed. If the secondary high-velocity peaks at |Δv|≳40 km/s drop by ~40× and fall below the observed excess absorption, the spiral-arm explanation for Hα is not robust; if the peaks persist or the two rate sets yield similar H 2s populations at T≈10^4 K, n_e≈10^8 cm^-3, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the spiral-arm morphology reproduces the observed high-velocity Na and Hα absorption depends critically on the survival of the H 2s population at the dense/attenuated arm interface. The paper's own §3.2.3 notes that adopting Osterbrock & Ferland (2006) collisional 2s↔2p rate coefficients instead of Janev et al. (2003) suppresses the synthetic Hα amplitude by a factor of ~40. Since the Hα spectra are already scaled by 0.8 to match the data (Fig. 9), a factor-40 reduction would leave the synthetic Hα secondary peaks at |Δv|≳40 km/s far too weak to explain the observed excess absorption. The paper does not justify the choice of Janev et al. rates in this regime, nor does it provide a sensitivity study; the Hα leg of the central claim therefore rests on an unverified atomic-rate assumption. This is not merely an amplitude calibration issue—it determines whether the Hα spiral-arm signal exists at all.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents 3D hydrodynamic simulations of the ultra-hot Jupiter WASP-121b, coupling non-equilibrium thermochemistry, ray-tracing radiative transfer, and hydrodynamics in the Kratos GPU framework. The fiducial model produces a transonic photoevaporative outflow sculpted by stellar gravity and Coriolis forces into two spiral arms. The authors argue that this morphology naturally reproduces the observed high-velocity Na and Hα absorption patterns (redshifted ~20 km/s at ingress, blueshifted ~35 km/s at egress, and secondary peaks at |Δv|≳40 km/s) without invoking strong super-rotation jets. A parameter study varies FUV, EUV, X-ray, optical flux, dust, and stellar wind properties. The central kinematic claim is supported by order-of-magnitude estimates (Eqs. 4, 5, 8). The paper also emphasizes that different species trace different atmospheric reservoirs: Fe the inner rotation-dominated layers, Na the dense spiral arms, Hα and He 10830 Å the extended ionized regions.","tokens_in":22390,"tokens_out":3655,"duration_ms":39680,"significance":"If the two-arm, Coriolis-sculpted outflow structure is real and sufficiently populated, it would provide a new and physically attractive explanation for the asymmetric morning–evening velocity signatures observed in transmission spectroscopy of ultra-hot Jupiters, replacing the need for extreme super-rotation jets. The paper's kinematic mechanism is backed by parameter-free order-of-magnitude estimates that match the simulated ~40 km/s arm velocities, and the parametric study (Table 2) is a useful step toward understanding how different stellar environments affect observable tracers. These are genuine strengths. However, the quantitative validation of the model is weakened by post-hoc amplitude scalings and, more critically, by the extreme sensitivity of the Hα signature to a single atomic-rate choice that the authors themselves identify as capable of suppressing the signal by ~40×. The paper is a valuable contribution if these gaps can be closed or the claims appropriately softened.","major_comments":[{"comment":"The claim that the spiral-arm morphology reproduces the observed high-velocity Hα secondary peaks is load-bearing for the abstract, but it rests on an unvalidated choice of collisional 2s↔2p redistribution rates. The paper itself states that adopting Osterbrock & Ferland (2006) rates instead of Janev et al. (2003) suppresses the Hα absorption amplitude by a factor of ~40. Since the synthetic Hα is already scaled by 0.8 to match the data (Fig. 9), a factor-40 reduction would make the predicted Hα spiral-arm signal far too weak to explain the observed excess absorption. No sensitivity test is provided, and the rate choice is not justified for the relevant plasma conditions (T∼10^4 K, n_e∼10^8 cm^-3). This is not a minor calibration detail; it determines whether the Hα leg of the central claim exists at all. The authors should either run a test case with the alternative rate set, or substan","section":"§3.2.3"},{"comment":"The quantitative spectral comparison is partly imposed rather than emergent. Synthetic Na absorption is multiplied by 0.5 (Fig. 8), Hα by 0.8 (Fig. 9), and the observed excess-absorption maps are shifted by +3 km/s 'for better fitting of the trend' (Fig. 6 caption). The +3 km/s shift directly affects the velocity comparison that is central to the paper's argument. While the scaling factors are disclosed, they are not derived from independent constraints (e.g., abundance measurements or a stellar UV characterization), and the paper later admits that the true FUV likely lies between Model 0 and FUV10 (§4.1). The velocity morphology itself may be robust, but the spectral comparisons should be framed as consistency checks, not validations, unless the scalings and the velocity offset are justified from first principles or independent data.","section":"§3.2.2, Figs. 8, 9, and 6 caption"},{"comment":"The parametric study is broad but some results are presented in a way that is hard to evaluate quantitatively. For example, Model FUV10 is said to produce Na secondary peaks that appear 'less prominent than the observed data in the fiducial model', and the paper suggests the real FUV lies between Models 0 and FUV10. This implies that the fiducial model's Na amplitude agreement (after ×0.5 scaling) may be partly an accident of choosing the lower FUV. A quantitative comparison of the scaled synthetic spectra and observed data for all models in Table 2, ideally with a simple χ² or residual metric, would help the reader judge which parameter combinations are truly preferred. Without such metrics, the parametric conclusions remain qualitative.","section":"§4.1, Table 2"}],"minor_comments":[{"comment":"The abstract says 'coupling hydrodynamics with consistent non-equilibrium thermochemistry, ray-tracing radiative transfer, and hydrodynamics' — 'hydrodynamics' appears twice. Please correct.","section":"Abstract"},{"comment":"'with theT eq = 2400 K equilibrium temperature' — missing space after 'the'.","section":"§2.1"},{"comment":"'dense spiral rams' should be 'dense spiral arms'.","section":"§3.1.2"},{"comment":"Typo: 'relatively large velociteis' should be 'velocities'.","section":"§3.2.2"},{"comment":"'bottowm row' should be 'bottom row'.","section":"Figure 4 caption"},{"comment":"'equitorial slices' should be 'equatorial slices'.","section":"§3.1"},{"comment":"The definition of the dimensionless equivalent width is garbled: it appears to have both an integral factor 1/Δλ and a division by F_c Δλ in the integrand. Please clarify the normalization.","section":"Eq. (10)"},{"comment":"The entry 'Arcangeli et al. 2018b' is identical to '2018a' (same journal, volume, page). Also, the entry 'Xu, S., Wang, L., Ho, L. C., Cen, R., & Xu, S. 2026' lists Xu, S. twice as an author; this is likely a typo.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses an important problem and the kinematic mechanism is appealing. However, the Hα leg of the central claim is critically sensitive to a single atomic-rate choice acknowledged in the text, and the amplitude/velocity comparisons involve disclosed but unphysical scaling and shifting. These are fixable in a revision, provided the authors either test the alternative rates, quantitatively assess the scaling factors, or appropriately limit the claims. I would not reject, but the current version overstates the level of validation. I recommend major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a serious simulation paper with a genuinely new interpretive claim. The central idea is that the two high-velocity features in WASP-121b's Na and Hα transmission spectra are the leading and trailing arms of a photoevaporative outflow, bent by Coriolis deflection, not a super-rotation jet. That is worth taking seriously.\n\nWhat is actually new: no one in the cited literature has put these observed 20–50 km/s features on the escape-flow geometry. The kinematic skeleton is the strongest part. Eq. 8 gives ~40 km/s from tidal acceleration and Coriolis deflection with no fitted constants, and the simulated LOS velocities land in the observed range. The species stratification (Fe inner, Na dense arms, Hα/He* interface) is physically plausible and gives a coherent picture of where the lines form. The parametric study is a real bonus: FUV, EUV, X-ray, and wind each change the morphology in ways that could be checked against future data. I believe the authors when they say the two-arm morphology emerges from the stated hydrodynamics; the velocity scale is parameter-free.\n\nWhere it gets softer. The amplitude agreement is not emergent. Na synthetic spectra are scaled by 0.5, Hα by 0.8, and the observed velocity zero-point is shifted by +3 km/s. The paper is transparent about this, but it means the 'reproduction' claim is about line shape and velocity, not absolute depth. The bigger soft spot is the Hα rate sensitivity: §3.2.3 admits that Osterbrock & Ferland (2006) rates would suppress the Hα amplitude by ~40×. Since the Hα comparison is one of the two observational pillars, this deserves a sensitivity run. It may not break the Na-based case, but as written the Hα leg rests on a single rate database. I also would have liked a control run with an imposed super-rotation jet to show that the data genuinely disfavor it; the dismissal is argumentative rather than demonstrated. The lack of public code/data is a reproducibility concern, though common for this kind of GPU framework paper.\n\nNet: the central kinematic claim—that escape-flow geometry can produce the high-velocity Na pattern—holds up. The Hα leg is conditional until the atomic physics is settled. This paper deserves a serious referee. I'd send it out and ask for: a sensitivity test on the 2s↔2p collision rates, a jet control run, release of at least the synthetic spectra and time-averaged fields, and a more honest presentation of the scaling factors. My verdict would be conditional accept, not reject.","headline":"A plausible new explanation for WASP-121b's high-velocity Na/Hα features—Coriolis-sculpted outflow arms instead of a super-rotation jet—with real kinematics but amplitude comparisons partly imposed.","tokens_in":23005,"tokens_out":2304,"would_cite":true,"duration_ms":23810,"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":"WASP-121b's escaping atmosphere forms two spiral arms that reproduce the planet's high-velocity sodium and hydrogen absorption without invoking super-rotating jet streams.","keywords":["ultra-hot Jupiter","WASP-121b","atmospheric escape","transmission spectroscopy","hydrodynamics","thermochemistry","sodium absorption","H-alpha"],"falsifier":"A high-cadence transmission spectrum of WASP-121b that resolves the sodium D lines and H-alpha during a single transit, with signal-to-noise high enough to measure the amplitudes of the secondary and tertiary velocity peaks without scaling factors. If the secondary peaks are absent or appear at different phases, or if the measured FUV flux implies a sodium photoionization rate that cannot be balanced by recombination in the arms, the spiral-arm explanation would fail.","tokens_in":21943,"feed_emoji":"🪐","tokens_out":2901,"duration_ms":33422,"temperature":0.7,"pith_summary":"This paper tries to establish that the ultra-hot Jupiter WASP-121b's outer atmosphere is not dominated by fast, uniform jet streams, but instead by a photoevaporative outflow sculpted into two spiral arms by stellar gravity and the Coriolis force. Those arms accelerate gas to roughly 40 km/s along the line of sight, naturally producing the observed asymmetric sodium and H-alpha velocity features: redshifted absorption at ingress, blueshifted at egress, with secondary peaks beyond 40 km/s. If correct, this reframes how we read transmission spectra of ultra-hot Jupiters: different atomic and ionic species trace distinct atmospheric reservoirs, and the same spiral structure can explain multiple tracers simultaneously. The paper supports this claim with coupled 3D hydrodynamics, non-equilibrium thermochemistry, and ray-traced synthetic spectra, and it explores how stellar UV flux, X-rays, and stellar winds would shift the observables.","feed_headline":"Spiral arms, not jet streams, drive WASP-121b's sodium signal","feed_subtitle":"A 3D simulation ties the planet's asymmetric ingress-egress absorption to Coriolis-shaped outflow arms, not to fast winds.","key_machinery":"The central mechanism is the two-armed spiral outflow: gas heated on the dayside overflows the Roche lobe near L1 and L2, then Coriolis deflection in the orbiting frame bends the streams into leading and trailing spiral arms that reach supersonic speeds. The argument is carried by a GPU-accelerated 3D simulation that couples hydrodynamics with a 32-species, ~185-reaction non-equilibrium chemical network and ray-tracing radiative transfer, which produces synthetic transmission spectra comparable to data. The key physical balances are: sodium neutrality maintained by recombination against FUV photoionization in dense arms; H 2s population via recombination and charge exchange at the dense-to-e","core_discovery":"The paper's central claim is that the observed high-velocity neutral sodium and H-alpha absorption features of WASP-121b arise from two dense spiral arms in the planet's escaping upper atmosphere. These arms form as gas spills over the L1 and L2 Lagrangian points and is deflected by the Coriolis force in the co-rotating frame; the leading arm points toward the star and is preferentially seen at ingress (redshifted ~20 km/s), while the trailing arm lags behind the orbit and is seen at egress (blueshifted ~35 km/s). Neutral sodium survives in the dense arms because recombination offsets FUV photoionization, while H-alpha (from the H 2s state) traces the interface between the dense arms and the","pith_inferences":["If the two-arm interpretation holds, it suggests that many ultra-hot Jupiters with similar orbital parameters could show analogous asymmetric absorption patterns, and that the amplitude of the asymmetry encodes the velocity field of escape rather than the strength of atmospheric circulation.","A natural testable extension would be to look for transit-to-transit variability in the secondary and tertiary sodium peaks, since the paper links such variability to Kelvin-Helmholtz instabilities in the shear layer between prograde deep winds and retrograde outflow.","The paper leaves open the possibility that FUV radiation itself is the main controller of where neutral sodium can survive; if so, measuring the true FUV spectrum of WASP-121 would sharpen all the predictions.","The model's reliance on scaled amplitudes suggests that a completely parameter-free match would require an independent determination of sodium abundance and stellar FUV flux, a likely next step for the authors or observers."],"forward_implications":["If the spiral-arm morphology is correct, the morning-evening asymmetry in sodium and H-alpha absorption is a direct diagnostic of outflow geometry, not of wind jets.","The same model predicts that H-alpha should show secondary absorption peaks at |Δv| ≳ 40 km/s near ingress and egress, a feature the paper finds in existing data.","The strong response of He 10830 A equivalent width to stellar wind compression means metastable helium can serve as a probe of wind-planet interactions.","Changes in FUV, EUV, and X-ray fluxes produce distinguishable changes in the sodium and helium signals, so multi-wavelength transit spectroscopy can constrain the high-energy stellar environment.","The paper implies that ground-based high-resolution observations of Na, H-alpha, and Fe during the same transit could test the layered interpretation directly."],"fun_headline_variants":["Coriolis arms, not jets, shape WASP-121b's sodium signal","3D model: Spiral arms drive WASP-121b's high-velocity sodium","WASP-121b's sodium comes from spiral arms, not jet streams","Spiral outflows reproduce WASP-121b's sodium and H-alpha","Simulation links WASP-121b's Na absorption to Coriolis arms"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The simulated sodium and H-alpha absorption amplitudes are rescaled by factors of 0.5 and 0.8 to match the observations, and the H-alpha signal depends on adopting specific collisional rate coefficients that, with some alternative rates, would make the predicted absorption roughly 40 times weaker than observed.","fun_headline_variants_meta":{"raw":{"variants":["Coriolis arms, not jets, shape WASP-121b's sodium signal","3D model: Spiral arms drive WASP-121b's high-velocity sodium","WASP-121b's sodium comes from spiral arms, not jet streams","Spiral outflows reproduce WASP-121b's sodium and H-alpha","Simulation links WASP-121b's Na absorption to Coriolis arms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000207,"raw_usage":{"total_tokens":1286,"prompt_tokens":843,"completion_tokens":443,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":587,"completion_tokens_details":{"reasoning_tokens":336}},"tokens_in":587,"tokens_out":443,"duration_ms":4713,"temperature":1.0,"reasoning_tokens":336,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T05:41:06.279370+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A high-cadence transmission spectrum of WASP-121b that resolves the sodium D lines and H-alpha during a single transit, with signal-to-noise high enough to measure the amplitudes of the secondary and tertiary velocity peaks without scaling factors. If the secondary peaks are absent or appear at different phases, or if the measured FUV flux implies a sodium photoionization rate that cannot be balanced by recombination in the arms, the spiral-arm explanation would fail.","supporting_citations":[],"review_version":1}