{"id":"29f132d4-688c-4c03-8897-c56878d31db3","arxiv_id":"2507.10110","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A collisional ion test-particle model driven by hybrid-simulation fields underestimates Rosetta plasma densities at 67P by 5-10 times, pointing to the adiabatic electron assumption as the cause.","lead":"This paper adapts a 3D test-particle model to simulate cometary ions at Rosetta's comet 67P, then compares the modeled plasma density with spacecraft measurements. The model underestimates the measured density by 5 to 10 times, which the authors attribute to overly strong electric fields coming from the hybrid simulation's simplified treatment of electrons.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The adiabatic-electron attribution for the 5-10x underestimate is not tested: one hybrid snapshot with fixed upstream conditions, plane-averaged and collisionless fields, and a uniform electron-impact boost can explain the gap without invoking electron closure; a controlled closure run is needed.","rationale":"The model is carefully built and validated in the collisionless limit against AMITIS (Appendix A2) and in the collisional limit against the 1D model of Lewis et al. (2024) (Appendix A1); the hybrid fields are publicly archived (Moeslinger & Gunell 2024). Those checks support the numerical implementation. What they do not support is the paper's central physical conclusion. The 5-10x underestimate in Figure 10a is a single-model-versus-aggregate-data comparison. To claim that the underestimate is caused specifically by the adiabatic electron closure, one must show that changing that closure, while leaving the rest of the model fixed, removes or substantially reduces the gap. No such experiment is reported. The paper's own Sections 5.3.1 and 5.3.2 list alternative and co-varying limitations: the uniform electron-impact boost, the collisionless field generation, and the fixed upstream state. These are not merely caveats; they are alternative explanations for the same observed mismatch. The reader's CONDITIONAL verdict is therefore appropriate, and I do not see grounds to strengthen or weaken it. The one check that would convert the conditional into a claim is a controlled rerun with a non-adiabatic or collisional electron closure. Without it, the paper should be read as a well-posed hypothesis, not an established result.","tokens_in":17330,"tokens_out":10219,"duration_ms":120211,"concrete_test":"Re-run the Section 4.2 comparison with the AMITIS electron pressure replaced by a non-adiabatic closure, e.g. the electron pressure equation of Koenders et al. (2015) or the collisional electron test-particle model of Stephenson et al. (2023), keeping all other inputs in Table 1 identical. If the terminator-plane density at 20-100 km does not rise by the 5-10x needed to reach the MIP/LAP median, the adiabatic electron closure is not the cause of the underestimate. If the discrepancy vanishes, the attribution is supported; if it persists, the paper's central conclusion requires a different explanation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4.2 compares a single AMITIS run (Table 1) to MIP/LAP medians binned over 19.85 < log10(Qnu) < 20.15. The model output is the entire terminator plane averaged into 20 km radial bins, not the actual Rosetta ephemeris; the data median merges many intervals with different heliocentric distances, upstream solar wind states, and spacecraft potentials. The quoted 5-10x underestimate is therefore not yet a robust estimate of what the hybrid model would predict for the specific plasma Rosetta sampled. More importantly, even if the discrepancy is taken at face value, the causal claim in the Conclusions - that Eqs. 12-13, the adiabatic electron closure, drives the over-fast transport - is not tested anywhere in the paper. No simulation varies gamma or the electron pressure treatment; no run shows that replacing Eqs. 12-13 with a collisional/cooled electron closure raises the densities by the required factor. The hybrid fields themselves are generated without self-consistent ion-neutral collisions; ion-neutral momentum transfer is added only in the test-particle post-processing, so the high ambipolar field reflects the absence of self-consistent ion-neutral friction as well as the electron closure. Section 5.3.1 also implements electron-impact ionisation as a uniform boost of the photoionisation frequency, although Stephenson et al. (2023) found near-nucleus enhancement; a radial production-rate error acts in the same direction as the reported density shortfall. The conclusion that kinetic, collisional electron modelling is required is an interpretation of one model-data gap, not a demonstrated mechanism.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper adapts a 3D collisional test-particle model, originally developed for electrons, to cometary ions (H2O+, H3O+, NH4+) at 67P, using electric and magnetic fields from an AMITIS hybrid simulation representative of 2.5-3 au and an outgassing rate Q = 5.4e26 s^-1 (Table 1). The model includes energy-dependent ion-neutral collisions and is validated against the authors' 1D ion acceleration model (Lewis et al. 2024) and against the collisionless AMITIS output (Appendix A). In Section 4.2, the modelled terminator-plane densities, averaged in 20 km radial bins, are compared with MIP/LAP electron densities binned over 19.85 < log10(Q nu) < 20.15. The paper finds that the model densities are 5-10 times lower than the spacecraft data, and attributes this discrepancy to the adiabatic electron closure in the hybrid fields (Eqs. 12-13), concluding that kinetic collisional modelling of electrons is necessary.","tokens_in":17708,"tokens_out":3476,"duration_ms":38177,"significance":"If the central result holds, the paper is valuable: it identifies a concrete limitation of hybrid simulations with adiabatic electron pressure for the inner coma of 67P at intermediate outgassing, and it quantifies the effect of ion-neutral collisions on density and bulk velocity. The study has clear strengths: energy-dependent ion-neutral cross sections, validation against both a 1D fluid model and a collisionless hybrid output, comparison with independent Rosetta MIP/LAP data, and the public availability of the AMITIS simulation data. However, the significance is tempered because the causal attribution to the adiabatic electron closure is not directly tested, and the model-data comparison rests on one hybrid snapshot with fixed upstream conditions.","major_comments":[{"comment":"The central 5-10x underestimate is not yet a robust estimate of what the hybrid model would predict for the plasma Rosetta actually sampled. The comparison uses one AMITIS snapshot with fixed upstream solar wind parameters (Table 1), and the model output is the full terminator plane averaged into 20 km radial bins, not the spacecraft ephemeris. The MIP/LAP median is taken over a broad bin 19.85 < log10(Q nu) < 20.15, which aggregates many intervals with different heliocentric distances, upstream conditions, and spacecraft potentials. Before attributing the gap to the electron closure, the paper should either degrade the model along the actual Rosetta trajectory or show that the discrepancy is robust to the choice of radial bins and to variations of the upstream parameters within the 2.5-3 au window.","section":"Section 4.2, Figure 10"},{"comment":"The causal claim that the adiabatic electron assumption (Eqs. 12-13) drives the over-fast ion transport and the resulting 5-10x density shortfall is not tested anywhere in the paper. No simulation varies the adiabatic index gamma or replaces Eqs. 12-13 with a collisional/cooled electron closure; no control run demonstrates that such a replacement raises the modelled densities by the required factor. Furthermore, the hybrid fields are generated without self-consistent ion-neutral collisions, so the high ambipolar field may reflect the absence of ion-neutral friction in the hybrid model as well as the electron closure. The conclusion should be reworded as a hypothesis supported by earlier work, or a controlled closure run should be added.","section":"Conclusions and Section 5.3.2"},{"comment":"The treatment of electron-impact ionisation as a uniform boost of the photoionisation frequency is a potential confound for the reported underestimate. Stephenson et al. (2023) found that electron-impact ionisation is enhanced in the first ~100 km near the nucleus; if the true production rate has such a radial enhancement, the test-particle model's production profile is too low near the nucleus, which acts in the same direction as the reported density shortfall. The paper should quantify the sensitivity of the 5-10x ratio to a radially peaked production profile, or implement a production profile consistent with the kinetic electron simulations, before attributing the entire discrepancy to transport.","section":"Section 5.3.1, Table 1"}],"minor_comments":[{"comment":"The text states that 'A value of un = 700 km s^-1 is assumed'; this should be 0.7 km s^-1 to be consistent with Table 1 and the surrounding discussion.","section":"Section 4.1, Eq. (11)"},{"comment":"In the sentence describing the field-free chemistry-free model, the parenthetical reads 'ui = ui = 0.7 km s^-1'; this should be 'ui = un = 0.7 km s^-1'.","section":"Section 5.3"},{"comment":"The sentence 'The Sun in the +X direction' is missing a verb and should read 'The Sun is in the +X direction'; similarly, 'the interplanetary magnetic field is oriented along the +Y axis' would improve readability.","section":"Section 2.2, page 4"},{"comment":"The notation vmean_n,th is not clearly defined; the text refers to the mean speed in the frame of the neutral gas, but the symbol vmean_n,th should be defined explicitly to avoid confusion with the most probable speed v_n,th defined just above.","section":"Eq. (9)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of MNRAS and the modelling work is generally careful. The main gap is that the headline causal attribution to the adiabatic electron closure is not tested by any controlled run; a revision that either performs a non-adiabatic closure experiment or substantially softens the causal claim would address my main concern. I would also encourage the authors to consider an ephemeris-matched comparison or a sensitivity study around the AMITIS snapshot."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a serious look. The paper adapts the Stephenson electron test-particle code to cometary ions, adds energy-dependent ion-neutral cross sections and exothermic protonation, and drives the particles with AMITIS hybrid fields at 67P for 2.5-3 au. That specific combination is new, and the collision treatment is careful: the model is validated against the authors' 1D fluid model and against collisionless AMITIS output, and the comparison reproduces the earlier results while adding new species-level detail. The key result, a 5-10x underestimate of MIP/LAP densities in the terminator plane, is clearly presented and robust for this particular run. I also thought the discussion of why H2O+ sits near photochemical equilibrium while H3O+ and NH4+ pile up in the inner coma was informative.\n\nThe soft spots are real but not fatal. The main one is that the causal claim in the conclusions - that adiabatic electron closure drives the strong ambipolar field and over-transport - is not directly tested anywhere. No run varies gamma or replaces Eqs. 12-13 with a cooled electron treatment; the paper argues from prior work and from the size of the gap. Plausible, but not demonstrated. Second, the data comparison uses one AMITIS snapshot with fixed upstream solar wind and a uniform electron-impact boost, then bins MIP/LAP over a window of Q-nu and compares plane-averaged radial profiles rather than the actual spacecraft ephemeris. The 5-10x figure is therefore a statement about this model run, not a robust statistical envelope over the full data interval. The uniform boost could also hide radial production-rate errors that act in the same direction. These are caveats to the headline claim, not reasons to dismiss the paper. The model code itself is not released, though the AMITIS fields are, and the authors are open about the limitations of the hybrid electron treatment.\n\nWho gets value: cometary plasma modelers and anyone interpreting Rosetta plasma densities at intermediate outgassing. The paper consolidates a lot of prior work and gives a concrete target for what a kinetic, collisional electron treatment should fix. I'd send it to peer review, with the expectation that the closure issue is tested or the claim softened. I'd also cite the density underestimate and the species-resolved collision analysis. Bring it to reading group only if you want to debate the evidential weight of an inferred mechanism.","headline":"Adapts a collisional test-particle model to cometary ions, gets a clear 5-10x density underestimate vs Rosetta data, but the adiabatic-electron attribution is an inferred interpretation rather than a tested mechanism.","tokens_in":18257,"tokens_out":1387,"would_cite":true,"duration_ms":18897,"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":"A collisional ion model still underestimates Rosetta's plasma density at 67P by 5–10 times, pointing to the adiabatic electron closure in hybrid simulations.","keywords":["comet 67P","Rosetta","test-particle model","cometary ionosphere","hybrid simulation","ambipolar electric field","adiabatic electrons","ion-neutral collisions"],"falsifier":"Run the hybrid simulation with a non-adiabatic, collisional electron closure (e.g., a full electron energy equation with cooling) and compare the resulting ion bulk velocities and densities to the MIP/LAP data at 2.5–3 au; alternatively, directly measure the ion bulk speed in the inner coma with an ion composition analyser. If the data show bulk speeds below about 2 km/s or densities matching MIP/LAP, the adiabatic closure is the cause; if not, the representativeness of the snapshot or the ion production profile is the limiting factor.","tokens_in":1862,"feed_emoji":"☄️","tokens_out":1857,"duration_ms":50707,"temperature":0.7,"pith_summary":"The paper tries to establish that at intermediate outgassing of comet 67P (about 5.4×$10^{26}$ $s^{-1}$, at 2.5–3 au from the Sun), a 3D collisional test-particle model of cometary ions, driven by electric and magnetic fields from a hybrid simulation, produces plasma densities 5–10 times lower than those measured by Rosetta's MIP and LAP instruments. Even when ion–neutral collisions are included, the modelled ions are accelerated too quickly by the simulated ambipolar electric field, so the model loses too many ions through transport. The authors attribute this to the hybrid simulation's assumption of adiabatic electrons, which overestimates the ambipolar electric field strength near the nucleus. If correct, hybrid simulations with adiabatic electron pressure are insufficient to model the inner coma of weakly to moderately outgassing comets, and kinetic, collisional treatment of electrons is required to capture ion dynamics.","feed_headline":"Ion model undercounts comet plasma by 5–10 times","feed_subtitle":"At 2.5–3 au, test-particle densities fall below Rosetta data; adiabatic electron closure is blamed.","key_machinery":"The central object is a 3D collisional test-particle model, adapted from an earlier electron test-particle code, which follows H2O+, H3O+, and NH4+ macroparticles through electric and magnetic fields provided by a hybrid simulation. Ions are created with photoionisation (boosted to include electron-impact ionisation), pushed by the Lorentz force, and undergo ion–neutral collisions (proton transfer, momentum transfer, electron transfer) using energy-dependent cross sections. The hybrid simulation treats ions kinetically but electrons as a fluid with an adiabatic closure, p_e = n_e k_B T_e = $\\alpha$ n_e^gamma, which sets the ambipolar electric field that accelerates the ions. The test-particle model solves the ion trajectories on a 1000 km cubed grid with 25 km field resolution, and the comparison with Rosetta data uses the terminator plane only.","core_discovery":"The central claim is that at 67P's intermediate outgassing regime, the collisional test-particle model, using fields from a hybrid simulation with adiabatic electrons, underestimates the measured plasma density by a factor of 5–10. The modelled total ion density in the terminator plane stays well below the MIP/LAP data for the same outgassing window, while the field-free, chemistry-free model (which assumes ions move radially at the neutral speed) actually agrees well with the data below about 30 km. The modelled ions are accelerated to bulk speeds of 11 km/s at 100 km despite collisions, far above the neutral expansion speed of 0.7 km/s. The authors conclude that the adiabatic electron closure (Equations 12 and 13) produces too strong an ambipolar electric field, leading to excessive ion transport and low densities, and that a self-consistent collisional electron treatment is needed.","pith_inferences":["If the adiabatic closure is indeed the culprit, replacing it with a collisional electron model that allows electron-neutral cooling would flatten the ambipolar potential well, reduce ion speeds, and raise densities, potentially reconciling the model with MIP/LAP data without invoking different ion production rates.","The same test-particle approach could predict spatial maps of ion composition ratios (e.g., NH4+/H3O+) as diagnostics of where chemistry vs transport dominates, which Rosetta's ROSINA-DFMS measurements could test.","The finding implies that hybrid simulations of weakly outgassing comets using adiabatic electrons may systematically overestimate ion loss; previous results from such simulations should be re-examined.","A natural next step is to couple the collisional ion model with the collisional electron test-particle model to compute the ambipolar field self-consistently rather than from an adiabatic closure."],"forward_implications":["Hybrid simulations with adiabatic electron pressure are unlikely to reproduce the inner-coma ion densities at intermediate outgassing; a kinetic or collisional electron closure is required.","Even with collisions, the modelled ion bulk speeds are too high (11 km/s at 100 km), so any model relying on adiabatic electrons may systematically overestimate ion loss through transport.","The field-free, chemistry-free model, which assumes ions travel radially at the neutral speed, remains a good predictor of plasma density at low outgassing and below about 30 km at this outgassing.","Energy-dependent 3D cross sections boost proton transfer and NH4+ production compared to 1D kinetic rates, but the product of density and bulk velocity (n_i u_i) is preserved, so earlier electric-field strength estimates remain valid."],"supporting_citations":[{"why":"Supplies the original test-particle model code, adapted here for cometary ions, and the moment-calculation method.","marker":"Stephenson et al. (2022)"},{"why":"Provides the AMITIS hybrid simulation run with fields and ion densities used as input and validation for the test-particle model.","marker":"Moeslinger et al. (2024)"},{"why":"Provides the 1D ion acceleration model and diamagnetic-cavity case study against which the 3D collisional model is validated.","marker":"Lewis et al. (2024)"},{"why":"Establishes the transition in the field-free chemistry-free model's performance and the role of ambipolar electric field in enhanced ion transport.","marker":"Vigren et al. (2019)"},{"why":"Demonstrates that at low outgassing the field-free model explains plasma density, the baseline the 3D model is compared against.","marker":"Galand et al. (2016)"},{"why":"Provides the error estimate for neglecting dissociative recombination and the data-model comparison at low outgassing and post-perihelion.","marker":"Heritier et al. (2018)"},{"why":"Supplies the ionisation rates and the electron test-particle modelling that shows collisional cooling and the flattening of the potential well.","marker":"Stephenson et al. (2023)"},{"why":"Supports the finding that solving the electron pressure equation instead of using an adiabatic closure reduces the additional force on ions in the inner coma.","marker":"Koenders et al. (2015)"}],"fun_headline_variants":["Comet plasma undercount traced to electron closure","10x density miss at 67P blamed on adiabatic electrons","Rosetta data exposes ion model's electron flaw","At 67P, test-particle model needs collisional electrons","Ion transport overestimated in comet plasma model"],"cache_read_input_tokens":20352,"weakest_assumption_plain":"The paper compares one fixed simulation snapshot—with fixed upstream solar wind parameters, a constant boosted photoionisation rate, and adiabatic electrons—against Rosetta data binned over a range of outgassing conditions; if the true electric fields, solar wind conditions, or radial ion production profile during those intervals differ from this snapshot, the 5–10 times density underestimate could be caused by these factors rather than by the adiabatic electron closure.","fun_headline_variants_meta":{"raw":{"variants":["Comet plasma undercount traced to electron closure","10x density miss at 67P blamed on adiabatic electrons","Rosetta data exposes ion model's electron flaw","At 67P, test-particle model needs collisional electrons","Ion transport overestimated in comet plasma model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000162,"raw_usage":{"total_tokens":1292,"prompt_tokens":1049,"completion_tokens":243,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":665,"completion_tokens_details":{"reasoning_tokens":165}},"tokens_in":665,"tokens_out":243,"duration_ms":3454,"temperature":1.0,"reasoning_tokens":165,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:39:22.373207+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the hybrid simulation with a non-adiabatic, collisional electron closure (e.g., a full electron energy equation with cooling) and compare the resulting ion bulk velocities and densities to the MIP/LAP data at 2.5–3 au; alternatively, directly measure the ion bulk speed in the inner coma with an ion composition analyser. If the data show bulk speeds below about 2 km/s or densities matching MIP/LAP, the adiabatic closure is the cause; if not, the representativeness of the snapshot or the ion production profile is the limiting factor.","supporting_citations":[],"review_version":1}