Niklas Götz
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Research Subjects
Niklas Götz has a total of 15 co-authors. Here are the top co-authors by number of shared publications:
Hannah Elfner
3 shared publications- Hannah Elfner
GSI Helmholtzzentrum für Schwerionenforschung
3 shared publications Anna Schäfer
1 shared publication- Carl B. Rosenkvist
Institute for Theoretical Physics, Goethe University
1 shared publication Charles Gale
1 shared publication- Hendrik Roch
Wayne State University
1 shared publication
Latest Publications
- Open Access
- 16 January 2026
- Published by Springer Nature in The European Physical Journal C
- Vol. 86 (1) , 27
SPARKX: a software package for analyzing relativistic kinematics in collision experiments
SPARKX is an open-source Python package developed to analyze simulation data from heavy-ion collision experiments. By offering a comprehensive suite of tools, SPARKX simplifies data analysis workflows, supports multiple formats such as OSCAR2013, and integrates seamlessly with SMASH and JETSCAPE/X-SCAPE. This paper describes SPARKX’s architecture, features, and applications and demonstrates its effectiveness through detailed examples and performance benchmarks. SPARKX enhances productivity and precision in relativistic kinematics studies.
- 16 January 2026
- Published by Springer Nature in The European Physical Journal C
- Vol. 86 (1) , 27
Open AccessSPARKX: a software package for analyzing relativistic kinematics in collision experiments
SPARKX is an open-source Python package developed to analyze simulation data from heavy-ion collision experiments. By offering a comprehensive suite of tools, SPARKX simplifies data analysis workflows, supports multiple formats such as OSCAR2013, and integrates seamlessly with SMASH and JETSCAPE/X-SCAPE. This paper describes SPARKX’s architecture, features, and applications and demonstrates its effectiveness through detailed examples and performance benchmarks. SPARKX enhances productivity and precision in relativistic kinematics studies.
- Open Access
- 30 September 2025
- Published by University Library J. C. Senckenberg
Extraction of transport coefficients in a hybrid approach for heavy-ion collisions
This thesis investigates transport coefficients in relativistic heavy-ion collisions using a hybrid approach. Relativistic heavy-ion collisions recreate extreme conditions of temperature and density, enabling the study of the quark-gluon plasma and the QCD phase diagram. To describe the full dynamical evolution, the SMASH-vHLLE hybrid model is employed, combining hadronic transport for the non-equilibrium initial and late stages with viscous hydrodynamics for the dense intermediate stage. A key focus is the determination of shear and bulk viscosities, η/s and ζ/s, as functions of temperature and baryochemical potential. The sensitivity of final-state observables to different initial state models (IP-Glasma, TRENTo, SMASH) is analyzed, demonstrating that initial transverse momentum distributions strongly influence flow anisotropies. To constrain transport coefficients, Bayesian inference is applied, combining experimental data with model predictions in a statistically consistent framework. Results indicate a strong temperature dependence of shear viscosity, vanishing near the crossover temperature, and substantial bulk viscosity values. The baryon-density dependence of shear viscosity remains less constrained, highlighting the need for more experimental data at finite μB. This work represents the first Bayesian extraction of temperature- and density-dependent viscosities in hybrid models with hadronic initial states. The findings improve our understanding of QCD matter under extreme conditions and have implications for future experiments at FAIR and astrophysical studies of neutron stars.
- 30 September 2025
- Published by University Library J. C. Senckenberg
Open AccessExtraction of transport coefficients in a hybrid approach for heavy-ion collisions
This thesis investigates transport coefficients in relativistic heavy-ion collisions using a hybrid approach. Relativistic heavy-ion collisions recreate extreme conditions of temperature and density, enabling the study of the quark-gluon plasma and the QCD phase diagram. To describe the full dynamical evolution, the SMASH-vHLLE hybrid model is employed, combining hadronic transport for the non-equilibrium initial and late stages with viscous hydrodynamics for the dense intermediate stage. A key focus is the determination of shear and bulk viscosities, η/s and ζ/s, as functions of temperature and baryochemical potential. The sensitivity of final-state observables to different initial state models (IP-Glasma, TRENTo, SMASH) is analyzed, demonstrating that initial transverse momentum distributions strongly influence flow anisotropies. To constrain transport coefficients, Bayesian inference is applied, combining experimental data with model predictions in a statistically consistent framework. Results indicate a strong temperature dependence of shear viscosity, vanishing near the crossover temperature, and substantial bulk viscosity values. The baryon-density dependence of shear viscosity remains less constrained, highlighting the need for more experimental data at finite μB. This work represents the first Bayesian extraction of temperature- and density-dependent viscosities in hybrid models with hadronic initial states. The findings improve our understanding of QCD matter under extreme conditions and have implications for future experiments at FAIR and astrophysical studies of neutron stars.
- 29 July 2025
- Published by American Physical Society (APS) in Physical Review C
- Vol. 112 (1)
Bayesian analysis of a hybrid approach with initial conditions from hadronic transport- 29 July 2025
- Published by American Physical Society (APS) in Physical Review C
- Vol. 112 (1)
Bayesian analysis of a hybrid approach with initial conditions from hadronic transport- 3 September 2024
- Published by American Physical Society (APS) in Physical Review C
- Vol. 110 (3) , 034901
Role of initial transverse momentum in a hybrid approach
Background: Significant theoretical uncertainties exist with respect to the initial condition of the hydrodynamic description of ultrarelativistic heavy-ion collisions. Several approaches exist, of which some contain initial momentum information. Its impact is commonly assumed to be small and final flow is seen as a linear response to the initial-state eccentricity.
- 3 September 2024
- Published by American Physical Society (APS) in Physical Review C
- Vol. 110 (3) , 034901
Role of initial transverse momentum in a hybrid approach
Background: Significant theoretical uncertainties exist with respect to the initial condition of the hydrodynamic description of ultrarelativistic heavy-ion collisions. Several approaches exist, of which some contain initial momentum information. Its impact is commonly assumed to be small and final flow is seen as a linear response to the initial-state eccentricity.
- Open Access
- 3 April 2024
- Published by American Physical Society (APS) in Physical Review C
- Vol. 109 (4) , 049901
Erratum: Out-of-equilibrium photon production in the late stages of relativistic heavy-ion collisions [Phys. Rev. C 105 , 044910 (2022)]
DOI: https://doi.org/10.1103/PhysRevC.109.049901
- 3 April 2024
- Published by American Physical Society (APS) in Physical Review C
- Vol. 109 (4) , 049901
Open AccessErratum: Out-of-equilibrium photon production in the late stages of relativistic heavy-ion collisions [Phys. Rev. C 105 , 044910 (2022)]
DOI: https://doi.org/10.1103/PhysRevC.109.049901
- Open Access
- 14 March 2024
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2024 (3) , 1-29
Accelerating HEP simulations with Neural Importance Sampling
Many high-energy-physics (HEP) simulations for the LHC rely on Monte Carlo using importance sampling by means of the VEGAS algorithm. However, complex high-precision calculations have become a challenge for the standard toolbox, as this approach suffers from poor performance in complex cases. As a result, there has been keen interest in HEP for modern machine learning to power adaptive sampling. While previous studies have shown the potential of normalizing-flow-powered neural importance sampling (NIS) over VEGAS, there remains a gap in accessible tools tailored for non-experts. In response, we introduce ZüNIS, a fully automated NIS library designed to bridge this divide, while at the same time providing the infrastructure to customise the algorithm for dealing with challenging tasks. After a general introduction on NIS, we first show how to extend the original formulation of NIS to reuse samples over multiple gradient steps while guaranteeing a stable training, yielding a significant improvement for slow functions. Next, we introduce the structure of the library, which can be used by non-experts with minimal effort and is extensivly documented, which is crucial to become a mature tool for the wider HEP public. We present systematic benchmark results on both toy and physics examples, and stress the benefit of providing different survey strategies, which allows higher performance in challenging cases. We show that ZüNIS shows high performance on a range of problems with limited fine-tuning.
- 14 March 2024
- Published by Springer Nature in Journal of High Energy Physics
- Vol. 2024 (3) , 1-29
Open AccessAccelerating HEP simulations with Neural Importance Sampling
Many high-energy-physics (HEP) simulations for the LHC rely on Monte Carlo using importance sampling by means of the VEGAS algorithm. However, complex high-precision calculations have become a challenge for the standard toolbox, as this approach suffers from poor performance in complex cases. As a result, there has been keen interest in HEP for modern machine learning to power adaptive sampling. While previous studies have shown the potential of normalizing-flow-powered neural importance sampling (NIS) over VEGAS, there remains a gap in accessible tools tailored for non-experts. In response, we introduce ZüNIS, a fully automated NIS library designed to bridge this divide, while at the same time providing the infrastructure to customise the algorithm for dealing with challenging tasks. After a general introduction on NIS, we first show how to extend the original formulation of NIS to reuse samples over multiple gradient steps while guaranteeing a stable training, yielding a significant improvement for slow functions. Next, we introduce the structure of the library, which can be used by non-experts with minimal effort and is extensivly documented, which is crucial to become a mature tool for the wider HEP public. We present systematic benchmark results on both toy and physics examples, and stress the benefit of providing different survey strategies, which allows higher performance in challenging cases. We show that ZüNIS shows high performance on a range of problems with limited fine-tuning.
- Open Access
Accelerating HEP simulations with Neural Importance Sampling
Many high-energy-physics (HEP) simulations for the LHC rely on Monte Carlo using importance sampling by means of the VEGAS algorithm. However, complex high-precision calculations have become a challenge for the standard toolbox, as this approach suffers from poor performance in complex cases. As a result, there has been keen interest in HEP for modern machine learning to power adaptive sampling. While previous studies have shown the potential of normalizing-flow-powered neural importance sampling (NIS) over VEGAS, there remains a gap in accessible tools tailored for non-experts. In response, we introduce ZüNIS, a fully automated NIS library designed to bridge this divide, while at the same time providing the infrastructure to customise the algorithm for dealing with challenging tasks. After a general introduction on NIS, we first show how to extend the original formulation of NIS to reuse samples over multiple gradient steps while guaranteeing a stable training, yielding a significant improvement for slow functions. Next, we introduce the structure of the library, which can be used by non-experts with minimal effort and is extensivly documented, which is crucial to become a mature tool for the wider HEP public. We present systematic benchmark results on both toy and physics examples, and stress the benefit of providing different survey strategies, which allows higher performance in challenging cases. We show that ZüNIS shows high performance on a range of problems with limited fine-tuning.
Open AccessAccelerating HEP simulations with Neural Importance Sampling
Many high-energy-physics (HEP) simulations for the LHC rely on Monte Carlo using importance sampling by means of the VEGAS algorithm. However, complex high-precision calculations have become a challenge for the standard toolbox, as this approach suffers from poor performance in complex cases. As a result, there has been keen interest in HEP for modern machine learning to power adaptive sampling. While previous studies have shown the potential of normalizing-flow-powered neural importance sampling (NIS) over VEGAS, there remains a gap in accessible tools tailored for non-experts. In response, we introduce ZüNIS, a fully automated NIS library designed to bridge this divide, while at the same time providing the infrastructure to customise the algorithm for dealing with challenging tasks. After a general introduction on NIS, we first show how to extend the original formulation of NIS to reuse samples over multiple gradient steps while guaranteeing a stable training, yielding a significant improvement for slow functions. Next, we introduce the structure of the library, which can be used by non-experts with minimal effort and is extensivly documented, which is crucial to become a mature tool for the wider HEP public. We present systematic benchmark results on both toy and physics examples, and stress the benefit of providing different survey strategies, which allows higher performance in challenging cases. We show that ZüNIS shows high performance on a range of problems with limited fine-tuning.
- Open Access
Momentum anisotropy generation in a hybrid approach
Anisotropic flow emerges in all three of hybrid approaches: initial conditions, viscous relativistic hydrodynamics as well as hadronic transport. Previous works focus mainly on a constant or temperature dependent shear viscosity $η/s$. Here instead, we study qualitatively the effect of a generalized $η/s(T,μ_B)$ in the hybrid approach SMASH-vHLLE-hybrid. The parameterization takes into account the constraints of matching to the transport coefficients in the hadronic phase, as well as of recent Bayesian analysis results. We compare the effect of the different parameterizations in the intermediate energy region of $\sqrt{s_{NN}}$=7.7 - 39.0 GeV. We observe that using the energy density dependent parameterization decreases the effect of the point of particlization. In addition, we quantify the uncertainty due to different initial state profiles, including the SMASH initial conditions as well as TrENTo and IP-Glasma profiles. It can be shown that the initial state transverse momentum impacts final state momentum anisotropy.Open AccessMomentum anisotropy generation in a hybrid approach
Anisotropic flow emerges in all three of hybrid approaches: initial conditions, viscous relativistic hydrodynamics as well as hadronic transport. Previous works focus mainly on a constant or temperature dependent shear viscosity $η/s$. Here instead, we study qualitatively the effect of a generalized $η/s(T,μ_B)$ in the hybrid approach SMASH-vHLLE-hybrid. The parameterization takes into account the constraints of matching to the transport coefficients in the hadronic phase, as well as of recent Bayesian analysis results. We compare the effect of the different parameterizations in the intermediate energy region of $\sqrt{s_{NN}}$=7.7 - 39.0 GeV. We observe that using the energy density dependent parameterization decreases the effect of the point of particlization. In addition, we quantify the uncertainty due to different initial state profiles, including the SMASH initial conditions as well as TrENTo and IP-Glasma profiles. It can be shown that the initial state transverse momentum impacts final state momentum anisotropy. - Open Access
- 10 November 2022
- Published by American Physical Society (APS) in Physical Review C
- Vol. 106 (5) , 054904
Temperature and net baryochemical potential dependence of in a hybrid approachBackground:The effect of a predicted nonconstant ) is largely unexplored in hydrodynamic simulations. Previous studies focus only on a temperature dependence or even only a constant effective shear viscosity.- 10 November 2022
- Published by American Physical Society (APS) in Physical Review C
- Vol. 106 (5) , 054904
Open AccessTemperature and net baryochemical potential dependence of in a hybrid approachBackground:The effect of a predicted nonconstant ) is largely unexplored in hydrodynamic simulations. Previous studies focus only on a temperature dependence or even only a constant effective shear viscosity. - Open AccessTemperature and net baryochemical potential dependence of $η/s$ in a hybrid approachIn this work, the qualitative impact of the net baryochemical potential dependence of the shear viscosity to entropy density ratio $\eta/s$ in hydrodynamical simulations is studied. The effect of a predicted non-constant $\eta/s$($\mu_B$) is largely unexplored in hydrodynamic simulations. Previous studies focus only on a temperature dependence or even only a constant effective shear viscosity. This work addresses this issue by studying qualitatively the effect of a generalized $\eta/s$($T,\mu_B$) in the hybrid approach SMASH-vHLLE, composed of the hadronic transport approach SMASH and the (3+1)d viscous hydrodynamic code vHLLE. In order to reduce the bias of the result on the equation of state used in the hydrodynamic part of the model, $\eta/s$ is parameterized directly in the energy density and net baryon number density. The parameterization takes into account the constraints of matching to the transport coefficients in the hadronic phase, as well as pQCD results. This work compares the impact of the density dependence for different system sizes and energies and compares the observables with experimental results in the RHIC - BES region $\sqrt{s_{NN}}$ =7.7 - 39.0 GeV, as the effect of this generalization is especially relevant for intermediate collision energies, for which the system is in equilibrium for a relevant amount of time, but the net baryochemical potential does not vanish. It is shown that the effect of an explicit net baryon number dependence on the elliptic flow is negligible and only relevant in the early stages of the collision. Additionally, we find that the proposed parameterization could be a good proxy for the shear viscosity in the non-equilibrium hadronic transport stage.Open AccessTemperature and net baryochemical potential dependence of $η/s$ in a hybrid approachIn this work, the qualitative impact of the net baryochemical potential dependence of the shear viscosity to entropy density ratio $\eta/s$ in hydrodynamical simulations is studied. The effect of a predicted non-constant $\eta/s$($\mu_B$) is largely unexplored in hydrodynamic simulations. Previous studies focus only on a temperature dependence or even only a constant effective shear viscosity. This work addresses this issue by studying qualitatively the effect of a generalized $\eta/s$($T,\mu_B$) in the hybrid approach SMASH-vHLLE, composed of the hadronic transport approach SMASH and the (3+1)d viscous hydrodynamic code vHLLE. In order to reduce the bias of the result on the equation of state used in the hydrodynamic part of the model, $\eta/s$ is parameterized directly in the energy density and net baryon number density. The parameterization takes into account the constraints of matching to the transport coefficients in the hadronic phase, as well as pQCD results. This work compares the impact of the density dependence for different system sizes and energies and compares the observables with experimental results in the RHIC - BES region $\sqrt{s_{NN}}$ =7.7 - 39.0 GeV, as the effect of this generalization is especially relevant for intermediate collision energies, for which the system is in equilibrium for a relevant amount of time, but the net baryochemical potential does not vanish. It is shown that the effect of an explicit net baryon number dependence on the elliptic flow is negligible and only relevant in the early stages of the collision. Additionally, we find that the proposed parameterization could be a good proxy for the shear viscosity in the non-equilibrium hadronic transport stage.