Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pythia”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Azimuthal correlation anisotropies in p + p collisions simulated using Pythia

Stimulated by a keen interest in possible collective behavior in high-energy proton-proton and proton-nucleus collisions, we study two-particle angular correlations in pseudorapidity and azimuthal differences in simulated p + p interactions using the Pythia 8 event generator. Multi-parton interactions and color connection are included in these simulations, which have been perceived to produce collectivity in final-state particles. Meanwhile, contributions from genuine few-body nonflow correlations, not of collective flow behavior, are known to be severe in these small-system collisions. We present our Pythia correlation studies pedagogically and report azimuthal harmonic anisotropies analyzed using several methods. We observe anisotropies in these Pythia simulated events qualitatively and semi-quantitatively, similar to experimental data. Furthermore, our findings highlight the delicate nature of azimuthal anisotropies in small-system collisions and provide a benchmark that can aid in improving data analysis and interpreting experimental measurements in small-system collisions.

Pythia↗
NASA NTRSDate not supplied

PlanetarY Telemetric Helicopter for Investigation and Analysis (PYTHIA): A Rotorcraft for Martian Lava Tube Exploration

The PlanetarY Telemetric Helicopter for Investigation and Analysis (PYTHIA) project presents an early feasibility study into future Mars rotorcraft. With the success of Ingenuity and the current development of the Mars Sample Retrieval Helicopters, there is motivation to explore additional vehicle concepts for Mars exploration. This work presents an early conceptual design of a lava tube exploring quadrotor. The nominal mission for the PYTHIA quadrotor includes a two-phase in-depth exploration of one or mutliple lava tubes in the area of interest, Arsia Mons. A sizing analysis was completed using the NASA Design and Analysis of Rotorcraft (NDARC) tool. An eight-bladed quadrotor was selected based on the NDARC sizing sweeps. With preliminary vehicle and blade sizing completed, a flow visualization study was conducted. Two simulations were completed in Rotorcraft Computational Fluid Dynamics (RotCFD): one of the quadrotor alone and another in a lava tube. This paper introduces a baseline reference mission, preliminary vehicle design, and initial sizing analysis of a rotorcraft to explore Martian lava tubes.

PlanetarY↗
DOE OSTI2025

Bayesian optimization of PYTHIA 8 tunes

A new tune (set of model parameters) is found for the six most important parameters of the PYTHIA 8 final state parton shower and hadronization model using Bayesian optimization. The tune fits the Large Electron-Positron collider (LEPI) data from ALEPH better than the default tune in PYTHIA 8. To the best of our knowledge, we present the most comprehensive application of Bayesian optimization to the tuning of a parton shower and hadronization model using the LEPI data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

MadGraph5_aMC@NLO + Pythia 8.3 event samples for cell resampling studies (13 TeV, Z+jets and top quark pair production)

This record contains two simulated proton–proton collision samples at a center-of-mass energy of √s = 13 TeV, each of 100,000 events, stored as complete truth-level particle records. The samples were produced for a study of optimal-transport-based cell resampling of negatively weighted Monte Carlo events, and are released so that the study can be reproduced and to support further work on negative-weight reduction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Strangeness enhancement at its extremes: multiple (multi-)strange hadron production in pp collisions at \(\sqrt{s}=5.02\) TeV

The probability to observe a specific number of strange and multi-strange hadrons (nS), denoted as P(nS), is measured by ALICE at midrapidity (|y| < 0.5) in $$\sqrt{s}=5.02$$ TeV proton-proton (pp) collisions, dividing events into several multiplicity-density classes. Exploiting, for the first time, a technique based on counting the number of strange-particle candidates event-by-event, this measurement allows one to extend the study of strangeness production beyond the mean of the distribution. This constitutes a new test bench for production mechanisms, probing events with a large imbalance between strange and non-strange content. The analysis of a large-statistics data sample makes it possible to extract P(nS) up to a maximum nS of 7 for $${\text{K}}_{\text{S}}^{0}$$, 5 for Λ and $$\overline{\Lambda }$$, 4 for Ξ− and $${\overline{\Xi } }^{+}$$, and 2 for Ω− and $${\overline{\Omega } }^{+}$$. From this, the probability of producing strange hadron multiplets per event is calculated, thereby enabling the extension of the study of strangeness enhancement to extreme situations where several strange quarks hadronize in a single event at midrapidity. Moreover, comparing hadron combinations with different u and d quark compositions and equal overall s quark content, the contribution to the enhancement pattern coming from non-strangeness related mechanisms is isolated. The results are compared with state-of-the-art phenomenological models implemented in commonly used Monte Carlo event generators, including PYTHIA 8 Monash 2013, PYTHIA 8 with QCD-based Color Reconnection and Rope Hadronization (QCD-CR + Ropes), and EPOS LHC, which incorporates both partonic interactions and hydrodynamic evolution. These comparisons show that the new approach dramatically enhances the sensitivity to the different underlying physics mechanisms modeled by each generator.

Abualrob, I J↗
DOE OSTI2025

Multiplicity dependent 𝐽/𝜓 and 𝜓⁡(2⁢𝑆) production at forward and backward rapidity in 𝑝 + 𝑝 collisions at $\sqrt{𝑠}$ = 200 GeV

Recent measurements of 𝐽/𝜓 production as a function of event charged-particle multiplicity at the collision energies of both the Large Hadron Collider (LHC) and the Relativistic Heavy Ion Collider (RHIC) show enhanced 𝐽/𝜓 production yields with increasing multiplicity. One potential explanation for this type of dependence is multiparton interactions (MPI). We present the first study of potential autocorrelations at RHIC energies and forward and backward rapidity of self-normalized 𝐽/𝜓 yields and 𝜓⁡(2⁢𝑆) to 𝐽/𝜓 ratio, as a function of self-normalized multiplicity in 𝑝 + 𝑝 collisions. In addition, detailed pythia studies tuned to RHIC energies were performed to investigate the MPI impacts. We find that the PHENIX data at RHIC are consistent with recent LHC measurements and can only be described by pythia calculations that include MPI effects. The forward and backward 𝜓⁡(2⁢𝑆) to 𝐽/𝜓 ratio is found to be less dependent on the charged-particle multiplicity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Understanding baryon stopping at the BNL Relativistic Heavy Ion Collider top energies

The nucleon exhibits a rich internal structure governed by quantum chromodynamics (QCD), where its electric charge arises from valence quarks, while its spin and mass emerge from complex interactions among valence quarks, sea (anti)quarks, and gluons. At the advent of QCD, an alternative hypothesis emerged suggesting, at high energies, the transport of a nucleon's baryon number could be traced by a nonperturbative configuration of gluon fields connecting its three valence quarks, forming a 𝑌-shaped topology known as the gluon junction. Recent measurements by the STAR experiment are compatible with this scenario. In light of these measurements, this study aims to explore the mechanisms of baryon transport in high-energy nuclear collisions using the pythia-8 framework, which incorporates a state-of-the-art hadronization model with advanced color flow (CF) and color reconnection (CR) mechanisms that mimic signatures of a baryon junction. Within this model setup, we investigate (i) the rapidity slope of the net-baryon distributions in photon-included processes (𝛾 + 𝑝) and (ii) baryon over charge transport in the isobaric (Ru + Ru and Zr + Zr) collisions. Our study highlights the importance of the CF and CR mechanisms in pythia-8, which play a crucial role in baryon transport. The results show that the CF and CR schemes significantly affect the isobaric baryon-to-charge ratio, leading to different predictions for baryon stopping and underscoring the need to account for CF and CR effects in comparisons with experimental measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Particle production as a function of charged-particle flattenicity in 𝑝⁢𝑝 collisions at $\sqrt{s}$ = 13 TeV

This paper reports the first measurement of the transverse momentum (𝑝 T ) spectra of primary charged pions, kaons, (anti)protons, and unidentified particles as a function of the charged-particle flattenicity in pp collisions at $\sqrt{s}$ =13 TeV. Flattenicity is a novel event shape observable that is measured in the pseudorapidity intervals covered by the V0 detector, 2.8 < 𝜂 < 5.1 and −3.7 < 𝜂 < −1.7. According to QCD-inspired phenomenological models, it shows sensitivity to multiparton interactions and is less affected by biases toward larger 𝑝 T due to local multiplicity fluctuations in the V0 acceptance than multiplicity. The analysis is performed in minimum-bias (MB) as well as in high-multiplicity events up to 𝑝 T = 20 GeV/⁢𝑐. The event selection requires at least one charged particle produced in the pseudorapidity interval |𝜂| < 1. The measured 𝑝 T distributions, average 𝑝 T , kaon-to-pion and proton-to-pion particle ratios, presented in this paper, are compared to model calculations using pythia 8 based on color strings and EPOS LHC. The modification of the 𝑝 T -spectral shapes in low-flattenicity events that have large event activity with respect to those measured in MB events develops a pronounced peak at intermediate 𝑝 T (2 < 𝑝 T < 8 GeV/𝑐), and approaches the vicinity of unity at higher 𝑝 T . The results are qualitatively described by pythia, and they show different behavior than those measured as a function of charged-particle multiplicity based on the V0M estimator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

$\overline{\Sigma }^{\pm }$ production in $\text {pp}$ and $\text {p}{-}\text{Pb}$ collisions at $\sqrt{s_{\textrm{NN}}} = 5.02$ TeV with ALICE

The transverse momentum spectra and integrated yields of anti-$Σ$ hyperons ($\overline{\Sigma}^{\pm}$) have been measured in and collisions at $\sqrt{s_{\textrm{NN}}} = 5.02$ TeV with the ALICE experiment. Measurements are performed via the newly accessed decay channel $\overline{\Sigma}^{\pm}$ → $\bar{\textrm{n}}π^±$. A new method of antineutron reconstruction with the PHOS electromagnetic spectrometer is developed and applied to this analysis. The p T spectra of $\overline{\Sigma}^{\pm}$ are measured in the range 0.5 < p T < 3 GeV/c and compared to predictions of the PYTHIA 8, DPMJET, PHOJET, EPOS LHC and EPOS4 models. The EPOS LHC and EPOS4 models provide the best descriptions of the measured spectra both in pp and p-Pb collisions, while models which do not account for multiparton interactions provide a considerably worse description at high p T . The total yields of $\overline{\Sigma }^{\pm }$ in both pp and p-Pb collisions are compared to predictions of the Thermal-FIST model and dynamical models PYTHIA 8, DPMJET, PHOJET, EPOS LHC and EPOS4. All models reproduce the total yields in both colliding systems within uncertainties. The nuclear modification factors R pPb for both $\overline{\Sigma}^{+}$ and $\overline{\Sigma}^{-}$ are evaluated and compared to those of protons, $Λ$ and $Ξ$ hyperons, and predictions of EPOS LHC and EPOS4 models. No deviations of R pPb for $\overline{\Sigma}^{\pm}$ from the model predictions or measurements for other hadrons are found within uncertainties.

Abualrob, I. J. [University of Houston] (ORCID:000↗
DOE OSTI2026

Pythia8 Quark and Gluon Jets (float32)

A float32 (single-precision) version of the quark and gluon jet dataset originally published by Komiske, Metodiev, and Thaler (Zenodo record 3164691). Only the 20-file subset without charm and bottom quark jets is included here. All simulation parameters and jet selection criteria are identical to the original: Pythia 8.226, √s = 14 TeV Quarks from WeakBosonAndParton:qg2gmZq, gluons from WeakBosonAndParton:qqbar2gmZg with the Z decaying to neutrinos FastJet 3.3.0, anti-k_t jets with R = 0.4 p_T^jet ∈ [500, 550] GeV, |y^jet| < 1.7 There are 20 files, each in compressed NumPy format (QG_jets_fp32_0.npz through QG_jets_fp32_19.npz). Each file contains two arrays: X: (100000, M, 4) — 50k quark and 50k gluon jets, randomly sorted, padded to max multiplicity M, with particle features (pt, rapidity, azimuthal angle, pdgid) stored as float32 y: (100000,) — jet labels, gluon = 0, quark = 1 The original dataset stores X in float64. Here X has been cast to float32, approximately halving file size. The y labels are unchanged. If you use this dataset, please cite the original Zenodo record and its associated paper: Komiske, Metodiev, Thaler, Energy Flow Networks: Deep Sets for Particle Jets, JHEP 01 (2019) 121, arXiv:1810.05165

energyflow↗
DOE OSTI2026

Pythia8 Quark and Gluon Jets (float16)

A float16 (half-precision) version of the quark and gluon jet dataset originally published by Komiske, Metodiev, and Thaler (Zenodo record 3164691). Only the 20-file subset without charm and bottom quark jets is included here. All simulation parameters and jet selection criteria are identical to the original: Pythia 8.226, √s = 14 TeV Quarks from WeakBosonAndParton:qg2gmZq, gluons from WeakBosonAndParton:qqbar2gmZg with the Z decaying to neutrinos FastJet 3.3.0, anti-k_t jets with R = 0.4 p_T^jet ∈ [500, 550] GeV, |y^jet| < 1.7 There are 20 files, each in compressed NumPy format (QG_jets_fp32_0.npz through QG_jets_fp32_19.npz). Each file contains two arrays: X: (100000, M, 4) — 50k quark and 50k gluon jets, randomly sorted, padded to max multiplicity M, with particle features (pt, rapidity, azimuthal angle, pdgid) stored as float32 y: (100000,) — jet labels, gluon = 0, quark = 1 The original dataset stores X in float64. Here X has been cast to float16, approximately halving file size. The y labels are unchanged. Users should be aware that float16 has limited dynamic range and precision. If you use this dataset, please cite the original Zenodo record and its associated paper: Komiske, Metodiev, Thaler, Energy Flow Networks: Deep Sets for Particle Jets, JHEP 01 (2019) 121, arXiv:1810.05165

energyflow↗
DOE OSTI2025

Heavy flavour jet substructure

We present a comprehensive study of energy correlation functions and jet angularities for heavy-flavour QCD jets. In particular, we discuss the possibility of using these observables to expose the dead cone effect, i.e. the suppression of collinear QCD radiation around massive quarks, and to investigate the sensitivity of different observable definitions to the presence of quark masses. Our calculations are presented as all-order resummed predictions at next-to-leading-logarithmic accuracy, matched to (partial) fixed-order results to obtain a better description of the transition around the dead cone threshold. We also compare our analytic results with Pythia, Herwig and Sherpa Monte Carlo predictions to estimate the impact of non-perturbative contributions such as hadronisation, underlying events and B-hadron decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

First observation of strange baryon enhancement with effective energy in pp collisions at the LHC

The production of (multi-)strange hadrons is measured at midrapidity in proton-proton collisions at $\sqrt{s}$ = 13 TeV as a function of the local charged-particle multiplicity in the pseudorapidity interval |η| < 0.5 and of the very-forward energy measured by the ALICE Zero-Degree Calorimeters. The latter provides information on the effective energy, i.e. the energy available for particle production in the collision once subtracted from the centre-of-mass energy. The yields of K$^{0}_{s}$, $Λ$ + $\overline{Λ}$, and $Ξ$ – + $\overline{Ξ}$ + per charged-particle increase with the effective energy. In addition, this work exploits a multi-differential approach to decouple the roles of local multiplicity and effective energy in such an enhancement. The results presented in this article provide new insights into the interplay between global properties of the collision, such as the initial available energy in the event, and the locally produced final hadronic state, connected to the charged-particle multiplicity at midrapidity. Notably, a strong increase of strange baryon production with effective energy is observed for fixed charged-particle multiplicity at midrapidity. These results are discussed within the context of existing phenomenological models of hadronisation implemented in different tunes of the PYTHIA 8 event generator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multiplicity-dependent jet modification from di-hadron correlations in pp collisions at $ \sqrt{s} $ = 13 TeV

Short-range correlations between charged particles are studied via two-particle angular correlations in pp collisions at $\sqrt{s}$ = 13 TeV. The correlation functions are measured as a function of the relative azimuthal angle ∆φ and the pseudorapidity separation ∆η for pairs of primary charged particles within the pseudorapidity interval |η| < 0.9 and the transverse-momentum range 1 < p T < 8 GeV/c. Near-side (|∆φ| < 1.3) peak widths are extracted from a generalised Gaussian fitted over the correlations in full pseudorapidity separation (|∆η| < 1.8), while the per-trigger associated near-side yields are extracted for the short-range correlations (|∆η| < 1.3). Both are evaluated as a function of charged-particle multiplicity obtained by two different event activity estimators. The width of the near-side peak decreases with increasing multiplicity, and this trend is reproduced qualitatively by the Monte Carlo event generators PYTHIA 8, AMPT, and EPOS. However, the models overestimate the width in the low transverse-momentum region (p T < 3 GeV/c). The per-trigger associated near-side yield increases with increasing multiplicity. Although this trend is also captured qualitatively by the considered event generators, the yield is mostly overestimated by the models in the considered kinematic range. The measurement of the shape and yield of the short-range correlation peak can help us understand the interplay between jet fragmentation and event activity, quantify the narrowing trend of the near-side peak as a function of transverse momentum and multiplicity selections in pp collisions, and search for final-state jet modification in small collision systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

J/ ψ -hadron correlations at midrapidity in pp collisions at $\sqrt{s}$ = 13 TeV

We report on the measurement of inclusive, non-prompt, and prompt J/ψ-hadron correlations by the ALICE Collaboration at the CERN Large Hadron Collider in pp collisions at a center-of-mass energy of 13 TeV. The correlations are studied at midrapidity (|y| < 0.9) in the transverse momentum ranges p T < 40 GeV/c for the J/ψ and 0.15 < p T < 10 GeV/c and |η| < 0.9 for the associated hadrons. The measurement is based on minimum bias and high multiplicity data samples corresponding to integrated luminosities of L int = 34 nb −1 and L int = 6.9 pb −1 , respectively. In addition, two more data samples are employed, requiring, on top of the minimum bias condition, a threshold on the tower energy of E = 4 and 9 GeV in the ALICE electromagnetic calorimeters, which correspond to integrated luminosities of L int = 0.9 pb −1 and L int = 8.4 pb −1 , respectively. The azimuthally integrated near and away side yields of associated charged hadrons per J/ψ trigger are presented as a function of the J/ψ and associated hadron transverse momentum. The measurements are discussed in comparison to PYTHIA calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Robust calibration of non-perturbative models with history matching

We apply, for the first time, Bayes Linear Emulation and History Matching to the calibration of non-perturbative models in Monte Carlo event generators. In contrast to the usual approach of "Monte Carlo tuning", History Matching does not result in best-fit plus ellipsoidal parameter uncertainty estimates but instead identifies all parameter space regions that are consistent with data. This approach leads to a systematic and robust quantification of parametric uncertainties in the models, especially in those challenging cases where different, possibly disjoint, regions of parameter space deliver similar results, which are usually not properly treated with current methodology. We highlight the power of this method with the hadronisation models available through Sherpa : the built-in cluster fragmentation Ahadic and string fragmentation through an interface to Pythia .

Iskauskas, Andrew [Durham U. (main)] (ORCID:000000↗
DOE OSTI2024

The Ghent Hybrid model in NuWro: a new neutrino single-pion production model in the GeV regime

Neutrino-induced single-pion production constitutes an essential interaction channel in modern neutrino oscillation experiments, with its products building up a significant fraction of the observable hadronic final states. Frameworks of oscillation analyses strongly rely on Monte Carlo neutrino event generators, which provide theoretical predictions of neutrino interactions on nuclear targets. Thus, it is crucial to integrate state-of-the-art single-pion production models with Monte Carlo simulations to prepare for the upcoming systematics-dominated landscape of neutrino measurements. In this work, we present the implementation of the Ghent Hybrid model for neutrino-induced single-pion production in the NuWro Monte Carlo event generator. The interaction dynamics includes coherently-added contributions from nucleon resonances and a non-resonant background, merged into the pythia branching predictions in the deep-inelastic regime, as instrumented by NuWro. This neutrino-nucleon interaction model is fully incorporated into the nuclear framework of the generator, allowing it to account for the influence of both initial- and final-state nuclear medium effects. We compare the predictions of this integrated implementation with recent pion production data from accelerator-based neutrino experiments. The results of the novel model show improved agreement of the generator predictions with the data and point to the significance of the refined treatment of the description of pion-production processes beyond the ∆ region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

TPCpp-10M: Simulated proton-proton collisions in a time projection chamber for AI foundation models

Scientific foundation models hold great promise for advancing nuclear and particle physics by improving analysis precision and accelerating discovery. Yet, progress in this field is often limited by the lack of openly available large scale datasets, as well as standardized evaluation tasks and metrics. Furthermore, the specialized knowledge and software typically required to process particle physics data pose significant barriers to interdisciplinary collaboration with the broader machine learning community. This work introduces a large, openly accessible dataset of 10 million simulated proton-proton collisions, designed to support self-supervised training of foundation models. To facilitate ease of use, the dataset is provided in a common NumPy format. In addition, it includes 70,000 labeled examples spanning three well defined downstream tasks: track finding, particle identification, and noise tagging, to enable systematic evaluation of the foundation model's adaptability. The simulated data are generated using the Pythia Monte Carlo event generator at a center of mass energy of $\sqrt{s}$ = 200 GeV and processed with Geant4 to include realistic detector conditions and signal emulation in the sPHENIX Time Projection Chamber at the Relativistic Heavy Ion Collider, located at Brookhaven National Laboratory. This dataset resource establishes a common ground for interdisciplinary research, enabling machine learning scientists and physicists alike to explore scaling behaviors, assess transferability, and accelerate progress toward foundation models in nuclear and high energy physics. The complete simulation and reconstruction chain is reproducible with the sPHENIX software stack. All data and code locations are provided under Data Accessibility.

Data Analysis, Statistics and Probability (physics↗