diff --git a/applications/solvers/combustion/plasmaReactingFoam/createFields.H b/applications/solvers/combustion/plasmaReactingFoam/createFields.H
index 06cd1407..8fe9826e 100644
--- a/applications/solvers/combustion/plasmaReactingFoam/createFields.H
+++ b/applications/solvers/combustion/plasmaReactingFoam/createFields.H
@@ -17,11 +17,6 @@ dimensionedScalar epsilon0
     physicalProperties.lookup("epsilon0")
 );
 
-dimensionedScalar k
-(
-    physicalProperties.lookup("k")
-);
-
 CSV<scalar> mueN("mueN", physicalProperties);
 CSV<scalar> DeN("DeN", physicalProperties);
 CSV<scalar> TeOfEn("TeOfEn", physicalProperties);
@@ -56,13 +51,21 @@ thermo.validate(args.executable(), "h", "e");
 volScalarField& qc = thermo.qc();
 basicMultiComponentMixture& composition = thermo.composition();
 
+//- Electron mass (default in [kg])
+const dimensionedScalar eMass = constant::atomic::me;
+
 //- Elementary charge (default in [C])
 const dimensionedScalar eCharge = constant::electromagnetic::e;
+
 //- Avogadro number (default in [1/mol])
 const dimensionedScalar NA = constant::physicoChemical::NA;
+
 //- Universal gas constant (default in [J/mol/K])
 const dimensionedScalar R = constant::physicoChemical::R;
 
+//- Pi
+const scalar pi = constant::mathematical::pi;
+
 PtrList<volScalarField>& Y = composition.Y();
 
 word inertSpecie(thermo.lookup("inertSpecie"));