diff --git a/documentation/scripts/plotting_scripts/profile_peng_qbar.py b/documentation/scripts/plotting_scripts/profile_peng_qbar.py
index 1e5f24e818..4c769050a7 100644
--- a/documentation/scripts/plotting_scripts/profile_peng_qbar.py
+++ b/documentation/scripts/plotting_scripts/profile_peng_qbar.py
@@ -6,7 +6,7 @@
from bokeh.plotting import figure, output_file, save
x = np.linspace(1.0, 3.0, 500)
-y1 = 5.0 * 1.3 * (1.0 - (1.0 / x) ** 0.6)
+y1 = 5.0 / (1.3 * (1.0 - 1.0 / x) ** 0.6)
y2 = 5.0 * (1.0 + 2.6 * (1.0 / x) ** 2.8) # Initial data for the second line
source = ColumnDataSource(data={"x": x, "y1": y1, "y2": y2})
@@ -50,7 +50,7 @@
const B = qbar.value;
const x = source.data['x'];
- const y1 = x.map(xi => A * 1.3 * (1.0 - (1.0 / xi) ** 0.6));
+ const y1 = x.map(xi => A / (1.3 * (1.0 - 1.0 / xi) ** 0.6));
const y2 = x.map(xi => B * (1.0 + 2.6 * (1.0 / xi) ** 2.8)); // Example transformation for the second line
source.data['y1'] = y1;
diff --git a/documentation/source/physics-models/plasma_current/plasma_current.md b/documentation/source/physics-models/plasma_current/plasma_current.md
index f7d31f95c9..ab384ac294 100644
--- a/documentation/source/physics-models/plasma_current/plasma_current.md
+++ b/documentation/source/physics-models/plasma_current/plasma_current.md
@@ -154,10 +154,12 @@ $$
$$
The STAR code document states that the $\bar{q}_0$ is normally equal to 3.0
-In PROCESS this is implemented differently as the function:
+In PROCESS this is implemented differently: $\bar{q}$ is obtained from the available $q_{95}$ parameter using the relation given in [^10] (Eq. 19 therein),
$$
-\bar{q} = q_{95} \times 1.3(1-\epsilon)^{0.6}
+q_{95} = 1.3\,\bar{q}\,(1-\epsilon)^{0.6}
+\quad\Rightarrow\quad
+\bar{q} = \frac{q_{95}}{1.3(1-\epsilon)^{0.6}}
$$
This is to allow the use of the available $q_{95}$ parameter as **the origin and definition of $\bar{q}_0$ is not fully known.**
@@ -205,7 +207,7 @@ The coefficient values $q_{95}$ and $\bar{q}_0$ for the PROCESS and STAR code im