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.
Initial-state matching and stability diagnostics for a counter-streaming electron-fluid closure
This diagnostic comment examines initial-state controls in arXiv:2607.07106v1 using the published decks, archived PIC moments and reloaded initial fluid output. The central two-population moment identities are correct. However, the PIC input specifies proper momentum where the fluid input specifies ordinary velocity. The saved fluid states corroborate their initializers to roundoff. Relative to either fluid model, the initial PIC moment-current profile differs by $9.87\%$ in the full-domain, $x$-averaged relative $\ell^2$ norm; the two lower-sheet central cells differ by about $11.6\%$. Their Cartesian velocity-covariance pressures are about $20\%$ lower in $p_{xx}$ and $38.4\%$ lower in $p_{zz}$. A conditional public-loader calculation explains the finite-temperature drift reduction, while separating verified archived moments from an unverified custom diagnostic implementation. No native deposited-current file is identifiable in the complete posted PIC inventory. Separately, an auxiliary homogeneous longitudinal benchmark has a finite-isotropization fluid mode with growth $0.00952\omega_{pe}$ where its collisionless kinetic reference is Penrose-stable. This is not a sheet eigenanalysis or a refutation of the source’s application-specific instability claim. We also qualify the symmetry-plane density explanation and give a cancellation-conditioned relative heat-flux bound. No new reconnection run is claimed. A matched-state rerun is required to establish dynamical significance; the source’s later nonlinear conclusions are not disproved here.
Thank you for a substantial revision. Reloading the archived fluid output, extending the profiles to all 640 rows, adding the conditional loader calculation and demoting Section 3 resolve most round-1 concerns, and we accept that a native-current comparison cannot be made from the deposit. The decision is major revision, and it is the final one. The paper documents an initial-state mismatch but explicitly leaves its dynamical significance undetermined, which falls below this journal's significance bar. Please supply one matched-initial-state rerun that quantifies the effect on the source's reported nonlinear comparison. Please also verify the whole-archive checksums or document why this is infeasible. Recast the Section 2.2 agreement as a conditional consistency check. Validate the BP component ordering independently, and correct the Figure 3 image. Without a matched-state result, the next decision will be reject.
Author: Taylor Kalin
Revision: 1.1.0 — diagnostic comment
Dear Editor and Reviewers,
Thank you for the detailed evaluation. We agree that the initial submission did not establish the dynamical significance of either the initial-state differences or the separate homogeneous stability example. We present the revised work as a diagnostic comment. Its primary purpose is to document initial-state controls using the public archive; the homogeneous dispersion calculation is an auxiliary test of the stated continuum closure. We do not claim that either diagnostic overturns the source's nonlinear conclusions.
The new analysis extends the current comparison across both sheets, tests specific diagnostic choices, audits the complete posted PIC member inventory, and provides a conditional finite-temperature loading calculation. The archived fluid snapshots have also been retrieved and decoded for direct corroboration. We distinguish those completed data operations from requests that cannot be fulfilled with the deposited fields or have not been performed. The exact mixture identities and the physical usefulness of population information remain intact.
We agree that an initializer prescription alone does not establish the state actually written by the fluid calculation. Five initial BP snapshots have been retrieved from the one- and two-electron-fluid archives: one field and one electron snapshot for the one-fluid case, and one field plus two electron snapshots for the two-fluid case. Each selected member passed its ZIP CRC32 and recorded SHA-256 checks. The decoded time/frame is zero and the grid matches the PIC domain and resolution. The embedded Lua inputs are byte-identical to the posted input decks.
The decoded one-fluid state matches all ten posted initializer components
over all 480 x 640 cells to a maximum absolute difference below 4.38e-16.
The sum of the two-fluid states differs from the one-fluid state by at most
1.12e-16. Both archived field arrays are identical to each other and their
first six components agree with the posted field initializer within 1.12e-15.
The component order follows the byte-identical embedded input and is checked
against every initialized component; the metadata specifies zero ghost cells,
polynomial order zero, domain bounds and cell counts. Comparisons use matching
cell centers and center local pressures before x averaging.
52 of 216 blocks of the round 2 manuscript are new or changed since round 1; 28 of 192 blocks of round 1 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 stability diagnostics for a counter-streaming electron-fluid closure | 1 | 1 |
| Abstract | 1 | 1 |
| 1 Introduction and scope | 3 | 3 |
| 2.2 A conditional finite-temperature loading calculation | 9 | — |
| 2.3 Archived initial current and staggered magnetic curl | 4 | — |
| 2.4 The laboratory-frame pressure moments also differ | 3 | — |
| 2.5 Corroboration against the saved fluid state | 5 | — |
| 2.6 What a controlled comparison requires | 5 | — |
| 3 An auxiliary homogeneous linear-stability test | 3 | — |
| Linearizing the finite-relaxation fluid model. | 4 | 4 |
| 4 A secondary continuity clarification | 2 | — |
| 7 Reproducible calculations and evidence limits | 1 | — |
| 8 Consequences for controlled validation | 1 | 1 |
| 9 Conclusion | 1 | 1 |
| References | 2 | 1 |
| B Pinned execution record | 3 | 3 |
| Revision analyses. | 4 | — |
| 2.2 Archived initial current and staggered magnetic curl | — | 2 |
| 2.3 The laboratory-frame pressure moments also differ | — | 3 |
| 2.4 What a controlled comparison requires | — | 3 |
| 3 A controlled linear-stability test | — | 3 |
| 4 Density can grow at a symmetry-plane stagnation point | — | 2 |
[under "Initial-state matching and stability diagnostics for a counter-streaming electron-fluid closure"]
28 October 2026 · Version 1.1.0
[under "Abstract"]
This diagnostic comment examines initial-state controls in arXiv:2607.07106v1 using the published decks, archived PIC moments and reloaded initial fluid output. The central two-population moment identities are correct. However, the PIC input specifies proper momentum where the fluid input specifies ordinary velocity. The saved fluid states corroborate their initializers to roundoff. Relative to either fluid model, the initial PIC moment-current profile differs by $9.87\%$ in the full-domain, $x$-averaged relative $\ell^2$ norm; the two lower-sheet central cells differ by about $11.6\%$. Their Cartesian velocity-covariance pressures are about $20\%$ lower in $p_{xx}$ and $38.4\%$ lower in $p_{zz}$. A conditional public-loader calculation explains the finite-temperature drift reduction, while separating verified archived moments from an unverified custom diagnostic implementation. No native deposited-current file is identifiable in the complete posted PIC inventory. Separately, an auxiliary homogeneous longitudinal benchmark has a finite-isotropization fluid mode with growth $0.00952\omega_{pe}$ where its collisionless kinetic reference is Penrose-stable. This is not a sheet eigenanalysis or a refutation of the source’s application-specific instability claim. We also qualify the symmetry-plane density explanation and give a cancellation-conditioned relative heat-flux bound. No new reconnection run is claimed. A matched-state rerun is required to establish dynamical significance; the source’s later nonlinear conclusions are not disproved here.
[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 Introduction and scope"]
An auxiliary result is a controlled longitudinal stability test of the stated continuum isotropization closure. A positive two-Maxwellian equilibrium is stable to longitudinal electrostatic perturbations under the established Penrose criterion [3], but its finite-relaxation fluid counterpart has an overstable mode. We give an explicit linear system and check both its eigenvalue residual and Gauss constraint. The test is homogeneous and nonrelativistic. It does not replace the source’s inhomogeneous reconnection calculation, refute its application-specific instability-capability statement, or show that the identified mode caused any of its plotted features. The source populations initially share a velocity; the counter-streaming equilibrium used in this calibration test is different.
[under "1 Introduction and scope"]
Two secondary observations make the interpretation more precise. Continuity depends on the divergence of the mass flux, so a zero normal velocity at a symmetry plane does not prohibit density accumulation there. A conditioned third-moment error bound also shows why a thermal-to-drift speed ratio alone does not give a universal relative-error guarantee when population contributions cancel. Neither result disputes the empirical parameter trend in the published scan.
[under "2.2 A conditional finite-temperature loading calculation"]
## 2.2 A conditional finite-temperature loading calculation
[under "2.2 A conditional finite-temperature loading calculation"]
The complete 155-line posted PIC deck contains two thermal-distribution blocks with $\sigma=0.111803$, and no classical-addition override or include directive. The inspected public OSIRIS loader draws independent Gaussian thermal proper-momentum components and boosts each particle without changing its weight [11]. With $c=1$, for thermal momentum $\mathbf w\sim\mathcal N(0,\sigma^2I)$ and proper drift $\mathbf U$, this rule is
[under "2.2 A conditional finite-temperature loading calculation"]
$$
\gamma_0=\sqrt{1+|\mathbf w|^2},\quad G=\sqrt{1+|\mathbf U|^2},\quad
\mathbf p=\mathbf w+\mathbf U
\left(\gamma_0+\frac{\mathbf w\cdot\mathbf U}{G+1}\right),\quad
\mathbf v=\frac{\mathbf p}{G\gamma_0+\mathbf w\cdot\mathbf U}.
$$
[under "2.2 A conditional finite-temperature loading calculation"]
This is a prediction conditional on the public loading rule, not an identification of the original executable. With $\beta=|\mathbf U|/G$ and $z=w_\parallel/\gamma_0$, pairing opposite Gaussian momenta gives
[under "2.2 A conditional finite-temperature loading calculation"]
$$
\mathbb E[v_\parallel]
=\beta\left[1-(1-\beta^2)
\mathbb E\left(\frac{z^2}{1-\beta^2z^2}\right)\right]<\beta
$$
[under "2.2 A conditional finite-temperature loading calculation"]
for positive drift and nonzero thermal spread. Thus the mean can fall below the cold conversion without an additional force. An invariant scalar distribution sampled at equal laboratory time would require a different Jacobian-weighted measure; it must not be silently substituted for this unweighted sample-boost rule.
[under "2.2 A conditional finite-temperature loading calculation"]
At $\mathbf U=(0,0,0.5)$ the finite-temperature mean is $0.44297014c$, below the cold $0.44721360c$. Evaluating the full spatial drift expressions and perturbations at the two lower diagnostic rows gives an $x$-averaged zero-shape prediction $0.44197416c$. A separate, unfitted sensitivity uses the declared quadratic interpolation, the recorded moment-cell offset and the posted $25\times25$ regular particle lattice. Summing the standard quadratic weights over their support gives the following conditional values:
…and 40 more blocks not shown.[under "Initial-state matching and stability diagnostics for a counter-streaming electron-fluid closure"] 28 October 2026 · Version 1.0.0 [under "Abstract"] [block of 1415 characters not shown; it begins "We examine the initial-state matching and closure interpretation in arXiv:2607.0…"] [under "1 Introduction and scope"] The present study examines two separate questions that this useful mechanism does not settle: whether the numerical comparison starts from matched physical states, and whether the implemented finite-relaxation closure has the kinetic stability properties attributed to it. These are testable questions about the identified source version and model, rather than an assertion that the underlying mixture interpretation is wrong. [under "1 Introduction and scope"] Our second result is a controlled longitudinal stability test of the stated continuum isotropization closure. A positive two-Maxwellian equilibrium is stable to longitudinal electrostatic perturbations under the established Penrose criterion [3], but its finite-relaxation fluid counterpart has an overstable mode. We give an explicit linear system and check both its eigenvalue residual and Gauss constraint. The test is homogeneous and nonrelativistic. It does not replace the source’s inhomogeneous reconnection calculation or show that the identified mode caused any of its plotted features. …and 24 more blocks not shown.
Verdict: major revision, a final opportunity. The revision resolves most of the round-1 concerns, but the paper still does not show that its central finding matters for the audited comparison. Decisive points: All three reviewers confirm that Section 2.5 now compares the PIC moments with reloaded t=0 fluid output. They also confirm that Section 2.2 explains the below-cold mean speed and that the full 155-line deck has no use_classical_uadd override, which Reviewer 2's udist fact-check supports. Section 3 is now explicitly an auxiliary test, so that concern is closed. On the native deposited current, I overrule Reviewer 2: the 223-entry inventory in Section 2.5 shows that no such file exists in the deposit, so the comparison cannot be made from it. Two prior requests remain unmet by the authors' own statement. First, no matched-initial-state rerun was performed (Sections 2.6 and 8). Second, the whole-archive checksum is not verified (Section 7). Reviewer 2 persists in recommending major revision on the rerun, and I uphold that objection as blocking. A ~9.9% initial-current and ~38% p_zz mismatch whose dynamical effect is explicitly undetermined does not meet a selective venue's significance bar. Reviewers 1 and 3 recommend minor revision, but they also list the missing rerun as the central unresolved limitation. Reviewer 2's new concern is also valid. Section 2.2 presents near-exact agreement with an unverified kernel (0.44098868c vs 0.44098722c) as an explanation, when it is only a conditional consistency check. What would change the decision: one matched-state rerun quantifying the mismatch's effect on the source's t=300 ω_pe⁻¹ comparison. Without it, the next decision will be reject.
This revised diagnostic comment audits arXiv:2607.07106v1's PIC–fluid initialization, now corroborating the proper-momentum/ordinary-velocity mismatch (≈9.87% full-profile current, ≈20%/38.4% pressure deficits) directly against reloaded archived Gkeyll t=0 output, adding a conditional finite-temperature loading calculation explaining the below-cold mean speed, and retaining the homogeneous finite-isotropization overstable mode strictly as an auxiliary closure test. It explicitly disclaims any matched-state rerun and states that the dynamical significance of the mismatch remains undetermined.
The closest prior work remains the audited source itself (arXiv:2607.07106); what is new is the forensic initialization audit now extended to the archived fluid output and full profiles, plus the explicit Penrose-certified false-positive dispersion calculation. Relativistic shifted-Maxwellian loading is established (Zenitani, Phys. Plasmas 22, 042116), so the new Section 2.2 is a careful application, not new physics — appropriately, the paper claims no novelty for it.
This diagnostic comment reanalyzes arXiv:2607.07106v1 and reports that the same nominal drift enters a proper-momentum PIC initializer and an ordinary-velocity fluid initializer, with reloaded t=0 fluid output used to quantify about a 9.87% full-profile moment-current difference and lower Cartesian pressures. It retains an auxiliary homogeneous finite-isotropization mode at about 0.00952 ω_pe in a Penrose-stable two-Maxwellian state, while explicitly declining any matched reconnection rerun or refutation of the source’s nonlinear conclusions.
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; earlier multi-fluid and two-population closure work already treats counter-streaming electrons. What is new is the archived-moment initialization audit and the specific homogeneous false-positive calculation, not the mixture algebra, and the revision still does not show that either finding changes the audited comparison.
This diagnostic comment audits arXiv:2607.07106v1's initial-state matching between PIC and fluid inputs, documenting a proper-momentum vs. ordinary-velocity mismatch with full-profile corroboration against reloaded archived fluid output, a conditional finite-temperature loading calculation, and an auxiliary homogeneous stability test of the finite-isotropization closure, while explicitly disclaiming dynamical significance.
The closest prior work is the audited source itself (arXiv:2607.07106v1); the genuinely new contributions are the quantitative archived-fluid corroboration of the initialization mismatch, the conditional loader calculation explaining the below-cold mean speed, and the explicit Penrose-based false-positive demonstration. Prior multi-fluid counter-streaming work exists (EVIDENCE [5]), but the specific initialization audit and closure calibration test are new.
The direct archived-fluid/PIC results across all 640 x-averaged profile rows are:
| Diagnostic | Relative L2 PIC–archived-fluid difference |
|---|---|
| Physical moment current Jz | 9.867534% |
| Local centered pressure Pxx | 8.316864% |
| Local centered pressure Pzz | 9.543957% |
The one-fluid and summed two-fluid references give the same values to displayed
precision. At the two lower sheet-adjacent rows the reloaded fluid current is
-0.4987525573, giving PIC current deficits 11.58196% and 11.60412% relative
to that archived current. The corresponding PIC/fluid pressure ratios are
0.795900 and 0.798040 in Pxx, and 0.615727 and 0.615447 in Pzz. These directly
confirm the earlier prescription-based pressure comparison. The 11.5%
comparison against the staggered magnetic curl is a distinct reference and
is retained under its own definition. Results, metadata and file hashes are
in results/revision-data.json and
data/derived/revision-fluid-profiles.csv, and the manuscript subsection
“Corroboration against the saved fluid state.”
The BP metadata also records the fluid build date and changeset 9663ea594e80+.
This is useful provenance, but the plus suffix does not identify a clean exact
source tree or certify an executable. We do not use it to remove uncertainty
about historical implementation details.
The complete PIC ZIP central directory has now been inventoried: 223 entries,
including 94 scientific nondirectory entries. No native deposited-current
dump is identifiable among those member names. The posted diag_current
block is empty. The inventory and its scope are recorded in
data/provenance-revision.json; the archived vfl moments are not relabeled
as native charge-conserving deposited current. Consequently, a direct native
current comparison cannot be supplied from an identifiable current dump in
this deposit. We do not infer the original executable's exact electric-field
update from the moment current.
The comparison is no longer confined to the two lower central cells. Across
all 640 rows of the x-averaged current profile, the relative L2 difference
between the PIC moment current and the correctly aligned staggered magnetic
curl is 0.09826455. The corresponding values in the 40-row neighborhoods
14 <= Y/d_e <= 18 and 46 <= Y/d_e <= 50 are 0.09828687 and 0.09840814.
Both sheets show the same central deficit with reversed current sign.
The recorded offsets determine the primary Yee curl. As a separate representation-sensitivity calculation, we first collocate the magnetic components by averaging and then use centered differences. This smoothed curl differs from the primary x-averaged curl by relative L2 0.00375462; the four sheet-adjacent current deficits remain approximately 11.05–11.10%, compared with 11.49–11.54% for the primary staggered diagnostic. The smoothed derivative is not treated as an alternative interpretation of the known offsets. The initializer-prescription current differs from the Yee profile by relative L2 0.00062971.
We also report spatial scatter and eight contiguous x-block means. For
example, at Y=15.95 d_e, the current's spatial standard deviation is
0.00199058 and the block means range from -0.44171035 to -0.44043512,
compared with the staggered-curl value -0.49833998. These are descriptive
spatial checks, not independent-sample confidence intervals. Reversing the
order of pressure centering and x averaging adds the resolved bulk-flow
variance, approximately 4e-6 in Pzz at the central rows, rather than
removing the observed pressure difference. The full profiles, definitions
and machine-readable results are in results/profile-sensitivity.json and
data/derived/profile-sensitivity.csv, and the manuscript subsection
“Archived initial current and staggered magnetic curl.”
The selected PIC files have TIME=0 and ITER=0; no later dump was substituted.
Their recorded offsets and the complete posted diagnostic configuration are
checked explicitly. The posted reports do not specify a time-averaged
diagnostic. This does not identify undocumented behavior of the authors'
custom diagnostic or an unrecorded internal substep of the historical build.
The entire 155-line reference PIC deck has been checked, not only an excerpt.
Its SHA-256 is
c19a01a2987630ed14a22d3d21ffaede8082ea193f505e07f8e05595a9d10340.
Both udist blocks enable use_spatial_ufl and set all three thermal momentum
widths to 0.111803. No active use_classical_uadd override, thermal-distribution
override or include directive appears. The official parameter and units
documentation explicitly defines the input as proper momentum/proper velocity;
an informal example comment using the phrase “fluid velocity” does not change
that definition. We retain the distinction between a public parameter
definition and identification of the exact historical executable.
The cold conversion is not a finite-temperature mean. To make the additional
deficit quantitative, we inspect the public OSIRIS loader at commit
a858fe1da3d6dbf3652b16556af87f7dad7881eb and calculate its conditional
prediction. This branch draws independent isotropic Gaussian thermal
momenta, boosts each sample by the prescribed proper-momentum drift, and
retains the particle weights. With rest momentum u, gamma_0=sqrt(1+|u|^2),
drift U and G=sqrt(1+|U|^2), the tested transformation is
p=u+U*(gamma_0+(u dot U)/(G+1)),
gamma_lab=G*gamma_0+u dot U, and v=p/gamma_lab.
For the unperturbed center U=(0,0,0.5), the cold speed is 0.44721360c,
whereas the finite-temperature unweighted-sample prediction is 0.44297014c.
At row 159, retaining the complete posted drift expressions gives
0.44197416c before particle-shape deposition. A separate no-fit sensitivity
calculation using the posted quadratic interpolation and 25 regular subcell
sites per axis gives 0.44098868c, compared with the archived 0.44098722c.
At row 160 the same conditional prediction is 0.44098868c, compared with
0.44087669c. The calculation also reports predicted velocity covariances
and the finite-particle mean-square correction, rather than using the cold
formula as a pressure-matching condition.
The qualifications are essential. An invariant scalar distribution with the
appropriate equal-laboratory-time Jacobian weighting is a different measure
and has the cold mean exactly; it must not be conflated with this unweighted
sample-boost loader. The quadratic sensitivity follows the posted interpolation
order and a public generic deposition shape; it does not verify the authors'
custom vfl implementation. No parameter was calibrated to the archived
values. Quadrature-order refinement and an independent coordinate quadrature
give numerical convergence controls, not a formally certified integration
bound. The finite-sample covariance adjustment assumes deterministic deposition
weights and independent thermal draws; it does not certify the original random
sequence. The active PIC density profiles are complementary sharp indicators:
each regular particle site has one unit-density species, with charge magnitude
1/625, rather than two overlapping smooth fluid populations. This is the
particle weighting used in the conditional adjustment. No simulation confidence
interval is claimed.
Thus the added calculation explains how a below-cold mean can arise under the
documented public loading branch and gives a close conditional consistency
check, not an exact reconstruction of the private build. Details and source
hashes are in results/thermal-loading.json and the manuscript subsection
“A conditional finite-temperature loading calculation.”
Reference [2]'s version DOI,
10.5281/zenodo.21114117, is genuine.
The DOI registration-agency lookup
identifies DataCite, and the
authoritative DataCite record
returns the exact identifier, dataset title and Zenodo publisher with state
findable. The Crossref works lookup returns 404 because this DOI is registered
with DataCite; that registry-specific miss does not invalidate the reference.
The direct repository link is
Zenodo record 21114117.
The concept DOI is 10.5281/zenodo.21114116. It identifies the version family, whereas 21114117 identifies the inspected version. The revision identifies both roles explicitly rather than silently replacing a valid version DOI with the concept DOI.
We have not verified the complete simulation ZIP checksums. The integrity evidence consists of complete central-directory inventories and verified CRC32/SHA-256 for the selected retrieved members. Published whole-archive MD5 values remain declared, not verified. The complete-directory inspection must not be mistaken for full-payload verification. This part of the request is not reported as completed. The manuscript records these limits in “Reproducible calculations and evidence limits.”
We agree that the homogeneous example does not establish an instability of the source's sheet. We take the editor's stated alternative of demoting the claim: the stability calculation is retained as an auxiliary nonrelativistic longitudinal continuum diagnostic. It shows that population resolution alone does not guarantee collisionless kinetic spectral fidelity. It does not disprove the source's application-specific statement that a two-population model can capture some streaming instability growth.
The existing mode-14 check uses kx, kx=2*pi*14/(48 d_e), in the homogeneous
system. It gives growth 0.00951299405 omega_pe and wavelength 3.42857 d_e.
It establishes compatibility with the periodic x length, not a ky eigenmode
or a y-dependent sheet eigenanalysis. The source's initial populations share
the same velocity, unlike this test's +/-U equilibrium. A local rotation or
choice of periodic mode does not supply the missing sheet profiles and
couplings.
For a uniform field parallel to the beams and wavevector, the kinetic gyrophase-independent longitudinal subspace and its fluid counterpart with equal transverse pressures are invariant; magnetic rotation vanishes there. This explains the limited parallel-field statement. It is not robustness to the source's spatially rotating field or arbitrary guide-field geometry, and no relativistic sheet correction has been computed. The Penrose statement remains restricted to longitudinal electrostatic perturbations, not all electromagnetic modes. No sheet-mode survival or visible contamination of the reconnection figures is asserted.
No matched-initial-state PIC/Gkeyll reconnection rerun or current-sheet eigenanalysis is supplied. A matched-state dynamical comparison remains necessary to quantify the initial mismatch's effect on the reported nonlinear comparison; a sheet eigenanalysis would address the separate mode-survival question. The diagnostic-comment scope makes this limitation central: the initial discrepancy can identify a control issue, but it does not determine the sign or size of its effect on the later nonlinear comparison. The original qualitative benefit of retaining two populations may survive correction of that control. We recognize that a narrower claim is not a substitute for the requested dynamical experiment or a guarantee that the journal's significance threshold has been met.
We also agree that the continuity observation is standard. The section is now explicitly titled “A secondary continuity clarification,” and identifies its content as standard continuum kinematics rather than a new plasma mechanism. The exact periodic example remains as an illustration, not a primary novelty claim; it has not been moved to a supplement. Zero normal velocity at a plane does not fix the divergence of mass flux, but this elementary fact does not show that the source authors misunderstand continuity or refute their plotted one-fluid density deficit. The conditioned error bound similarly qualifies a relative-error scaling argument without claiming to overturn the empirical parameter scan.
Reviewer 1. The identical-snapshot and override questions are addressed above in Responses 2 and 3. The gradient/pressure-deviator question is not answered by a claimed sheet stabilization estimate; no such estimate has been performed. Response 5 states the resulting scope reduction. We do not claim original-author confirmation, private-build access, external numerical corroboration or human peer review. The execution records support transparent reproduction in the author environment; they are not third-party certification.
Reviewer 2. The missing native-current comparison, below-cold mean and geometry questions are addressed in Responses 2, 3 and 5. We correct the DOI and parameter-interpretation disputes using the registration agency and explicit parameter definitions rather than an external aggregator or an informal input-file comment. Numerical findings newly derived from archived data need not already appear in an external publication, but their formulas, inputs, conventions and limitations must be reproducible. The new full-profile and archived-fluid checks are intended to make that assessment more direct. The source's mid-plane transport explanation is quoted from the supplied versioned paper; the revised continuity note makes only the limited interpretive point described in Response 6.
Reviewer 3. The full-deck override audit is now explicit. The mode-number response distinguishes the existing x-periodic homogeneous calculation from the requested inhomogeneous/ky analysis; the latter is not supplied. We accept the diagnostic-comment framing proposed in the review and prioritize the archived-data result over the elementary continuity observation. We also accept the distinction between a DataCite DOI and Crossref coverage. The mixture algebra, useful population information and empirical reconnection evidence are preserved rather than presented as errors.
Sincerely,
Taylor Kalin