{"id":"146d4565-7795-4cb9-83ea-8ae10237ca3b","arxiv_id":"2502.04169","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"New 3D simulations of the eccentric hot Jupiter TOI-150b show that a locally calculated magnetic drag weakens the equatorial jet and produces inter-orbit variability in high-resolution emission spectra.","lead":"This paper models the 3D atmosphere of the eccentric hot Jupiter TOI-150b with and without magnetic drag, and predicts that the combination of an eccentric orbit and a temperature-dependent magnetic force changes wind patterns and creates orbit-to-orbit variations in its spectrum. A general reader might care because these are concrete, testable predictions for JWST and ground-based high-resolution spectrographs.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Monotonic scaling with field strength is contradicted by the paper's own Figures 5 and 13: the 3G case strengthens the upper-atmosphere jet relative to 0G, and the 10G case shows less Doppler broadening than 30G.","rationale":"The reader's weakest_assumption focused on the unconstrained pseudo-synchronous rotation period, which is a legitimate limitation but is explicitly acknowledged in Section 4.3. The reader's rationale also flagged the non-monotonic magnetic scaling, identifying it as a reason for the conditional verdict. I regard the scaling inconsistency as the most load-bearing concern because it is an internal contradiction in the manuscript's central claim: the abstract and conclusions assert monotonic dependence on field strength, while Figures 5 and 13 show clear exceptions (3G stronger than 0G for the jet; 10G less broadening than 30G). This is not a question of uncertain external parameters but a direct mismatch between the claims and the presented model output. The concrete test is purely reanalysis of existing figures, so it can settle the issue immediately. The rotation-period concern, while real, would require new simulations and does not contradict any internal statement of results. Therefore, the reader's CONDITIONAL verdict remains appropriate: the paper should revise or qualify the scaling claims before the central conclusions are accepted.","tokens_in":21077,"tokens_out":8914,"duration_ms":88875,"concrete_test":"Quantitatively reproduce Figure 13 by computing the full-width-at-80%-maximum of the cross-correlation between Doppler-on and Doppler-off spectra for the 10G and 30G models at all phases; if 10G has smaller width than 30G, the stated monotonic decrease with field strength is false. Similarly, recompute the zonal wind profile comparison in Figure 5 at pressures below 0.01 bar for the 0G and 3G cases; if the 3G jet is stronger, the conclusion that increasing field strength weakens the jet must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim in the abstract and conclusions is that the strength and magnitude of the eccentric-orbit circulation and spectroscopic effects 'scale with the chosen global magnetic field strength' and that Doppler broadening decreases 'as magnetic field strengths increase.' This monotonic scaling is not supported by the paper's own results. Section 3.2 states that 'the jet in the 3G case is stronger in the upper atmosphere than the 0G case' (Figure 5), directly contradicting the conclusion that increasing field strength weakens the equatorial jet. More importantly, Section 3.5.1 and Figure 13 show that the 10G model has the least Doppler broadening at every phase, even less than the 30G model except at phase=0. Thus the claimed monotonic decrease in broadening with field strength fails at the 10G vs 30G comparison. These inconsistencies are internal to the manuscript, not a matter of external consensus, and they undercut the headline scaling statement that is a core part of the paper's contribution. The pseudo-synchronous rotation assumption is an acknowledged limitation, but the scaling inconsistency is a verifiable contradiction within the presented data.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"Beltz et al. present the first 3D general circulation model study of an eccentric hot Jupiter (TOI-150b, e=0.26) with kinematic MHD drag. They run eccentric and circular models with dipole field strengths B = 0, 3, 10, and 30 G, using a locally calculated, temperature-dependent drag timescale from Perna et al. (2010). They find that eccentricity alone narrows and strengthens the equatorial jet, produces high-latitude westward winds, and creates a phase-dependent thermal inversion, and that MHD drag weakens and narrows the jet, alters the dayside circulation, and reduces Doppler broadening and net Doppler shifts in simulated high-resolution emission spectra. The paper further explores 1D cloud models and low-resolution spectra, and concludes that the eccentric-orbit effects and their spectroscopic signatures scale with the chosen magnetic field strength.","tokens_in":21301,"tokens_out":6011,"duration_ms":62164,"significance":"The study is a useful extension of the kinematic MHD approach to a physically important regime: eccentric hot Jupiters with periodically varying irradiation and thermal ionization. The drag prescription is locally calculated and not fitted to the target results; the magnetic field strength is varied as a free parameter, and the reported jet weakening, circulation changes, and spectral signatures emerge from the simulations. The work also connects 3D dynamics to concrete, testable predictions, such as phase-dependent thermal inversions, the 4.2 micron CO2 feature, and the 6.2 micron H2O emission/absorption switch. If the claims are appropriately qualified, the paper will be a valuable reference for observers targeting eccentric hot Jupiters and for modelers incorporating magnetic effects beyond simple uniform Rayleigh drag.","major_comments":[{"comment":"The abstract states that the strength and magnitude of the eccentric-orbit circulation effects 'scale with the chosen global magnetic field strength,' and Section 5 states that 'Increasing the magnetic field strength weakens the equatorial jet.' This monotonic scaling is contradicted by the paper's own results: Section 3.2 states that 'the jet in the 3G case is stronger in the upper atmosphere than the 0G case' (Figure 5), and the Figure 5 caption claims a weaker and narrower jet with increasing field strength. Because the scaling claim is part of the headline contribution, the abstract, conclusions, and figure caption should be revised to describe a non-monotonic or phase-dependent dependence, or to explicitly separate the 3G exception from the strong-field (10G/30G) trend.","section":"Abstract, Section 3.2, Figure 5, Section 5"},{"comment":"The claim of 'decreased Doppler broadening as magnetic field strengths increase' is not supported by the presented data. Section 3.5.1 states that 'the 10G model has the least amount of broadening at each phase, showing even less broadening than the 30G model except for at phase=0,' and Figure 13 plots this directly. The broadening is therefore not monotonically decreasing with field strength. The abstract and conclusions should be reworded to state that active magnetic drag reduces Doppler broadening relative to the drag-free case, with a non-monotonic dependence across the 3G, 10G, and 30G models, and the authors should either explain the 10G versus 30G behavior or present it as an apparently nonlinear response.","section":"Abstract, Section 3.5.1, Figure 13, Section 5"},{"comment":"The acknowledged limitation about pseudo-synchronous rotation should be given more weight in the presentation of the Doppler-shift predictions. As written, Section 4.3 notes that the 4.14-day rotation period is observationally unconstrained and that different rotation rates would alter the circulation and Doppler shifts, but the abstract presents the spectroscopic variability and Doppler-broadening trends as robust eccentric-orbit signatures. Either the abstract should carry a brief caveat, or the discussion should include a sensitivity estimate, such as a comparison against a synchronously rotating eccentric model or an order-of-magnitude estimate of the expected change in net Doppler shift.","section":"Section 4.3"}],"minor_comments":[{"comment":"The footnote contains a typo: 'timsecale' should be 'timescale.'","section":"Section 3.3, footnote 4"},{"comment":"The opening sentence, 'Our work is builds upon Kataria et al. (2013),' should be corrected to 'Our work builds upon Kataria et al. (2013).'","section":"Section 4.2"},{"comment":"The phrase 'the eccentric model has a faster rotation rate then the circular model' should read 'than the circular model.'","section":"Section 3.1, Figure 3 caption"},{"comment":"The quantity 'Full-Width 80%-Max' is not defined in the text; a sentence explaining the normalization and the cross-correlation width measure would improve reproducibility.","section":"Section 3.5.1, Figure 13"},{"comment":"The discussion attributes the circulation changes to 'the inclusion of eccentricity,' but the eccentric model also uses a faster pseudo-synchronous rotation period. The text acknowledges this on the previous page, but a one-sentence reminder in Section 3.1 would prevent readers from misinterpreting the comparison as isolating eccentricity alone.","section":"Section 3.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of ApJ and the underlying simulations appear internally consistent, with a physically motivated drag prescription and no obvious circularity. The main obstacle to acceptance is the overstated monotonic-scaling language in the abstract and conclusions, which is contradicted by the paper's own Figures 5 and 13. This is fixable with careful rewording and additional discussion; I do not see a need to rerun the simulations. The pseudo-synchronous rotation caveat is real but already acknowledged; it should be elevated slightly rather than treated as a reason for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid, clearly written modeling study that does something new — the first kinematic MHD GCM of an eccentric hot Jupiter — and it produces concrete, testable predictions for TOI-150b. The eccentricity-only effects (narrowed equatorial jet, westward mid-latitude winds, phase-dependent thermal inversion) reproduce Kataria et al. 2013 and Lewis et al. 2014, which is a good check. The new contribution is the locally calculated, temperature-dependent magnetic drag applied in an eccentric orbit, plus the resulting spectral predictions: inter-orbit Doppler variability that peaks around 10G, and reduced line broadening at higher field. Those are worth taking seriously.\n\nThe main soft spot is real and internal. The abstract and the conclusions say the effects “scale with the chosen global magnetic field strength” and that Doppler broadening decreases as field strength increases. But Section 3.2 reports that the 3G jet is stronger than the 0G jet in the upper atmosphere (Fig 5), and Section 3.5.1 reports that the 10G model has less broadening than the 30G model at nearly all phases (Fig 13). The monotonic framing is not supported by their own output. They do flag the 3G exception in the text, but the abstract and the conclusion bullets still state the monotonic version. That needs to be fixed, by dropping “scale with” or by reframing as “generally weaker, with exceptions at intermediate field strengths.”\n\nThe pseudo-synchronous rotation assumption is a genuine limitation, and they acknowledge it (Section 4.3, unconstrained observationally). That means the specific Doppler-shift predictions are less robust than the thermal structure results. A rotation-period sensitivity test would strengthen the paper; its absence is a caveat, not a fatal flaw.\n\nMethods come from established, published codes (RM-GCM, picket-fence radiative transfer, ExoMol opacities, EGP+ for clouds), and the parameter choices are transparent. The drag prescription is physical rather than fitted to the headline outputs, so the circularity burden is low. Citation pattern is fair; self-citations are to their own published prescription.\n\nRecommendation: send to peer review. The scaling claim needs revision, and the rotation sensitivity would help, but the core modeling and the spectroscopic predictions deserve referee time.","headline":"Solid first kinematic-MHD GCM of an eccentric hot Jupiter with testable spectral predictions, but the abstract's monotonic scaling claim does not survive their own Figures 5 and 13.","tokens_in":21856,"tokens_out":2785,"would_cite":true,"duration_ms":24448,"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 first 3D general circulation models of an eccentric hot Jupiter with temperature-dependent magnetic drag predict that the planet's equatorial jet weakens and narrows as field strength grows, that a dayside thermal inversion appears…","keywords":["hot Jupiters","exoplanet atmospheres","general circulation models","kinematic magnetohydrodynamics","magnetic drag","eccentric orbits","thermal inversions","high-resolution emission spectroscopy"],"falsifier":"A JWST emission spectrum of TOI-150b near periastron and another near apastron would settle the phase-dependent-inversion claim: the models predict the roughly 6.2-micron water feature in emission at periastron (thermal inversion present) and in absorption at apastron (no inversion), with the periastron CO2 feature at roughly 4.2 microns weakening as field strength rises to 10 Gauss. Observing no such phase flip in the water feature would contradict the central mechanism.","tokens_in":20905,"feed_emoji":"🪐","tokens_out":7980,"duration_ms":68705,"temperature":0.7,"pith_summary":"This paper uses 3D general circulation models to ask what happens to the atmosphere of an eccentric hot Jupiter when magnetic drag is calculated from local conditions rather than assumed as a uniform frictional timescale. The target is TOI-150b, a 1.75-Jupiter-mass planet on an e=0.26 orbit whose equilibrium temperature swings between roughly 1300 K and 1700 K over each 5.86-day orbit. The authors find that the eccentric orbit alone narrows and strengthens the equatorial jet, adds westward mid- and high-latitude winds, and makes thermal inversions appear only near periastron; adding kinematic MHD drag weakens the jet and deepens the westward flow as field strength rises from 0 to 30 Gauss. Post-processed emission spectra then predict that stronger magnetic fields reduce Doppler broadening, and that a moderate 10-Gauss field produces the largest orbit-to-orbit variability in net Doppler shifts. A reader cares because these spectral signatures give observers a way to infer an exoplanet's magnetic field strength from ground-based high-resolution spectroscopy.","feed_headline":"Magnetic drag reshapes winds of an eccentric hot Jupiter","feed_subtitle":"New 3D models of TOI-150b tie field strength to jet width, thermal inversions, and orbit-to-orbit Doppler shifts.","key_machinery":"The central mechanism is a locally calculated magnetic drag timescale, applied in the momentum equation as a Rayleigh drag $-u/\\tau_{\\mathrm{mag}}$ with a matching Ohmic dissipation term in the energy equation. The timescale is $\\tau_{\\mathrm{mag}} = 4\\pi\\rho\\eta/(B^2|\\sin\\phi|)$, where $B$ is the assumed dipole field strength, $\\phi$ the latitude, $\\rho$ the density, and $\\eta = 230\\sqrt{T}/x_e$ the magnetic resistivity; the ionization fraction $x_e$ comes from the Saha equation summed over the first 28 elements. Because resistivity depends exponentially on temperature through ionization, the drag varies by orders of magnitude between the hot dayside and cooler nightside, so the same field strength slows winds strongly near periastron and barely at all near apastron. This temperature dependence is what couples orbital phase, magnetic field strength, and the predicted Doppler variability.","core_discovery":"The paper establishes that the first application of kinematic (temperature-dependent) magnetic drag to an eccentric hot Jupiter changes both the circulation and the observable spectra in a field-strength-dependent way. Compared with a circular-orbit simulation of the same planet, the eccentric model develops a narrower, stronger equatorial superrotating jet and westward flow at high latitudes that spreads toward mid-latitudes as the dipole field strength is increased to 10-30 Gauss. The thermal structure becomes phase-dependent: a dayside temperature inversion exists only near periastron, where irradiation peaks, and the day-night temperature contrast grows with field strength. When the 3D structures are post-processed into R=100,000 emission spectra, higher field strengths produce less Doppler broadening, and the 10-Gauss model shows the largest inter-orbit variation in net Doppler shift when the same hemisphere is viewed at different orbital phases, a signature the authors attribute to the local drag responding to the changing temperature structure.","pith_inferences":["Extending the paper's logic, hotter eccentric planets whose equilibrium temperature stays above the thermal ionization threshold all orbit would keep magnetic drag active year-round, so their inter-orbit variability could be even larger than what is predicted here for TOI-150b.","The prediction that 10 Gauss maximizes variability while 30 Gauss suppresses it implies a non-monotonic mapping between field strength and observable variability; if that mapping holds, single-epoch Doppler measurements cannot rank field strengths without multi-epoch, phase-resolved data.","Because the drag timescale depends on local temperature, any cloud or chemistry process that alters the temperature profile would shift where and when the drag acts; coupling the 3D model with the cloud prescription used here in 1D would test whether clouds mute the predicted Doppler variability.","The paper's estimate that the orbit has only a 22 percent geometric chance of showing a secondary eclipse means that a nondetection would not falsify the models, while a detection would constrain the orbital geometry and sharpen the phase-curve predictions."],"forward_implications":["If the models are right, TOI-150b's equatorial jet is weaker and narrower for stronger assumed fields, and westward flow extends from high latitudes toward mid-latitudes at 10-30 Gauss.","Thermal inversions on this planet should appear only near periastron and disappear near apastron, so emission spectra should flip from showing water in emission to water in absorption across the orbit.","Stronger magnetic fields should yield measurably less Doppler broadening in high-resolution emission spectra, and the 10-Gauss model predicts the largest orbit-to-orbit Doppler variability for a fixed viewing geometry.","Phase curves should peak before periastron, with the peak occurring closer to periastron as field strength increases.","Clouds, modeled here in 1D, should be more abundant near apastron and may vanish at periastron, with only small (no more than 5 percent) spectral effects."],"supporting_citations":[{"why":"Supplies the 4.14-day pseudo-synchronous rotation period assumed for the eccentric models, on which the circulation and Doppler predictions depend.","marker":"Hut (1981)"},{"why":"Derives the local magnetic drag timescale and the resistivity-ionization prescription that the kinematic MHD scheme implements.","marker":"Perna et al. (2010a)"},{"why":"First implementation of the locally calculated magnetic drag timescale in a hot Jupiter GCM; the method this paper applies to an eccentric planet.","marker":"Rauscher & Menou (2013)"},{"why":"Previous 3D GCM study of eccentric hot Jupiters that found the narrowed jet and westward mid-latitude winds; the eccentricity baseline this work extends with magnetic drag.","marker":"Kataria et al. (2013)"},{"why":"Applied kinematic MHD to an ultra-hot Jupiter and produced high-resolution emission spectra; supplies the comparison for spectra and the finding that drag lowers continuum flux.","marker":"Beltz et al. (2022b)"},{"why":"1D models of eccentric hot Jupiters that predicted the post-periastron peak-flux offset and the time-dependent thermal inversion that this paper's 3D models reproduce.","marker":"Mayorga et al. (2021)"},{"why":"Provides the eccentric orbital fit, semimajor axis, stellar age, and the secondary-eclipse upper limit used for model parameters and the eclipse probability calculation.","marker":"Kossakowski et al. (2019)"},{"why":"Provides the mass and radius used in the simulations.","marker":"Cañas et al. (2019)"},{"why":"Ray-tracing radiative transfer scheme used to post-process the GCM structures into emission spectra.","marker":"Zhang et al. (2017)"}],"fun_headline_variants":["Kinematic MHD narrows jet on eccentric hot Jupiter TOI-150b","Magnetic drag controls phase-dependent thermal inversions","First kinematic MHD study of eccentric hot Jupiter shows wind shifts","Field strength scales Doppler shifts in hot Jupiter emission"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes a pseudo-synchronous rotation period of 4.14 days for the planet, and the authors state plainly that this rotation rate is unconstrained observationally; all the predicted circulation changes, thermal-inversion timing, and Doppler shifts are computed against that assumed spin, so a different real rotation rate would change the wind pattern and the spectroscopic signatures.","fun_headline_variants_meta":{"raw":{"variants":["Kinematic MHD narrows jet on eccentric hot Jupiter TOI-150b","Magnetic drag controls phase-dependent thermal inversions","First kinematic MHD study of eccentric hot Jupiter shows wind shifts","Field strength scales Doppler shifts in hot Jupiter emission"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000644,"raw_usage":{"total_tokens":3003,"prompt_tokens":1030,"completion_tokens":1973,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":646,"completion_tokens_details":{"reasoning_tokens":1904}},"tokens_in":646,"tokens_out":1973,"duration_ms":13536,"temperature":1.0,"reasoning_tokens":1904,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T23:16:41.166330+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A JWST emission spectrum of TOI-150b near periastron and another near apastron would settle the phase-dependent-inversion claim: the models predict the roughly 6.2-micron water feature in emission at periastron (thermal inversion present) and in absorption at apastron (no inversion), with the periastron CO2 feature at roughly 4.2 microns weakening as field strength rises to 10 Gauss. Observing no such phase flip in the water feature would contradict the central mechanism.","supporting_citations":[],"review_version":1}