{"id":"9c2e5285-32f5-492f-aaef-7a95797c80fe","arxiv_id":"2602.07452","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A new general-relativistic particle-in-cell code, FPIC, is validated on black-hole magnetosphere tests and reproduces the Blandford–Znajek power with a one-parameter fit.","lead":"This paper presents FPIC, a new computer code for simulating plasmas around black holes using particle-in-cell methods in curved spacetime, including a hybrid particle-mover that switches between fast and accurate integrators. It reproduces known black-hole magnetosphere effects and matches the predicted Blandford–Znajek jet power, making it a resource for modeling black-hole jets.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Simplified GJ-based pair injection could contaminate the split-monopole BZ luminosity; an injection-parameter sweep would settle it.","rationale":"The reader's weakest_assumption is the simplified pair-plasma injection scheme, and this is indeed the most load-bearing concern. The BZ luminosity is the headline validation result, and it is computed from electromagnetic fields that the injection directly affects. The paper is transparent about the approximation but does not provide a sensitivity study or evidence that results are independent of injection choices. The fitted κ strengthens the concern because it can absorb a systematic normalization error, weakening the claim of 'very good agreement with analytical predictions.' I considered other issues: the one-parameter fit is a real limitation but is mitigated by the fact that the spin dependence of the curve is still tested; the Penrose-process claim from negative-energy particles lacks an energy-flux budget but is secondary to the BZ validation; the lack of code release is a reproducibility concern but not a correctness flaw. Injection sensitivity is the most direct threat to the central claim. The recommended verdict remains CONDITIONAL (unchanged), because the concern is significant enough to require additional evidence but does not warrant rejection absent evidence of an actual failure. The proposed test is a standard parameter sweep that can be performed with the existing code and would settle whether the injection is load-bearing.","tokens_in":30358,"tokens_out":4612,"duration_ms":51126,"concrete_test":"Re-run the split-monopole simulation at a*=0.9 (the resolution-test case) with all parameters fixed except the injection prescription: vary multiplicity M over {1, 3, 10, 30, 100}, temperature Θ over {0.1, 0.5, 2.0}, and injection cadence Δt_inj over {0.005M, 0.01M, 0.02M}. For each run, compute the time-averaged P_BZ/P_maxBZ0 over 5≤t̄/M≤40 and the fitted κ. If the BZ power and κ vary by less than ~6% (the reported resolution error) across this sweep, the injection scheme is not determinative; if the power shifts by more than that, the central claim is injection-dependent and the BZ validation is not yet robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that FPIC reproduces the BZ luminosity in the split-monopole configuration—depends on the pair-plasma injection scheme (§3.4). Particles are injected only where n < M n_GJ, with M=10, temperature Θ=0.5, and a fixed injection cadence Δt_inj=0.01M; the Wald case uses M=3 and n0=5n_GJ. This prescription is an ad hoc source that supplies plasma to screen electric fields, but it is not derived from microphysical pair-production or transport. The paper itself concedes in §4 that 'more accurate particle-injection strategies need to be developed and tested.' If the injection rate or multiplicity were too high or too low, the magnetosphere could become artificially force-free or charge-starved, changing the toroidal field and Poynting flux that are integrated to compute P_BZ in Eq. (49). Moreover, the reported agreement with the analytic formula (51) uses a one-parameter fit for κ=0.041; a constant multiplicative error in the simulated Poynting flux would be absorbed by κ, so the shape agreement does not validate the absolute luminosity. The claim in §4 that results are 'consistent' with more detailed injection schemes (Parfrey et al. 2019; El Mellah et al. 2022; Chen et al. 2025) is asserted but not demonstrated in this paper. Thus the central BZ validation is load-bearing on the unverified injection prescription.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents FPIC, a new general-relativistic particle-in-cell (GRPIC) code for stationary, axisymmetric black-hole spacetimes. The code solves Maxwell's equations on a spherical Kerr-Schild grid with a Yee/FDTD scheme, deposits charges/currents with metric-weighted volume interpolation, and evolves particles with RK4, implicit midpoint, Hamiltonian, and a new hybrid integrator. The numerical methods are described in detail and validated on neutral and charged particle orbits, the vacuum Wald solution, the plasma-filled Wald configuration, and a split-monopole magnetosphere. The central claims are that the hybrid pusher improves energy conservation at reduced cost, that the plasma-filled Wald simulation shows evidence for the Penrose process through negative-energy-at-infinity particles, and that the split-monopole simulations reproduce the Blandford-Znajek luminosity in agreement with analytic predictions.","tokens_in":30721,"tokens_out":5004,"duration_ms":55737,"significance":"If fully substantiated, FPIC would be a valuable new GRPIC tool for axisymmetric black-hole magnetospheres. The paper contains genuine strengths: independent test-particle benchmarks against known geodesics and Wald-field orbits, a converged vacuum Wald comparison, a resolution study for the split-monopole Poynting flux, and a very detailed description of the numerical algorithms that will help reproducibility. The hybrid integrator idea is interesting and appears to give useful speedups in the neutral-particle tests. However, the headline BZ and Penrose claims are currently supported by weaker evidence than the abstract suggests: the BZ comparison uses a fitted normalization, the pair-plasma injection is an ad hoc prescription with no demonstrated insensitivity, and the Penrose claim is based on negative-energy particles rather than a measured energy flux. These are load-bearing issues for the manuscript's central validation narrative.","major_comments":[{"comment":"The abstract and Sec. 4 state that FPIC 'successfully reproduce[s] the Blandford-Znajek luminosity.' The comparison is made after a one-parameter fit of κ=0.041 in Eq. (50) to the simulation data. A fitted normalization validates the spin dependence F(Ω_h) but cannot independently validate the absolute luminosity; any overall multiplicative error in the computed Poynting flux is absorbed by κ. Please either constrain κ from the magnetic topology, report the fit uncertainty and explicitly limit the claim to shape agreement, or provide an independent normalization test. As it stands, the 'very good agreement with analytical predictions' is overstated.","section":"Sec. 3.4, Eqs. (50)-(51), Fig. 10"},{"comment":"The plasma-filled simulations use an ad hoc pair-injection prescription: particles are added where n < M n_GJ, with M=10 (split monopole) or M=3 (Wald), a Maxwell-Jüttner temperature Θ=0.5, and a fixed injection cadence Δt_inj=0.01M. No parameter sweep is presented, and Sec. 4 itself concedes that 'more accurate particle-injection strategies need to be developed and tested.' Because the BZ power and current-sheet structure depend on the plasma supply, an injection rate that is too high or too low could change the toroidal field and hence P_BZ in Eq. (49). The asserted consistency with more detailed injection models (Parfrey et al. 2019; El Mellah et al. 2022; Chen et al. 2025) is not demonstrated in this paper. A sensitivity study over M, Θ, and injection cadence, or a direct comparison with a more physical injection model, is needed before the BZ agreement can be attributed to the magne","section":"Secs. 3.3-3.4, injection scheme"},{"comment":"The presence of electrons with negative energy at infinity, ⟨e_∞⟩<0, inside the ergosphere is necessary but not sufficient to establish that the Penrose process is active. The paper does not show that these negative-energy particles actually cross the event horizon, nor does it measure a net outward energy flux associated with their absorption. A flux diagnostic, e.g., the horizon-integrated energy flux carried by particles with e_∞<0, should be provided before claiming that 'the Penrose process is active' in the plasma-filled Wald simulation. Without this, the statement remains suggestive rather than demonstrated.","section":"Sec. 3.3, Penrose process claim"},{"comment":"The hybrid integrator is validated only for neutral test particles in pure geodesic motion. In a full PIC code, interpolation of electromagnetic fields is expected to dominate the energy error, as the charged-particle tests in Sec. 3.1.3 indicate. To support the claim that the hybrid approach 'guarantees high precision at comparatively small computational costs' in the intended application, it should be demonstrated on charged-particle trajectories with the same field-interpolation scheme used in the PIC runs, or at least on a representative GRPIC configuration. Without such a test, the practical benefit of the hybrid method for FPIC is not established.","section":"Sec. 3.1.2, Figs. 4-5"}],"minor_comments":[{"comment":"Grammar: 'a code built in this way, i.e., FPIC is to reproduce' appears to be missing 'able'; please revise to 'is able to reproduce.'","section":"Abstract and Sec. 4"},{"comment":"The caption says 'iθ in the radial direction'; this should be the polar/angular direction.","section":"Fig. 1 caption"},{"comment":"The divergence-cleaning procedure is described as 'periodically' solving the Poisson equation every 25 timesteps with 500 Jacobi iterations. Please state the convergence criterion used for the Jacobi solver, or clarify why a fixed iteration count is sufficient.","section":"Sec. 2.4"},{"comment":"The text says adaptive timesteps are employed in the hybrid schemes, but Table C.2 lists fixed values of Δt for RK4 and Hamiltonian steps. Please clarify how the adaptive step is determined and whether the listed values are upper bounds.","section":"Table C.2 and Sec. 3.1.2"},{"comment":"Reference 'Pierre Jacques et al. 2025' appears incomplete (missing journal/page data). Please check the reference list entries for consistency.","section":"Sec. 1"},{"comment":"Given the stated goal of reproducibility, consider adding a statement on code availability or a link to a public repository, even if only a limited distribution is planned.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a solid code-description paper with useful validation tests, but the abstract and conclusions claim more than the current evidence supports. The BZ luminosity comparison is a one-parameter fit to the simulation data, and the pair-plasma injection is an uncontrolled ad hoc source in the same simulations. The Penrose-process claim similarly rests on a necessary-but-not-sufficient diagnostic. These issues are fixable with additional targeted tests and softened claims, so I do not recommend rejection, but the paper should not be accepted in its present form. The companion paper Meringolo et al. (2025) may already contain some of the missing injection or flux analysis; if so, the present paper should either include those results or cite them explicitly for the central claims."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Good news first: this is a real code paper with a genuinely new piece of algorithm. The hybrid pusher that switches between RK4 and the Hamiltonian integrator based on a monitored energy error is clever and clearly demonstrated on the geodesic testbed; it buys better energy conservation at lower cost, and the paper documents the thresholds and timings. The description of the grid, deposition, divergence cleaning, and boundary treatments is detailed enough to be useful for anyone building or extending a GRPIC code. The validation is generally honest: known orbits, the Wald electrovacuum solution, and a resolution test on the BZ power all behave as they should.\n\nThe soft spots are in the astrophysical claims, not the numerics. The BZ comparison in Fig. 10 fits one global constant kappa. With that, the shape match over spin is a good consistency test, but it is not a parameter-free prediction of the luminosity. The paper calls it a 'very strong validation'; I would say it is a solid check, not a proof. The pair-plasma injection is simplified GJ-based with M=10, and the paper itself concedes more accurate injection is needed. The claim that results are consistent with Parfrey et al., El Mellah et al., and Chen et al. is asserted, not shown. Since the toroidal field and Poynting flux can depend on the plasma supply, an injection-parameter sweep would settle whether the BZ power is robust. I'd like to see that before the BZ result is used as a benchmark. The Penrose-process statement is also a bit ahead of the evidence: negative-energy particles inside the ergosphere 'indicate' the process is active, but no net energy flux across the horizon is computed. Finally, the hybrid pusher is only tested for neutral particles; its behavior for charged particles in a PIC environment is untested. And despite the stated reproducibility goal, no code or data is released, which is a real gap for a methods paper.\n\nAll that said, none of this is fatal. The code is carefully built, the benchmarks are appropriate, and the hybrid integrator is a worthwhile contribution. I'd send this to a referee. The revision should add a short injection study, soften the BZ and Penrose language, and either release the code or explain why not. This is a paper for people building GRPIC codes and for black-hole magnetosphere simulators who want a second independent code. Worth engaging, not the last word.","headline":"Serious code paper with a genuinely new hybrid pusher; the BZ validation is solid but not parameter-free, and the Penrose-process claim outruns the evidence.","tokens_in":31141,"tokens_out":2568,"would_cite":true,"duration_ms":27608,"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 new general-relativistic particle-in-cell code, FPIC, tracks collisionless plasma near spinning black holes and reproduces the analytic Blandford-Znajek luminosity for the split-monopole magnetosphere.","keywords":["GRPIC","particle-in-cell","Kerr black holes","Blandford-Znajek process","Penrose process","Kerr-Schild coordinates","black-hole magnetospheres","kinetic plasma"],"falsifier":"Run the same split-monopole setup with a different plasma injection prescription—for example, a pair-cascade-based source rate, a different multiplicity, or a different temperature—and check whether the measured Blandford-Znajek power as a function of spin still collapses onto the analytic curve; alternatively, stop particle injection after a steady state and see whether the negative-energy electrons inside the ergosphere persist or decay away.","tokens_in":30238,"feed_emoji":"🕳️","tokens_out":3347,"duration_ms":37220,"temperature":0.7,"pith_summary":"FPIC is a new general-relativistic particle-in-cell code that models collisionless plasma around stationary, axisymmetric black holes. The paper argues that FPIC can reproduce fully nonlinear kinetic plasma dynamics near rotating black holes, supported by two astrophysical tests. In a plasma-filled Wald magnetosphere it finds electrons with negative energy at infinity inside the ergosphere, a signature that the Penrose process is active. In a split-monopole magnetosphere it measures the Blandford-Znajek power across black-hole spins and matches a high-order analytic expression using a single topology-dependent coefficient. The paper also introduces a hybrid particle pusher that switches between RK4 and an energy-conserving Hamiltonian integrator, improving energy conservation at lower computational cost.","feed_headline":"Plasma code reproduces black-hole jet power","feed_subtitle":"FPIC simulates the plasma around spinning black holes and matches the analytic Blandford-Znajek luminosity.","key_machinery":"The argument rests on a few components: spherical Kerr-Schild coordinates, which keep the metric regular at the event horizon and encode the spherical topology; a finite-difference time-domain Yee-grid Maxwell solver with divergence cleaning; a set of particle pushers, including a novel hybrid scheme that monitors the violation of the Hamiltonian energy and dynamically switches between the fast RK4 integrator and the energy-conserving Hamiltonian integrator, using adaptive timesteps; and a volume-weighted charge and current deposition that dominates the interpolation error. For the astrophysical tests, the key comparison is the measured Poynting flux through a sphere near the horizon against","core_discovery":"The central claim is that FPIC is capable of reproducing correctly and accurately fully nonlinear plasma dynamics near rotating black holes. The evidence is a series of validations: neutral and charged particle trajectories conserve Hamiltonian energy well; the vacuum Wald solution relaxes to the expected rotating configuration with the Meissner expulsion of magnetic-field lines; plasma-filled Wald simulations show negative-energy-at-infinity electrons inside the ergosphere, indicating an active Penrose process; and split-monopole simulations produce a Blandford-Znajek luminosity that agrees with analytic high-order predictions after fitting one spin-independent coefficient tied to the magne","pith_inferences":["The simplified plasma injection scheme, based on the local Goldreich-Julian density with chosen multiplicity and temperature, is the main uncertainty; varying the injection prescription could test whether the reported Penrose particles and Blandford-Znajek power are robust or artifacts of injection.","Because the code is axisymmetric (2.5D), three-dimensional instabilities and non-axisymmetric modes are excluded; extending to 3D could change plasmoid dynamics and the quantitative Blandford-Znajek agreement.","The hybrid-integrator principle, selecting a scheme based on a monitored Hamiltonian error, could be applied to other Hamiltonian systems, including those with radiation reaction or pair-production terms.","The successful one-parameter match to the perturbative Blandford-Znajek formula suggests that PBZ versus spin could serve as a standard quantitative cross-code benchmark in GRPIC."],"forward_implications":["If the central claim holds, FPIC provides a validated path to kinetic, first-principles modeling of black-hole magnetospheres, capturing microphysics that GRMHD cannot.","The hybrid integrator offers a practical way to improve energy conservation in particle pushers without paying the full cost of an implicit Hamiltonian scheme.","The measured Blandford-Znajek power for the split monopole, matching analytics over a range of spins, provides a strong benchmark for future GRPIC codes.","The plasma-filled Wald simulations indicate that Penrose-process negative-energy particles can be produced self-consistently in an ergospheric current sheet.","The detailed method description supports reproducibility and enables others to build or verify similar codes."],"fun_headline_variants":["Black-hole jets simulated to match theory","New PIC code reproduces black-hole jet power","FPIC: plasma code matches black-hole jet luminosity","Simulating black-hole plasma: FPIC hits jet-power target"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The plasma-filled results assume that injecting electron-positron pairs according to the local Goldreich-Julian density, with a fixed multiplicity and temperature, faithfully represents the physical particle supply near the black hole; if this injection is not realistic, the negative-energy Penrose particles and the Blandford-Znajek power could be artifacts of the injection rather than genuine magnetospheric physics.","fun_headline_variants_meta":{"raw":{"variants":["Black-hole jets simulated to match theory","New PIC code reproduces black-hole jet power","FPIC: plasma code matches black-hole jet luminosity","Simulating black-hole plasma: FPIC hits jet-power target"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000298,"raw_usage":{"total_tokens":1581,"prompt_tokens":786,"completion_tokens":795,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":530,"completion_tokens_details":{"reasoning_tokens":734}},"tokens_in":530,"tokens_out":795,"duration_ms":8165,"temperature":1.0,"reasoning_tokens":734,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T03:34:56.371052+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same split-monopole setup with a different plasma injection prescription—for example, a pair-cascade-based source rate, a different multiplicity, or a different temperature—and check whether the measured Blandford-Znajek power as a function of spin still collapses onto the analytic curve; alternatively, stop particle injection after a steady state and see whether the negative-energy electrons inside the ergosphere persist or decay away.","supporting_citations":[],"review_version":1}