{"id":"fcefad50-82b2-45a5-a857-517480f402f6","arxiv_id":"2411.14255","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"PACIAE 4.0 is a new public release of an event generator that adds PYTHIA 8.3 support, heavy-quark mass corrections, an improved coalescence model, and new hadron channels, but no validation is shown.","lead":"This paper describes PACIAE 4.0, a rewritten version of a Monte Carlo event generator for high-energy particle collisions. It adds a C++ interface to PYTHIA 8.3 and new heavy-quark, coalescence, and hadron-cascade options, but reports no validation results.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"PACIAE 4.0's claim of being a working rewrite with a PYTHIA 8.3 interface rests entirely on an unverified Fortran-C++ bridge and absent physics validation; no build log, test output, or cross-check with PACIAE 3.0 is provided.","rationale":"I agree with the reader's weakest assumption: the operational correctness of the Fortran-C++ interface is the most load-bearing unverified premise. The paper is honest about features still under development (e.g., medium-induced radiation in Section 3.2.2 and PYTHIA 8 decay handling in Section 3.4.2), and those limitations do not threaten the central claim because they are explicitly deferred. The standard physics formulas, such as the LO-pQCD cross sections and the Donnachie-Landshoff parameterization, appear reasonable. The decisive gap is that no evidence is shown that the released code actually compiles, runs, and reproduces known results. For a code-release paper in Computer Physics Communications, this is not a matter of style but of the central claim: a generator that cannot be built or that produces incorrect output would not 'surpass' PACIAE 3.0. The proposed test directly checks both compilation and physical correctness, including the specific risk of double counting or normalization loss when PYTHIA 8 events are fed into the cascade stages. Since the reader's conditional verdict already captures this, I keep the verdict unchanged.","tokens_in":15405,"tokens_out":3979,"duration_ms":37884,"concrete_test":"Clone the linked repository (https://github.com/ArcsaberHep/PACIAE4), follow the documented build with a standard gfortran/g++ toolchain and PYTHIA 8.3, and confirm that the example compiles and runs. Then run the same physics setup (e.g., pp at 13 TeV, and a PbPb configuration) in both C_PY6 and C_PY8 modes, and compare charged-hadron multiplicity, transverse-momentum spectra, and yields against PACIAE 3.0 and against ALICE data; the PYTHIA8 mode should match the PYTHIA6 mode within statistical uncertainties after accounting for expected differences. Also inspect the interface module to verify that PYTHIA's event record is filled and read back correctly for each NN collision in pA/AB runs, with no re-initialization or double counting.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central assertion is that PACIAE 4.0 is a functioning successor to PACIAE 3.0, with a modern Fortran/C++ codebase and a working PYTHIA 8.3 interface. The load-bearing premise is operational correctness: the modern Fortran INTERFACE block described in Section 2.1 must instantiate and delete PYTHIA 8 objects, exchange event records through the C++ API, and feed the PYTHIA-8-generated initial partonic state into the existing parton and hadron cascades without double counting NN collisions or losing normalization. The paper provides no build log, no test output, no example executable, and no comparison with PACIAE 3.0 or with experimental data. The claimed physics improvements (heavy-quark mass corrections in Section 3.2.1, the improved coalescence model in Section 3.3.2, and new heavy-hadron reaction channels in Section 3.4.1) are presented only as formulas or option lists; there is no evidence that they are actually implemented and tested in the released repository. Because the entire value of a code-release paper depends on the code being usable and correct, this unverified interface and missing validation is the most fragile point in the argument.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript announces PACIAE 4.0, a new release of the PACIAE Monte Carlo event generator. The paper describes a rewrite of the former FORTRAN 77 code in modern Fortran mixed with C++, the addition of a C++-based PYTHIA 8.3 interface alongside the retained PYTHIA 6.4 support, and several physics improvements: heavy-quark mass corrections in parton-parton scattering cross sections, an improved coalescence hadronization model, and additional hadron-hadron reaction channels involving heavy hadrons. The paper also sketches the program flow, lists seven simulation frameworks, and states that the code is publicly available on GitHub and Gitee. The central claim is that PACIAE 4.0 is a working successor to PACIAE 3.0 that 'surpasses' it.","tokens_in":15702,"tokens_out":3110,"duration_ms":30683,"significance":"If the software actually builds and runs as described, PACIAE 4.0 would be a valuable public resource for the heavy-ion and high-energy physics community, connecting the modern PYTHIA 8.3 initial-state generator to the PACIAE parton and hadron cascade stages, with heavy-quark corrections and an expanded hadronic channel set. The availability of the source code on public repositories is a real strength, and the paper is written in the established style of the PACIAE series. However, the paper provides no evidence that the new interface and physics implementations are operational: there are no build instructions, no test outputs, no comparisons with PACIAE 3.0 or experimental data, and no sample event output. The value of a code-release paper rests on the code being usable and correct, and that premise is currently unsupported by the manuscript.","major_comments":[{"comment":"The claim that PACIAE 4.0 is 'surpassing the version 3.0' is load-bearing but is not supported anywhere in the manuscript. No comparison with PACIAE 3.0 results, no validation against experimental data, and no benchmark of the new physics modules are presented. I recommend either removing 'surpassing' and replacing it with a neutral description of the upgrade, or adding a section with quantitative evidence (e.g., comparisons of particle yields, pT spectra, or cross-section checks).","section":"Abstract and Summary"},{"comment":"The Fortran-C++ interface is the most fragile technical component: the paper states that the Fortran main program 'instantiates and deletes PYTHIA 8 objects on the heap' via the modern Fortran INTERFACE block, but no build log, compiler/toolchain version, or minimal working example is provided. I ask the authors to include a short 'Installation and smoke test' subsection reporting the tested compilers (e.g., gfortran and g++ versions), the required PYTHIA 8 version, and the output of a simple run (e.g., a few pp events at a fixed energy) to demonstrate that the described interface actually compiles and runs.","section":"Section 2.1, program flow"},{"comment":"The claimed physics improvements — heavy-quark mass corrections in the parton cascade, the improved coalescence model with gluon splitting and quark deexcitation, and the new heavy-hadron reaction channels — are presented only as formulas or reaction lists. There is no evidence that these are implemented in the released code and that they produce reasonable results. For a code paper, at least one demonstrative output is expected, such as a plot of the heavy-quark scattering cross section versus energy, a coalescence-model yield for a known system, or a comparison of a few hadronic observables with PACIAE 3.0 or PYTHIA 8 alone. Without such a check, the reader cannot verify the central claim that the new version works as described.","section":"Sections 3.2.1, 3.3.2, and 3.4.1"}],"minor_comments":[{"comment":"The text says 'the objective-oriented C++'; this should be 'object-oriented C++'.","section":"Section 1, Introduction"},{"comment":"The description of the decoupled input files ('usu.dat', 'pythia6_extra.cfg', 'pythia8_extra.cfg') would be clearer with a short example of each file type or a pointer to the repository for their default contents.","section":"Section 2.3, input files"},{"comment":"Equation (2) defines the nuclear density normalization ρ0, but it is not explicitly stated that this is normalized such that the integral over the Woods-Saxon profile equals A; stating this explicitly would remove ambiguity.","section":"Section 3.1.2, pA and AB collisions"},{"comment":"The two free parameters a and b in the Lund deexcitation function, Eq. (13), are introduced without any discussion of their default values or sensitivity; the paper should at least quote typical values and cite the reference for them.","section":"Section 3.3.2, coalescence model"},{"comment":"Since the code is publicly hosted, the paper should give a specific version tag or commit hash used for the described results, so that readers can reproduce the exact state of the code that is being introduced.","section":"Section 4, Summary"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a standard software-description paper for Computer Physics Communications, and the lack of any validation or build verification is unusual for this venue. The authors may argue that the code repository itself is the artifact and that the paper is only a 'brief introduction'. However, the abstract's 'surpassing version 3.0' claim raises the bar: a quantitative comparison or at least a basic smoke-test output should be mandatory before acceptance. The listed minor issues are easy to fix. I do not see any signs of circular reasoning or unsupported external claims beyond the missing verification of the new code's functionality."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a legitimate code-release paper for a generator people actually use, and the new capabilities are real, but the central claim—PACIAE 4.0 works and supersedes 3.0—is currently taken on faith. There are no validation plots, no comparison with 3.0, no build log or test output, and the Fortran-C++ interface described in Sec. 2.1 is exactly the kind of thing that breaks in practice. The stress-test note lands where the load-bearing fragility is.\n\nWhat's genuinely new: the first PACIAE version with a PYTHIA 8.3 interface, a rewrite from fixed-format FORTRAN 77 to modern Fortran mixed with C++, mass-corrected heavy-quark matrix elements for the parton cascade (Eqs. 9 and 10), an improved coalescence model with Altarelli-Parisi-based gluon splitting and Lund deexcitation, an Angantyr mode, and a set of new heavy-hadron reaction channels. These are incremental but real additions to a code line with a long history. The code is on GitHub and Gitee under GPLv2, which is the right way to release it. The authors are also honest about limits: medium-induced radiation is explicitly flagged as under development.\n\nThe soft spots are about evidence, not physics. The formulas are standard and cited from Combridge, Altarelli-Parisi, PYTHIA, and the additive quark model, and there is no circularity: the self-citations are historical context, and nothing is fitted to data. But the paper never shows one event, one distribution, or one cross-check that PYTHIA 8 initial states feed into the PACIAE cascade without double-counting NN collisions or losing normalization. Key inputs—Debye screening cutoff, Lund a and b, constant cross sections for channels without data—are left unspecified. For a 'brief introduction' one could argue the repository is the validation, but the abstract says 'surpassing' version 3.0, and the text provides no evidence for that. A serious referee should ask for at least a smoke test: compile instructions, one pp and one Pb-Pb run, and a comparison with PACIAE 3.0 under identical settings.\n\nWho it is for: heavy-ion phenomenologists who want a PYTHIA 8 based cascade option, and the existing PACIAE user base. The paper deserves peer review at CPC—software papers are the right venue for this—but it should come back with required validation. If the code runs as described, it is a useful contribution.","headline":"PACIAE 4.0 is a plausible, honest code-release paper with real new capabilities, but the core operational claim rests on an unverified Fortran-C++ bridge and no validation, so it needs revision rather than acceptance as is.","tokens_in":16272,"tokens_out":3383,"would_cite":false,"duration_ms":30107,"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":"PACIAE 4.0 is a rewritten, open-source event generator that joins a modern Fortran core with C++-based PYTHIA 8.3 and adds heavy-quark physics to its parton and hadron cascades.","keywords":["event generator","high-energy collisions","transport cascade model","partonic rescattering","hadronic rescattering","PYTHIA 8.3 interface","heavy-quark cross sections","coalescence hadronization"],"falsifier":"A reader could settle the central claim by compiling the released source, running the same pp collision through both the PYTHIA 6.4 and PYTHIA 8.3 modes with identical PACIAE settings, and checking that final multiplicities and spectra agree within generator differences.","tokens_in":15214,"feed_emoji":"⚛️","tokens_out":11061,"duration_ms":84246,"temperature":0.7,"pith_summary":"This paper introduces PACIAE 4.0, a new generation of the PACIAE Monte Carlo event generator used to simulate high-energy particle and nuclear collisions. The central claim is that the code has been rewritten from fixed-format FORTRAN 77 into free-format modern Fortran mixed with C++, and that it now interfaces with the C++-based PYTHIA 8.3 while keeping the older PYTHIA 6.4 as an option. The authors report three physics upgrades in this version: mass-corrected leading-order cross sections for heavy-quark scattering, an improved coalescence hadronization model with gluon splitting and quark deexcitation, and dozens of new inelastic hadron-hadron channels involving strangeness and charm. If the code performs as described, the heavy-ion physics community gains an open-source generator that can start from PYTHIA 8.3 initial states and then evolve partonic and hadronic cascades with better treatment of heavy quarks.","feed_headline":"PACIAE 4.0 plugs PYTHIA 8.3 into a Fortran-C++ event generator","feed_subtitle":"A rewritten heavy-ion cascade code adds heavy-quark mass corrections and new hadronic channels.","key_machinery":"The load-bearing mechanism is the Fortran-C++ interface module, built on the modern Fortran INTERFACE block, which lets a Fortran main program create and destroy PYTHIA 8.3 objects on the heap and exchange event information with them. This interface is what makes the whole upgrade possible: it is the bridge that lets PACIAE keep its Fortran skeleton while drawing initial states and string fragmentation from the C++ generator. On top of that bridge sit the physics additions that the paper emphasizes: the mass-corrected heavy-quark matrix elements of Eqs. (9) and (10), the gluon-splitting and quark-deexcitation rules of Eqs. (11)-(13) in the coalescence model, and the list of new hadronic reaction channels in Sec. 3.4.1.","core_discovery":"The authors assert that PACIAE 4.0 is the first version of PACIAE to leave FORTRAN 77 and to connect to the C++-based PYTHIA 8.3, using a modern Fortran INTERFACE block to instantiate and delete PYTHIA 8 objects on the heap while a Fortran main program steers the simulation. The release also splits input into a PACIAE control file and separate PYTHIA 6/8 card files, and adds new simulation modes based on PYTHIA 8 and on the Angantyr heavy-ion initial-state model. Physically, PACIAE 4.0 introduces mass-dependent matrix elements for heavy-quark scattering off light quarks and gluons (Qq -> Qq and Qg -> Qg), an improved coalescence model in which gluons split to quark-antiquark pairs sharing the gluon's forward light-cone momentum and energetic quarks deexcite by exciting pairs from the vacuum, and an extended hadronic rescattering stage that includes charmed and strange reaction channels such as J/psi+N, D+pi, and K+Xi. The paper's claim is that these changes make PACIAE 4.0 able to model a wider range of collisions, including heavy-flavor observables, than its predecessor.","pith_inferences":["A natural next step, which the paper itself flags, is to replace the PYTHIA 6 decayer with the PYTHIA 8 decayer; until then, the hadron-decay stage remains the only part of the chain still tied to the older generator.","Because PYTHIA parameters are now in separate card files, systematic studies of generator tuning are easier to automate; one could, for instance, scan over color-reconnection settings in PYTHIA 8 without touching PACIAE's own parameters.","The decoupling of input files also suggests a testable extension: using the same PACIAE cascade settings with PYTHIA 6.4 and PYTHIA 8.3 initial states would isolate how much of any final-state difference comes from the initial-state generator rather than from the cascade physics.","If the Fortran-C++ interface proves stable, the same pattern could let PACIAE call other C++-based generators or modules, though the paper does not claim this."],"forward_implications":["Users can now generate pp, pA, and AB initial states with PYTHIA 8.3 and feed them into PACIAE's parton and hadron cascades, something earlier PACIAE versions could not do because they were tied to PYTHIA 6.4.","The new mass-dependent matrix elements mean charm and bottom quarks are no longer treated as massless in the partonic rescattering stage, which should change heavy-quark energy loss and momentum spectra in simulated events.","The improved coalescence hadronization gives an alternative to string fragmentation that produces hadrons from quark-antiquark pair coalescence with phase-space constraints, allowing the model to address heavy-baryon-to-meson ratios.","The expanded hadronic rescattering channels cover charmonium dissociation (e.g., J/psi + N, J/psi + pi) and open-charm reactions, which are relevant for interpreting quarkonium suppression and heavy-flavor hadron production in nuclear collisions.","The code is released on open-source repositories, so the community can run the PYTHIA 8 mode and compare its output against the PYTHIA 6 mode as a cross-check."],"supporting_citations":[{"why":"PYTHIA 8.3 manual and generator description; the C++-based engine that PACIAE 4.0 newly interfaces with and whose card files it now accepts.","marker":"[23]"},{"why":"PYTHIA 6.4 manual; the Fortran-based engine used by earlier PACIAE versions and retained as an option in 4.0.","marker":"[20]"},{"why":"PACIAE 3.0 paper; the previous version whose cascade frameworks, coalescence model, and rescattering stages PACIAE 4.0 inherits and extends.","marker":"[19]"},{"why":"Combridge's QCD matrix elements for heavy-quark processes; the source of the mass-corrected Qq and Qg cross sections introduced in Sec. 3.2.1.","marker":"[60]"},{"why":"Altarelli-Parisi splitting functions; the kernels used for gluon splitting and quark deexcitation in the improved coalescence model.","marker":"[65]"},{"why":"Compilation of hadronic total cross-section data; the parameterization used for the new inelastic hadron-hadron channels.","marker":"[76]"},{"why":"Angantyr model paper; the heavy-ion initial-state model from PYTHIA 8 used by the new B_ANG and C_ANG simulation modes.","marker":"[35]"}],"fun_headline_variants":["PACIAE 4.0 drops Fortran 77 for PYTHIA 8.3 and C++","Heavy-quark mass terms and new hadronic channels debut in PACIAE 4.0","PACIAE 4.0 rewires cascade physics with PYTHIA 8.3 backend","C++-native PACIAE 4.0 adds charmed and strange hadron channels","PACIAE 4.0 moves beyond PYTHIA 6 with heavy-flavor upgrades"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper's load-bearing premise is that the new Fortran-C++ interface actually compiles and runs with standard toolchains, and that PYTHIA 8.3 events flow into the cascade stages without double-counting nucleon-nucleon collisions or losing normalization.","fun_headline_variants_meta":{"raw":{"variants":["PACIAE 4.0 drops Fortran 77 for PYTHIA 8.3 and C++","Heavy-quark mass terms and new hadronic channels debut in PACIAE 4.0","PACIAE 4.0 rewires cascade physics with PYTHIA 8.3 backend","C++-native PACIAE 4.0 adds charmed and strange hadron channels","PACIAE 4.0 moves beyond PYTHIA 6 with heavy-flavor upgrades"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000965,"raw_usage":{"total_tokens":4150,"prompt_tokens":1029,"completion_tokens":3121,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":645,"completion_tokens_details":{"reasoning_tokens":2994}},"tokens_in":645,"tokens_out":3121,"duration_ms":20036,"temperature":1.0,"reasoning_tokens":2994,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:22:22.649822+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A reader could settle the central claim by compiling the released source, running the same pp collision through both the PYTHIA 6.4 and PYTHIA 8.3 modes with identical PACIAE settings, and checking that final multiplicities and spectra agree within generator differences.","supporting_citations":[{"cited_title":"Baldini, V","cited_arxiv_id":null,"evidence_quote":"Compilation of hadronic total cross-section data; the parameterization used for the new inelastic hadron-hadron channels."}],"review_version":1}