SMASH simulation supports two sets of potentials:
1) Skyrme and/or Symmetry potentials
2) VDF (vector density functional) model potentials, https://arxiv.org/pdf/2011.06635.pdf
Coulomb potentials can be enabled additionally.

Skyrme and VDF potentials both describe the behavior of symmetric nuclear matter. The symmetry potential can adjust the Skyrme potential (but not the VDF potential) to include effects due to isospin. The Skyrme and Symmetry potentials are semi-relativistic, while the VDF potential is fully relativistic.

Example: Configuring Skyrme Potentials
The following extract from the configuration file configures SMASH such that the Skyrme as well as the Symmetry potential are activated for the simulation. There is however no requirement to include both simultaneously. They can be switched on and off individually.

Potentials:
    Skyrme:
        Skyrme_A: -209.2
        Skyrme_B: 156.4
        Skyrme_Tau: 1.35
    Symmetry:
        S_Pot: 18.0
    Coulomb:
        R_Cut: 5.0

Note that the Coulomb potential requires a lattice while for the other potentials it can be used as an optimisation.

Example: Configuring VDF Potentials
The following extracts from the configuration file configure SMASH such that the VDF potential is activated for the simulation. In the first example, VDF potentials are configured to reproduce the default SMASH Skyrme potentials (without the symmetry potential, as it is not described within the VDF model):

Potentials:
    VDF:
        Sat_rhoB: 0.168
        Powers: [2.0, 2.35]
        Coeffs: [-209.2, 156.5]

In the second example, VDF potentials are configured to describe nuclear matter with saturation density of \(\rho_0 =\) 0.160 fm \(^{-3}\), binding energy of \(B_0 = -16.3\) MeV, the critical point of the ordinary nuclear liquid-gas phase transition at \(T_c^{(N)} = 18\) MeV and \(\rho_c^{(N)} = 0.375 \rho_0\), the critical point of the conjectured "QGP-like" phase transition at \(T_c^{(Q)} = 100\) MeV and \(\rho_c^{(Q)} = 3.0\rho_0\), and the boundaries of the spinodal region of the "QGP-like" phase transition at \(\eta_L = 2.50 \rho_0\) and \(\eta_R = 3.315 \rho_0\):

Potentials:
    VDF:
        Sat_rhoB: 0.160
        Powers: [1.7681391, 3.5293515, 5.4352788, 6.3809822]
        Coeffs: [-8.450948e+01, 3.843139e+01, -7.958557e+00, 1.552594e+00]