diff --git a/applications/solvers/combustion/plasmaReactingFoam/createFields.H b/applications/solvers/combustion/plasmaReactingFoam/createFields.H
index ea969cd1..50a4bb1d 100644
--- a/applications/solvers/combustion/plasmaReactingFoam/createFields.H
+++ b/applications/solvers/combustion/plasmaReactingFoam/createFields.H
@@ -54,7 +54,7 @@ volScalarField& qc = thermo.qc();
 basicMultiComponentMixture& composition = thermo.composition();
 
 //- Elementary charge (default in [C])
-const dimensionedScalar& eCharge = constant::electromagnetic::e;
+const dimensionedScalar eCharge = constant::electromagnetic::e;
 
 PtrList<volScalarField>& Y = composition.Y();
 
@@ -162,3 +162,78 @@ volScalarField ne
     ),
     mesh
 );
+
+Info<< "Creating field gas number density\n" << endl;
+volScalarField ng
+(
+    IOobject
+    (
+        "ng",
+        runTime.timeName(),
+        mesh,
+        IOobject::NO_READ,
+        IOobject::NO_WRITE
+    ),
+    mesh,
+    dimensionedScalar("ng", ne.dimensions(), 1e-25)
+);
+
+Info<< "Creating field reduced electric field\n" << endl;
+volScalarField En ("En", mag(E) / (ne+ng));
+
+Info<< "Creating field electron mobility\n" << endl;
+volScalarField mue
+(
+    IOobject
+    (
+        "mue",
+        runTime.timeName(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
+    ),
+    mesh
+);
+
+Info<< "Creating field electron diffusivity\n" << endl;
+volScalarField De
+(
+    IOobject
+    (
+        "De",
+        runTime.timeName(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
+    ),
+    mesh
+);
+
+Info<< "Creating field electron temperature\n" << endl;
+volScalarField Te
+(
+    IOobject
+    (
+        "Te",
+        runTime.timeName(),
+        mesh,
+        IOobject::MUST_READ,
+        IOobject::AUTO_WRITE
+    ),
+    mesh
+);
+
+Info<< "Calculating face flux field ve\n" << endl;
+surfaceScalarField ve
+(
+    IOobject
+    (
+        "ve",
+        runTime.timeName(),
+        mesh,
+        IOobject::NO_READ,
+        IOobject::NO_WRITE
+    ),
+    linearInterpolate(mue*E) & mesh.Sf()
+);
+