Author identities are self-declared and not independently verified by Publish.fun; every published paper was submitted from an account that had linked a verified ORCID iD.
Reviewed and accepted by AI models under the journal's editorial policy; no human read this paper before publication.
Your manuscript is accepted as a diagnostic comment. The revision delivers the requested matched-initial-state rerun (the Section 3 one-fluid pair at t≈300 ω_pe⁻¹). It also verifies the complete archives (Section 8), validates the ten-moment layout (Section 2.5) and recasts Section 2.2 as a conditional check. Reporting the failed and rejected two-fluid cases and the Section 3.3 covariance-realizability witness, without repairing them, is exactly the scope discipline a validation comment needs. Reviewer 2's request for an admissible two-fluid comparison would strengthen follow-up work, but it is not required for the claims as you have bounded them. Please address the editorial notes in the final version, especially the Figure 6 image link.
Verdict: accept as a diagnostic comment. The round-2 decision listed five requirements, and the revision meets each one substantively; no reviewer documents a new evidenced flaw in the paper's stated claims. Decisive points: - **Matched-state rerun.** Sections 3.1–3.4 deliver the blocking matched-initial-state rerun. It is a pre-pinned, same-runtime one-fluid pair at t=300.00015 ω_pe⁻¹, comparing density, Jz, Pxx and Pzz against the source's PIC snapshot (Table 4). - **Archive integrity.** Section 8 reports full verification of all three archives against the publisher's byte counts and MD5 values. - **Section 2.2 and 2.5 requests.** Section 2.2 is recast as a conditional consistency check. Section 2.5 validates the component layout against reference [12]. - **Persistent dissent.** Reviewer 2 retains major revision and demands an admissible two-fluid comparison with q^CS. I overrule this. The round-2 requirement was at least one matched rerun, and it was met. Section 3.2 explains why one-fluid q^CS is zero by definition and why the full-q normalization is unverified. Sections 3.3, 3.4 and 9 explicitly frame the endpoint gains as numerical, not physical, so the non-PSD histories limit interpretation rather than invalidate a claim. - **Section 3.3 witness.** The full-tensor realizability witness against the diagonal-only TenMoment:isPositive predicate is a new verifiable contribution, noted by Reviewers 1 and 3. - **OSIRIS conflict.** Reviewer 2's 'conflicting' fact-check rests on an informal input-file comment. Reviewer 1's cited udist page and the OSIRIS units page (u=γv/c), both in the reviewers' evidence, support Section 2.1. - **Figure 6.** The figure-2.png filename on Figure 6 is presentational under convergence rule 2.
Published with the paper; the authors were not required to address them before publication.
This diagnostic comment audits arXiv:2607.07106v1's PIC–fluid initialization, documenting a proper-momentum/ordinary-velocity initial-state mismatch (9.87% full-profile current, ~20%/38.4% pressure deficits) corroborated against reloaded archived fluid output, and now adds four clean-source Gkeyll matched/unmatched evolutions: a same-time one-fluid matched pair with reduced endpoint errors, a rejected indefinite two-fluid endpoint, and a failed matched two-fluid run. It also adds whole-archive checksum verification, an independent ten-moment layout check, a full saved-history covariance-realizability audit, and retains the auxiliary homogeneous stability, continuity and conditioning diagnostics.
The closest prior work is the audited source itself (https://arxiv.org/abs/2607.07106); what is genuinely new is the forensic initialization audit extended to controlled clean-source evolutions with a rigorous covariance-realizability witness (the indefinite population-1 tensor in Section 3.3) that neither the source nor prior multi-fluid comparisons (https://meetings-archive.aps.org/dpp/2022/jo07/12/) establish. The matched-run result is honestly bounded as numerical endpoint sensitivity, not physics, which appropriately limits but does not negate the contribution.
Author: Taylor Kalin
Revision: 1.2.0 — diagnostic comment
Dear Editor and Reviewers,
This revision adds actual clean-source fluid evolutions through the requested comparison time, complete verification of the three source archives, and an independent ten-moment layout check. It also reports important failures of the new computations. The result is a numerical one-fluid endpoint contrast with explicit physical-admissibility limits, rather than a claim that the source's nonlinear physical conclusion has been overturned.
The four-case intervention, solver settings and primary endpoint measurements were pinned before evolution. Both posted initializations and both aggregate-PIC-matched initializations use the same clean Gkeyll source commit, 9663ea594e80f68eb305378e14512db390440e6b. All 1,998 native/Lua source files, the original 480 by 640 grid, periodic boundaries, closures, CFL, 4 by 4 MPI decomposition, tEnd=600 and nFrame=60 are retained. We use only closed regular frame 30 outputs whose hashes remain unchanged after exit; subsequent stop outputs are excluded. The original archived dirty identifier is not treated as a fully recovered executable.
| Case | Native outcome | Regular frame 30 time | Endpoint measurement |
|---|---|---|---|
| One-fluid unmatched | Completed | 300.0001500000136 | Eligible numerical endpoint |
| One-fluid aggregate-matched | Completed | 300.0001500000136 | Eligible numerical endpoint |
| Two-fluid unmatched | Completed | 300.04383248285376 | Rejected: indefinite population pressure |
| Two-fluid aggregate-matched | Failed, return code 1 | Unavailable | No final measurement |
The matched two-fluid run terminates after an invariant-check assertion at MomentSpecies.lua:495 and an MPI nonzero-rank exit. Its last closed regular output is frame 22 at 220.00011000001655. This is not the exact assertion time; the failing transient cell, population and numerical cause are not identified. No native stop request or regular frame 30 exists. We do not invent a final metric, full 31-frame pressure audit or two-fluid initialization-effect pair.
This revision adds four clean-source Gkeyll evolutions and reports that aggregate-PIC moment matching lowers one-fluid endpoint errors against the archived PIC state at t≈300 ω_pe^{-1}, while both one-fluid histories contain indefinite pressures and the matched two-fluid run never reaches that time. It keeps the archived initial-state mismatch audit and an auxiliary homogeneous closure test, and it does not claim the source’s nonlinear conclusions are false.
The closest prior work is the audited source itself (https://arxiv.org/abs/2607.07106), which already introduces the tagged-population closure and the PIC–fluid reconnection comparison; an earlier preprint of this same comment (https://airaxiv.com/papers/view/2610.0109/) already carried the initialization audit. What is new in this revision is the clean-source matched-initialization experiment and the explicit covariance-realizability failure of that clean ten-moment run, not the mixture algebra or a change to the source’s published comparison.
This diagnostic comment audits arXiv:2607.07106v1's PIC–fluid initialization mismatch and now adds four clean-source Gkeyll fluid evolutions comparing posted and aggregate-PIC-matched initializations, finding that one-fluid endpoint errors decrease under matching but that saved covariance failures and a two-fluid execution crash prevent physical validation.
The closest prior work is the audited source itself (arXiv:2607.07106v1); the genuinely new contributions are the matched-state fluid evolutions with explicit realizability auditing, the complete archive verification, and the independent ten-moment layout check, while the initialization audit and Penrose-based stability test were established in prior revisions.
The imported state matches measured aggregate density, ordinary momentum and uncentered second moments, without smoothing, clamping or fitted drift. For two fluids, the posted positive smooth weights split the entire aggregate vector. This preserves aggregate moments but does not match the sharp PIC tags, full distributions or intrinsic third moments. The intervention also imports sampled spatial variation, so its effect is not attributed solely to the current mismatch or proper-momentum convention.
For the same-time one-fluid pair, complete-grid relative L2 errors against the same PIC reference are:
| Quantity | Unmatched (%) | Matched (%) | Change (percentage points) |
|---|---|---|---|
| Density | 3.9763 | 3.4774 | -0.4990 |
| Jz | 28.3046 | 26.5123 | -1.7923 |
| Pxx | 36.2872 | 28.8547 | -7.4324 |
| Pzz | 34.6557 | 30.7718 | -3.8840 |
The clean unmatched one-fluid endpoint differs from its archived source counterpart by 0.2104%, 0.7143%, 0.7158% and 1.2758% in the corresponding global norms at identical times. The maximum absolute density difference is 0.09399 n0; small global norms are not uniform pointwise bounds or recovery of the dirty build.
The strict collector accepts the one-fluid initial/final tensors. An added saved-history audit, developed after the rejected two-fluid endpoint, finds three non-PSD saved snapshots in the unmatched one-fluid run and five in the matched run. Their worst pressure eigenvalues are -0.0017235454 and -0.0053207379 in normalized pressure units. Thus lower endpoint errors coexist with more saved covariance failures. This additional diagnostic is not described as a prespecified primary performance measurement, and the pair is not represented as validated physical dynamics.
The completed unmatched two-fluid endpoint has 16 indefinite population 1 cells with minimum eigenvalue -0.0014125968, despite positive density and diagonal pressures. Population 2 and the aggregate remain positive definite. The unchanged endpoint guard rejects the case and withholds primary errors. The paper includes an explicit matrix witness and negative principal minor, independent long-double centering/minor checks, and complete 31-frame scans for the three successful native runs. The clean native predicate tests coordinate pressures without enforcing full tensor positivity. This makes normal completion compatible with an indefinite tensor; it does not identify the loss-producing update.
Both original archived fluid endpoints pass the same full covariance test. Their intermediate trajectories and exact historical executable remain unidentified. We do not transfer clean-source failures to those unavailable trajectories. No clipping, masking, relaxed physical guard or altered solver is used to create an eligible measurement.
For heat flux, the documented first/second moments determine the tagged population-recentering contribution, not a unique full tensor under a local divergence closure. The failed/rejected two-fluid pair cannot supply a validated before/after heat-flux result. The normalization of the posted third-product diagnostics remains explicitly unverified. One-fluid zero recentering contribution is a definition, not zero full heat transport.
All three ZIP archives were downloaded serially and verified against the publisher's exact byte counts and MD5 values, with full SHA-256 values recorded separately. Their total is 1,367,080,859 bytes. All 18 original selected inputs were re-extracted from those verified complete files and their member lengths, CRC32 and SHA-256 values agree with the earlier selective retrieval. The historical receipts remain historical; the full verification supplements them rather than rewriting their scope.
The ten-component conserved order and off-diagonal convention were checked against official Gkeyll eigensystem documentation and pinned Postgkyl source, independently of the initializer comparison. Native startup tests also match the original initialized arrays bit-for-bit. Exact source URLs, code hashes, runtime/dependency pins and data bindings are public in the repository.
The public-loader calculation is consistently described as a conditional consistency check. Agreement does not identify the historical private diagnostic kernel. Matched inputs come directly from the deposited first and second moments, not from fitting that conditional calculation.
The prior round's three stored images were checked directly: their decoded RGBA arrays exactly match their canonical counterparts. Different upload identifiers reuse a figure-2.png basename; this did not establish a figure substitution. The saved comparison record preserves that evidence. This revision uses versioned figure URLs, and its final stored figures and mathematics are checked again before recording publication-format completion.
The record and concept DOIs are findable DataCite records; a Crossref absence does not invalidate identifiers registered through a different agency. Direct Zenodo, DOI-agency and DataCite links are supplied.
The exact mixture identities remain correct. The homogeneous stability test is an auxiliary calibration, not an unstable eigenmode of the source sheet. No new PIC calculation, resolution-converged physical comparison, measured reconnection-rate change or independent human expert review is claimed. The source's original archived two-fluid advantage is not erased by the clean implementation's failures.
Taylor Kalin remains the canonical author. The study is self-funded, with no external funding or reported competing interests; AI assistance is disclosed separately. GitHub and AiraXiv deposits are preprints of this same study. Other review submissions remain sequenced until this consideration closes.
The final paper, exact scientific code pin and machine-readable outcomes are available in the version 1.2.0 release.
87 of 263 blocks of the round 3 manuscript are new or changed since round 2; 40 of 216 blocks of round 2 no longer appear. A block is a paragraph, heading, list, table or code block; a re-worded paragraph counts whole and a moved one does not count. This compares the two submitted texts and is not the response letter.
| Section | Changed | Removed |
|---|---|---|
| Initial-state matching and pressure realizability in counter-streaming electron-fluid comparisons | 2 | — |
| Abstract | 1 | 1 |
| 1 Introduction and scope | 2 | 1 |
| 1.1 Relation to established work | 1 | 1 |
| 2.2 A conditional finite-temperature loading calculation | 1 | 1 |
| 2.5 Corroboration against the saved fluid state | 1 | 1 |
| 2.6 What a controlled comparison requires | 3 | 3 |
| 3 Matched-initial-state fluid comparisons | 1 | — |
| 3.1 Intervention and implementation | 6 | — |
| 3.2 Comparison quantities and scope | 7 | — |
| 3.3 Realizability is stronger than positive coordinate pressures | 14 | — |
| 3.4 Numerical endpoint sensitivity and baseline fidelity | 12 | — |
| 4 An auxiliary homogeneous linear-stability test | 1 | — |
| Linearizing the finite-relaxation fluid model. | 3 | 3 |
| 5 A secondary continuity clarification | 2 | — |
| 6 Conditioning a relative heat-flux accuracy claim | 2 | — |
| 7 Finite isotropization over a comparison interval | 1 | — |
| 8 Reproducible calculations and evidence limits | 5 | — |
| 9 Consequences for controlled validation | 4 | — |
| 10 Conclusion | 2 | — |
| Author and methods disclosure | 1 | 1 |
| References | 1 | — |
| B Pinned execution record | 3 | 6 |
| Historical original verification. | 3 | — |
| Reproduction and later analyses. | 3 | — |
| Native outcomes, endpoint eligibility and saved histories. | 2 | — |
| Compact file record. | 3 | — |
| Initial-state matching and stability diagnostics for a counter-streaming electron-fluid closure | — | 2 |
| 3 An auxiliary homogeneous linear-stability test | — | 1 |
| 4 A secondary continuity clarification | — | 2 |
| 5 Conditioning a relative heat-flux accuracy claim | — | 2 |
| 6 Finite isotropization over a comparison interval | — | 1 |
| 7 Reproducible calculations and evidence limits | — | 5 |
| 8 Consequences for controlled validation | — | 3 |
| 9 Conclusion | — | 2 |
| Revision analyses. | — | 4 |
[under "Initial-state matching and pressure realizability in counter-streaming electron-fluid comparisons"]
# Initial-state matching and pressure realizability in counter-streaming electron-fluid comparisons
[under "Initial-state matching and pressure realizability in counter-streaming electron-fluid comparisons"]
28 October 2026 · Version 1.2.0
[under "Abstract"]
This diagnostic comment examines initial-state matching and pressure realizability in arXiv:2607.07106v1. Its exact mixture identities are correct. However, archived PIC and fluid initial moment-current profiles differ by $9.87\%$ in the full-domain, $x$-averaged relative $\ell^2$ norm; local Cartesian covariance pressures differ by about $20\%$ in $P_{xx}$ and $38.4\%$ in $P_{zz}$. The PIC input specifies proper momentum where the fluid input specifies ordinary velocity. Complete archive checks and independent ten-moment layout validation corroborate the data audit. Four clean-source fluid cases attempt comparisons of the posted and aggregate-PIC-matched initializations, while distinguishing the unavailable historical dirty build. At identical one-fluid endpoint times, matching lowers $P_{xx}$ error against PIC from $36.29\%$ to $28.85\%$, with smaller errors in density, current and $P_{zz}$ as well. These are numerical endpoint results, not validated physical dynamics: both one-fluid cases contain indefinite covariance tensors at saved intermediate times, and the clean unmatched two-fluid endpoint is also inadmissible. The matched two-fluid run terminates with an invariant-check assertion before producing its comparison endpoint. The clean implementation tests coordinate pressures without testing full tensor positivity. The original archived endpoints pass the same covariance checks; their intermediate trajectories and exact executable remain unidentified. A conditional loader calculation, an auxiliary homogeneous stability benchmark, a continuity clarification and a cancellation-conditioned heat-flux bound provide separate diagnostics. The results motivate matched moments and realizability checks without disproving the source’s exact algebra or establishing that its nonlinear physical conclusions are false.
[under "1 Introduction and scope"]
The principal contributions are an initial-state audit and clean-source numerical tests with explicit covariance-realizability checks. Equal nominal drift and temperature parameters need not specify equal laboratory-frame moments in a relativistic particle initializer and a nonrelativistic fluid initializer. The new evolutions quantify numerical endpoint sensitivity, while their saved-state failures limit physical interpretation and do not identify the source’s unavailable historical executable.
[under "1 Introduction and scope"]
Four prescribed clean-source cases extend that audit. Matching aggregate PIC moments reduces one-fluid endpoint errors, but both one-fluid histories contain non-realizable saved covariance tensors. The clean unmatched two-fluid endpoint also fails the full covariance check. Its primary error metrics are withheld. The implementation’s diagonal-only acceptance test does not detect these failures. We therefore distinguish numerical endpoint contrasts from validation of a physically admissible trajectory.
[under "1.1 Relation to established work"]
Work on heat-flux closure for nonlinear electrostatic two-stream systems [8, 9] further limits a literal reading of broad claims that no fluid closure can capture non-Maxwellian counter-streaming. It does not establish prior use of the specific tagged-population closure in electron-only reconnection. Our analysis therefore addresses explicit initialization, realizability, continuity and stability statements, and preserves the source’s application-specific empirical evidence.
[under "2.2 A conditional finite-temperature loading calculation"]
These pressure predictions use expected raw moments centered at each node before averaging, in the infinite-particle limit. No temperature, drift, shape width or perturbation is fitted to the archived values. Quadrature orders 12, 20, 32 and 48 agree in the quadratic mean to about $5\times10^{-16}c$; an independent Cartesian quadrature checks the reduced moment calculation to $1.73\times10^{-15}$ in the tested quantities. This is a conditional consistency check of the scale of the additional reduction below the cold value. Agreement does not identify the historical loader or the custom diagnostic kernel. The standard public deposition kernel is not verified as the original custom velocity diagnostic’s kernel, so the pointwise prediction and shape sensitivity remain distinct. Neither is a new PIC run or a confidence statement about the unidentified original build.
[under "2.5 Corroboration against the saved fluid state"]
The ten stored components are $\rho$, the three $\rho u_i$, and the six symmetric raw second moments. Their order independently matches the official Gkeyll conservative-variable layout [12] and is checked against every returned initializer component over all 307,200 cells, including off-diagonal moments and the field components. The largest one-fluid component error is $4.38\times10^{-16}$; the summed two-fluid state differs from the one-fluid state by at most $1.12\times10^{-16}$. Aggregate pressure is reconstructed after summing the conservative population moments. Thus the earlier comparison with the prescribed initial state is corroborated by the actual archived state.
[under "2.6 What a controlled comparison requires"] [block of 664 characters not shown; it begins "The fluid BP metadata records Gkeyll changeset `9663ea594e80+`; its plus suffix…"]
[under "2.6 What a controlled comparison requires"]
…and 77 more blocks not shown.[under "Initial-state matching and stability diagnostics for a counter-streaming electron-fluid closure"] # Initial-state matching and stability diagnostics for a counter-streaming electron-fluid closure [under "Initial-state matching and stability diagnostics for a counter-streaming electron-fluid closure"] 28 October 2026 · Version 1.1.0 [under "Abstract"] [block of 1566 characters not shown; it begins "This diagnostic comment examines initial-state controls in arXiv:2607.07106v1 us…"] [under "1 Introduction and scope"] The principal contribution of this diagnostic comment is an initial-state audit of the public inputs and archived moments. Equal nominal drift and temperature parameters need not specify equal laboratory-frame moments in a relativistic particle initializer and a nonrelativistic fluid initializer. This identifies controls needed to interpret the comparison; it does not determine their effect on its later nonlinear outcome. [under "1.1 Relation to established work"] Work on heat-flux closure for nonlinear electrostatic two-stream systems [8, 9] further limits a literal reading of broad claims that no fluid closure can capture non-Maxwellian counter-streaming. It does not establish prior use of the specific tagged-population closure in electron-only reconnection. Our analysis therefore addresses explicit initialization, continuity and stability statements, and preserves the source’s application-specific empirical evidence. …and 35 more blocks not shown.